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ABSTRACT 


This report summarizes the results of studies performed under NASA 
Grant NSG-5011 and entitled Satellite Support to Weather Modification 
in the Western U. S. Region. These studies were completed during the 
period 1 April, 1974 to the present. The applications of meteorological 
satellite data to both summer and winter weather modification programs 
were addressed in this research. These studies included the appraisal 
of the capability of satellites to assess seedability, to provide real- 
time operational support, and to assist in the post-experiment analysis 
of a seeding experiment. This research has led to the incorporation of 
satellite observing systems as a major component in the Bureau of 
Reclamations weather modification activities. 

Satellite observations acquired under this grant are an integral 
part of the South Park Area Cumulus Experiment (SPACE) . The objective 
of the SPACE research program is to formulate a quantitative hypothesis 
for enhancing precipitation from orographically induced summertime 
mesoscale convective systems (orogenic mesoscale systems). Satellite 
observations assist in classifying the important mesoscale systems, 
defining their frequency and coverage, defining a potential area of 
effect. Satellite studies of severe storms were investigated under 
this grant. This document addresses the progress of those studies. 
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FINAL REPORT: Studies of Satellite Support to Weather Modification 

in the Western U. S. Region 
NASA/GLAS/ Severe Storms Research Program 
Grant NSG-5011 

1.0 INTRODUCTION 

This report summarizes the results of studies of satellite support 
to weather modification in the Western U. S, region for the period 
1 April 1974 to the present. The early stages of this research were 
focused on the collection and analysis of satellite data for four major 
weather modification programs : 

1) Meso/Synoptic Storm Modifications (National Oceanic and 
Atmospheric Administration) 

2) National Hail Research Experiment (NCAR/NSF, RANN) 

3) High Plains Precipitation Enhancement (HIPLEX; U. S. Dept, 
of Interior, Bureau of Reclamation) 

4) Wintertime Orographic Cloud Modification (U. S. Dept, of 
Interior, CSU/Extra-Area Effect (NSF/RANN)). 

During the period of this research, the first program never 
materialized, and the second program was rapidly phased out, while 
the HIPLEX program and wintertime orographic cloud modification studies 
have become increasingly functional. The Bureau of Reclamations 
weather modification programs have become increasingly dependent upon 
satellite support to the extent that they have provided funding for 
investigations of satellite support to both their summer and winter 
programs under Contract 6-07- DR -20020. The results of this research 
are described in section 2.0. The CSU/Extra Area Effect Experiment 
has relied on satellite data for important input. 

Colorado State University has operated a series of field experiments 
in the Colorado Mountains (the South Park Area Cumulus Experiment - 
(SPACE) which are designed to investigate the ^relationship between 
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mountain-generated cumulus convection and precipitating mesoscale 
systems over the High Plains. This experiment provided an excellent 
opportunity for testing satellite analysis techniques because of the 
extensive ground truth data collected during experiments. Therefore 
the later stages of the research covered by this grant emphasized the 
SPACE experiment. The ultimate objective of this research is to 
formulate a quantitative hypothesis for the modification of precipitating 
mesoscale systems over the High Plains via modification of orogenic 
cumuli. The role of satellite meteorology in these investigations is 
described in section 3.0. 

The intimate relationship between severe storm research and 
convective cloud modification is also pointed out in section 3.0. In 
section 4.0 recent progress in severe storm analysis using satellite- 
derived data is discussed. 

2.0 SUMMARY OF SCIENTIFIC PROGRESS DURING THE PERIOD 1 APRIL 1974 

TO 14 NOVEMBER 1976 . 

A thorough summary of the results of this research is given by 
Reynolds, VonderHaar and Grant (1978; Appendix Al) . This paper surveys the 
applications of meteorological satellite data to both summer and winter 
weather modification programs. Included in the discussion is the 
appraisal of the capability of satellites to assess the potential for 
precipitation enhancement (seedability) , to provide real-time operational 
support, and to assist in the post-experiment analysis of a seeding 
experiment. This paper demonstrates that present satellite observing 
systems have become an essential component in weather modification 
programs. Moreover, anticipated future advances in satellite observing 
systems will make the weather modification community increasingly dependent 
upon satellite technology. 
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3.0 THE SOUTH PARK AREA CUMULUS EXPERIMENT (SPACE) 

The SPACE has been an evolving program whose ultimate objective 
is to develop a capability of precipitation enhancement from summertime 

r 

convection over the Colorado Mountains and High Plains. Since the bulk 

of precipitation from summertime convection occurs in organized mesoscale 

*■ 

systems, the experiment has naturally evolved into investigations into 
the factors contributing to the organization and intensity of such 
systems . 

Initial satellite studies concentrated on the determination of satellite- 
derived cloud climatology (Stodt and Grant (1976); A2, Stodt (1978)) the 
determination of satellite 'ground truth' regarding cloud brightness 
(Reynolds, McKee and Danielson (1976, 1978); A3, A4) and the factors 
contributing to initial convection (Breed (1975); A5) . Cotton, George 
and Pielke (1976; A6) demonstrated the analogy between sea-breeze 
generated mesoscale convective systems and mountain-generated or orogenic 
mesoscale convective systems. 

SPACE-1977 was a major experiment which utilized a host of observation 
systems including triple-Doppler radar, convectional radar, cloud physics 
and turbulence equipped aircraft, Lidar, rawinsondes, micromet towers, 
boundary layer profilers, an acoustic sounder, recording mesomet stations, 
fhe NCAR PAM system, stereo photography and satellite observations. 

A complete description of the SPACE-77 is given m SPACE LOG (Danielson 
and Cotton (1977)). The SPACE-77 data set provides comprehensive 
observations of a' wide variety of cumulus-scale and convective mesoscale 
meteorological events ranging from locally suppressed cumulus convection, 
heavy precipitating and hail-producing cumulonimbus systems exhibiting 
little motion relative to the mountain source region, and heavy 
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precipitating, tornado-producing, eastward propagating mesoscale 
systems . 

July 19, 1977 has been selected as the first day for comprehensive 
case study analysis. George and Cotton (1978; A7) have described the 
characteristics of the evolving mesoscale system on that day. The day 
was characterized by locally heavily precipating cumulonimbi and several 
hail falls. Since the clouds exhibited little motion relative too the 
observational system, this day provides an excellent opportunity to 
test out the analysis techniques and to perform three-dimensional model 
simulations of storm-scale and mesoscale events. Knupp, Danielson 
and Cotton (1978; A8) have described the morphological characteristics 
of the radar observations of several cumulonimbus cells that developed 
on that day. The analyses performed on this day will provide a number 
of opportunities for ground truth-interpretation of satellite-derived 
data. This is especially enhanced by the use of 3D model simulations 
as a part of the overall analysis. Because of the relative spatial 
stationarity of the mesoscale systems that formed on July 19, it provides 
a contrasting observational model to the days characterized by eastward 
propagating systems. 

August 4, August 8 and August 10, 1977 represent three days in 
which convection first occurred over the mountains and then propagated 
eastward as a major mesoscale convective system. These are not the only 
such events that occurred during the observational period but are the most 
active in terms of precipitation, and severe weather events. August 4 
has been selected as a day for comprehensive case study analysis. It 
is intended that the analysis of this day should serve as a descriptive 
model of eastward propagating orogenic mesoscale systems . The detailed 
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observations over South Park are complemented by 3 minute-interval 
satellite data on this particular day. Figure 3.1 illustrates the 
scenario of events observed by satellite on August 4, 1977. Deep 
convection first formed over the mountains. By 02Z (August 5) a line 
of very deep cells extended along the entire Colorado eastern border 

i 

into the Texas panhandle with the most intense cells east of South Park. 
By 18Z the entire complex had migrated well into western Missouri. This 
system produced a tornado near Kiowa, Colorado and produced in excess 
of 50mm of precipitation over widespread areas of northeast and Central 
Kansas. That amount of precipitation represents roughly 15% of the total 
annual precipitation in that region. 

While a thorough climatic study has not been performed, it appears 
that eastward propagating orogenic mesoscale systems contribute about 
75% of the annual precipitation over eastern Colorado and the Kansas 
High Plains. Thus, a major portion of the weather modification potential 
over the High Plains must be tied to the modification of such systems. 
Because of the scale of these systems (a north-south dimension of several 
hundred kilometers) and the total extent of propagation (greater than 
1500 km) , satellites are an excellent platform for defining the charac- 
teristics of such systems and the potential ’’area of effect" of any 
implemented mesoscale modification hypotheses. As can be seen in 
Figure 3.1, the potential area of effect of a modification capability 
is a large one. 

Another factor of serious concern is that the same systems that are 
dominant contributors to precipitation are also sporadic generators of 
tornadoes and hail. These systems do not produce major outbreaks of 
severe weather but instead often a few small tornadoes or moderate hail 
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Figure 3.1. GOES-1 visible and infrared satellite imagery on August 4, 
1977 for the period 1800 GMT to 1100 GMT (August 5, 1977). 
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Figure 3.1 Continued. 
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events such as the Kiowa tornado. Thus any methodology for enhancing 
precipitation from nesoscale convective systems must also induce the 
capability for recognizing severe weather potential ana ae fining 
immediate response strategies. Thus satellite support is essential 
to the definition of local severe weatner indices. In aadition, the 
implementation and testing of a convective mesoscale modification 
hypothesis will require active satellite support for the routine 
evaluations of these severe weather indices ana for the activation of 
such response strategies as immediate shutdown of seeding operations, 
avoidance of a potential severe weather zone or implementation of alternate 
modification techniques. 

Thus satellite-support is essential to the development of a 
quantitative mesoscale modification hypothesis and to the actual 
implementation and testing of such a hypothesis. Moreover one can see 
that the formulation and implementation of such a hypothesis is intimately 
tied to severe storm research. 

4.0 SEVERE STORMS RESEARCH IN TEE CENTRAL PLAINS USING SATELLITE DATA 


Funding from NSG-501L also supported a variety of severe convective 
storm research efforts that utilized satellite data. The particular 
day studied was April 24, 1975; a day when a number of severe storms 
occurred over the Central Plains and Central Mississippi Valley 
(see Fig. 4.1). A destructive tornado, spasmed by an isolated 
supercell thunderstorm, struck Neosho, Missouri, just before sunset. 

An unusually comprehensive set of conventional and special satellite 
data were available for this particular day. 

An upper-air sounding taken at Monett, Missouri, (approximately 
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25 nm east-southeast of Neosho) less than 4 h prior to the tornado 
indicated that significant lifting and/or air mass modification would 
be required to initiate deep convection. This sounding is shown in 
Fig. 4.2. Note the deep layer of negative buoyancy between the 
LCL (lifted condensation level) and LFC (level of free convection) . 
Wilson (1976) computed large scale vertical motion fields using 2100 
GMT soundings and his 700 mb oj field is shown in Fig. 4.3 with 
regions in which storms developed during the following several hours 
indicated. Since the storms developed in regions of weak upward 
motion, or in areas of subsidence, mesoscale features and forcing 
mechanisms must have played a dominant role in triggering storm 
development. Indeed, during late afternoon and evening when the 
storms x?ere most intense they were relatively unorganized (e.g. no 
distinct squall or line orientation) and isolated. 

Surface analyses (Fig. 4.4) for 18, 21, 00, and 03 GMT indicate 

that a sub-synoptic scale low (see Tegtmeir, 1974) and an attendant 
narrow zone of dry air moved eastward along the Kansas /Oklahoma border. 

An intense supercell storm developed just ahead of the dryline and moved 
eastward — becoming tornadic when it interacted with an old thunder- 
storm produced thermal boundary in extreme northeast Oklahoma and 
southwest Missouri. 

Analyses of the 500 mb level (Fig. 4.5) indicate that the meso- 
scale surface low and the intense afternoon and evening severe storms 
were likely associated with a weak short-wave trough. This feature was 
moving rapidly eastward and weakening as a stronger short-wave had taken 
over as the dominant feature and helped to organize the storm activity 
into a well defined squall-line. 
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Figure 4.2. Monett, Mo. upper-air sounding for 24 April 1975, 
2100 GMT. 



Figure 4.3. 700 mb w field (u bars s , from Wilson, 1977) for 

2100 GMT, 24 April 1975, regions where severe thunder 
storms developed are cross-hatched. 



Figure 4.4a. Surface Analysis for 24 April 1975, 1800 GMT. Winds 

are in knots with full barb equal to 10 kt. Temperatures 
are in °F and 10 = 1010 mb. Surface dryline is shown 
with open frontal barbs. Thunderstorm outflow boundary 
is shown as squall line with frontal barbs added to 
indicate direction of movement. 







Figure 4.4b. Same as 4.4a except for 2100 GMT. 
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Figure 4.4c. Same as 4.4a except for 25 April 1975, '0000 GMT. 







Figure 4.4d. Same as 4.4a. except for 25 April 1975, 0300 G^fT. 
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Figure 4.5a. 500 mb analysis for 24 April 1975, 2100 GMT. Short- 

wave troughs are indicated by heavy dashed lines and 
winds are m kt (flag = 50 kt.). 
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Figure 4.5b. Same as 4.5a except: for 25 April 1975, 0600 GMT. 
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Many interesting and significant interactions between features 

on different scales appeared to occur on this particular day. Some 

of the more important effects seem to be upscale energy inputs from 

storm and mesoscale convective features that acted to alter the 

larger scale jetstream level feature and magnitudes. Jetstream charts 

~ 1 

(maximum \dnds aloft are plotted m ms at each station) are shown in 
Fig. 4.6. Satellite imagery indicated thunderstorm areas are also 
shown. At 25/0000 GMT (Fig. i.6a) several severe thunderstorms had 
developed in the diffluence zone between' the strong subtropical 
jetstream (y 130 kt maximum) and a moderate polar jet (y 110 kt 
maximum) . Only 6 h later dramatic changes were evident m the large 
scale flow. The subtropical jet had weakened to about 110 kt while 
the polar jet had reformed northward over the plains and undergone 
remarkable intensification over the lower Great Lakes region. A jet; 
maximum of about 140 kt was now present,. The concurrent satellite 
photo indicated an extensive canopy of thunderstorm generated cirrus spanned 
the region between the two jetstream maxima (note that radar data indicated 
that the active storms at 0600 GMT were oriented along the southern- 
third of the cirrus region) . These changes were undoubtedly due m 
part to the approaching second short-wave trough, hut it certainly 
appears that convective scale feedbacks also played a role in the 
modification of the large scale jet structure. 

The satellite data used in the various studies undertaken were 
in digital format on magnetic tape. A number of cloud growth rate 
computations and radar/satellite cloud top height comparisons were 
made. Low-level winds were generated by tracking clouds on the NASA 
AOIPS system. Dynamic parameters were calculated for these wind" 
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Figure 4-6a. 


Maximum winds aloft for 25 April 1975, 0000 GMT. 
Winds are m ms~l with flag - 50 ms“l. Satellite 
indicated thunderstorm regions are hatched. 




21 



Figure 4.6b. Same as 4.6a except for 25 April 1975, 0600 G^IT. 
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fields- The most significant results indicated that satellite -jind 
fields might be combined with a moisture analysis (either conven- 
tionally or remotely sensed) to enable computation of moisture 
divergence fields- Several such fields were generated, and it was 
shown that regions characterized by high values of moisture conver- 
gence at 2100 GMT were generally the areas of subsequent severe storm 
development. Real time computation of such fields could provide 
valuable input to forecast and warning systems. 

Papers detailing the specific research completed under NSG-5011 
are included as appendices to this report. 
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5.0 CONCLUDING REMARKS 

This research has elucidated a variety of satellite applications 
to weather modification programs in the Western U. S. Region. Several 
of these satellite applications have become a routine part of the decision 
making process and analysis procedures of the Bureau of Reclamation's 
weather modification programs. 

In addition, satellite support has been an integral part of the 
CSU/SPACE program; a program operating on the frontiers of weather 
modification research. Satellite observations have also become an 
integral part of severe storm analysis. 
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Modification 
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and Lewis O. Grant 
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Colorado State University 
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Abstract 

During the past several years, many weather modification 
programs have been incorporating meteorological satellite 
data into both the operations and the analysis phase of these 
projects. This has occurred because of the advancement of 
the satellite as a mesoscale measurement platform, both 
temporally and spatially, and as the availability of high 
quality data has increased. This paper surveys the applica- 
tions of meteorological satellite data to both summer and 
winter weather modification programs. A description of the 
types of observations needed by the programs is given, and 
an assessment of how accurately satellites can determine these 
necessary parameters is made. 

1. Introduction 

The advancement in meteorological satellite technology 
has been rapid over the past decade. In the 1960s we 
saw the development of sun synchronous polar orbiter 
satellites with both visible and IR wavelength sensors 
that could monitor “weather” twice daily. These were 
very useful in monitoring synoptic scale disturbances but 
lacked the temporal frequency necessary to monitor 
short-term weather features. In 1966 the first geosyn- 
chronous meteorological satellite was launched; it pro- 
vided visible imagery (4 km resolution at satellite sub- 
point) twice hourly and opened completely new avenues 
for satellite applications (Suomi and Vonder Haar, 
1969). The time domain allowed such measurements as 
wind from cloud motion, cloud growth rates, identifica- 
tion of mesoscale disturbances causing short-term weather 
phenomena, and many others. In the 1970s the NOAA 
and DMSP (Defense Meteorological Satellite Program) 
series of polar orbiter satellites were developed, carrying 
sophisticated visible and IR radiometers and achieving 
very high resolution visible and IR imagery (1 km with 
the NOAA VHRR (Very High Resolution Radiometer) 
and 0.6 km with the DMSP VHR (very high resolution) 
sensor). The satellites also carried the first vertical tem- 
perature sounders, which allowed remote sensing of the 
temperature structure of the environment. NASA’s Nim- 
bus satellites also provided sounding capabilities includ- 
ing the microwave region and increased both the vertical 
and the horizontal resolution of the temperature and 
moisture measurements. The second generation geosyn- 
chronous satellites called the SMS-GOES (Synchronous 
Meteorological Satellite-Geosynchronous Operational 
Environmental Satellite) (Fordyce et al., 1974) have been 
developed and are providing both visible (1 km resolu- 
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tion) and IR (8 km resolution) data twice hourly. The 
IR capability allows 24 h a day coverage of weather 
features for the first time. There are two GOES satellites 
presently in operation, GOES-E (east) positioned at 
75°\V and GOES-W (west) positioned at 135°W. Dur- 
ing the First GARP Global Experiment (FGGE), the 
world will be girdled by such satellites (two from the 
United States, one from Europe, and one from Japan). 

With the ATS satellites came a data dissemination 
network called the Satellite Field Service Station (SFSS). 
These were expanded and received more emphasis once 
the SMS-GOES satellites were launched. There are six 
of these SFSSs, which act as hubs or distribution points 
for photographic data users. Thus, near real-time im- 
agery (20 min after scan) is now available from the GOES 
satellites to anyone tied into this distribution network. 
The quantitative digital data from the GOES satellites 
can be even more valuable both in research analysis and 
in operations if a suitable system is available to process 
the data. This paper will discuss how the digital satellite 
data can be used for weather modification applications 
and will describe some interactive time domain image- 
processing systems that have been and are being de- 
veloped for applications use. 

Dennis et al. (1973) investigated the usefulness of 
meteorological satellites in weather modification pro- 
grams. They concluded that, with the observational re- 
quirements needed in a weather modification program, 
presently available satellite information (as of 1973) 
could not by itself satisfy these requirements. Their 
strongest recommendation was that imagery from a geo- 
synchronous satellite with improved spatial resolution 
and the capability of accurately measuring cloud top 
temperature (CTT) be provided in real time to the field 
users. As we have stated, the SMS-GOES satellite sys- 
tems have satisfied these major requirements and have 
provided even more in the way of support. We will dis- 
cuss these present capabilities, as well as point out what 
satellites may offer in the next decade. 

2. Application of satellite information to 
weather modification programs 

a. Summertime cumulus modification experiments 

1) Cloud climatologies 

One of the first objectives in beginning a summer cumu- 
lus weather modification program in a given area is 
determining whether enough convective clouds are 
present naturally for seeding to take place and what 
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their natural characteristics are, i.e., isolated cells, well- 
developed squall lines, etc. The next objective is to 
determine whether these cumulus clouds are developing 
with the proper microphysical characteristics so that 
seeding would change their natural growth pattern in 
such a way as to increase precipitation. 

There are several ways in which one could try to de- 
termine the first objective. One would be to take a long 
history of surface cloud observations to determine how 
many days convective clouds were observed and the de- 
gree of intensity of the day’s activity. This is appropriate 
when a limited area is to be chosen or where a number 
of surface weather stations are available in a given area 
to allow complete coverage. This lias been attempted by 
the authors and others and constitutes a very lengthy 
undertaking. Limited visibility due to rain, low clouds, 
etc., will severely limit the quality of these types of ob- 
servations. Also, surface observations could not accurately 
provide numbers and sizes of convective clouds observed. 
Another procedure would be to use radar plan position 
indicator (PPI) data to determine the number of pre- 
cipitating convective cells, sizes, intensity, rate of growth, 
and lifetimes (Bark, 1975). Unfortunately, for many of 
the areas that are being considered for summer weather 
modification programs, no radar data are available, or 
the data have not been archived in a format that could be 
easily processed. Use of radar data is, by its very nature, 
limited to precipitating clouds. Satellite visible data can 
sense all convective type clouds above a certain size 
threshold both before and after they reach the precipita- 
tion stage or even if they never precipitate. A compara- 
tive study between radar and satellite climatology per- 
formed simultaneously may show the ratio of the number 
of clouds observed (satellite) versus the number of pre- 
cipitating clouds, which may indicate the potential of 
a locality for cloud seeding. 

The third alternative in determining the seeding po- 
tential of a given area is through the use of satellite data. 
The most useful satellite data in such a study would be 
geosynchronous satellite data, preferably with both 
visible and IR measuring capabilities. This type of in- 
formation has only been available since 1975 when 
SMS-1 was first positioned at 75°W with the VISSR 
(Visible and Spin Scan Radiometer) on board, providing 
high-resolution (1 km visible, 8 km IR) data on a half- 
hour basis. This type of data is the best source of in- 
formation since the geosynchronous data provide the 
temporal frequency of the convective activity. The high- 
resolution visible data allow detection and measurement 
of very small cumulus clouds, and the IR data allow 
measurement of CTTs for the larger size clouds. The 
latter is needed to determine if the cloud has reached a 
temperature level where supercooled water may be 
present and seeding can be preformed. 1 Work is now 
beginning in the archiving of this information in digital 
form so that a true cloud climatology can be developed 

i This technique would obviously not be applicable to 
tropical cumulus undergoing the warm rain precipitation 
process. 
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Fic. 1. (Top) Gray shade display of a small sector of SMS 
visible data showing the clouds as the white areas. (Bottom) 
Computer output giving statistics on the clouds observed in 
top photo. The size is in satellite element numbers. 


(Reynolds and Vonder Haar, 1976). Figure 1 is an ex- 
ample of a computer-derived cloud count using SMS 
data. The image is reproduced in a gray shade format, 
and the number and size of each cloud are given. 

Prior to the launch of SMS-GOES, the ATS series of 
satellites provided geosynchronous visible satellite data 
(4 km resolution at subpoint) on a half-hour basis, so 
that at least cloud sizes, numbers, locations, and tem- 
poral variations could be determined during daylight 
hours. Several years of ATS data have been used to 
determine cloud populations and active areas of con- 
vection for the three High Plains Cooperative Proj- 
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Fig. 2. Cloud size distribution histograms for Miles City, 
Mont., area as observed from the ATS-3 satellite. The satel- 
lite resolution is 113 km s . 


ect (HIPLEX) field sites (Miles City, Mont.; Good- 
land, Kans.; and Big Spring, Tex.) being used by the 
Bureau of Reclamation (see Reynolds and Vonder Haar, 
1975). Not only can information be derived about the 
seeding potential of the sites, but decisions on the opti- 
mum location of rain gage networks can also be aided 
by having this type of information. Figure 2 shows re- 
sults from a cloud population study done over Miles 
City, Mont., using ATS-3 data. 

Other types of satellite data have also been used for 
these types of studies. Stock and Grant (1976) have used 
data from the low-orbiting, very high resolution (0.6 km) 
DMSP satellite to obtain cloud numbers and sizes for 
weather modification sites over the western United 
States. Although this information provides only a snap- 
shot of the clouds at local noon, even the very small 
cumulus can be detected with these data (see Table 1). 

Another aspect of using satellite information for de- 
termining cloud climatologies is that the general num- 
bers, sizes, growth patterns, and also cirrus extent will 
be known for both the target site and areas downwind. 
This will aid in determining to what extent seeding may 
affect the normal convective activity both at the site 
and downwind of the site. This will be discussed in 
more detail in a later section. 


There are problems in using satellite data for de- 
veloping cloud climatologies as there are with other 
techniques. Cirrus clouds can obscure lower clouds, and 
ground resolution problems may cause the smaller 
clouds to be missed. Computational time is extensive in 
processing the digital satellite data. Determination of 
cloud thresholds is also difficult. A man-computer in- 
teractive processing system is necessary for this type of 
study for accurate results to be obtained. However, we 
believe that the satellite provides a valuable tool to help 
assess the seeding possibilities in a given area. 

2) Real-time operational support 

As was mentioned earlier, near real-time high-resolution 
SMS-GOES satellite photo imagery is available to opera- 
tional field sites at half-hour intervals. This allows quali- 
tative analysis of both synoptic scale cloud features, as 
well as smaller mesoscale convective cloud developments 
to aid in the "nowcasting" of target site conditions for 
deploying aircraft or starting seeding operations. 
Reynolds and Matthews (1976) and Reynolds and Von- 
der Haar (1976) have described how geosynchronous 
satellite information has aided in supporting operations 
for HIPLEX (Bureau of Reclamation). It was noted that 
such convective triggering mechanisms as dry lines, 
vorticity centers, and fronts could be defined using the 
imagery and used for making the day’s forecast of con- 
vective activity (Fig. 3; see also Purdom (1974)). Since 
all of the HIPLEX sites are located east of mountainous 
terrain, some of the convection passing over the sites is 
generated over the mountains. The generation and 
movement of these cells and whether ot not they are in- 
creasing or decreasing in intensity can easily be moni- 
tored using the satellite imagery. Enhanced IR imagery 
can provide this type of information by monitoring 
CTTs to see if they are increasing or decreasing. Figure 
4 is an example of this type of imagery, and the change 
in CTTs can be noted by the change in gray level of the 
tops (see Corbell et al. (1976) for details). 

3) Use of digital and image satellite data in 

POSTEXPERIMENT ANALYSIS 

Operations debriefing and evaluation. Much of the 
postanalysis work to date has been in using imagery to 
help locate certain synoptic or mesoscale features that 
may have influenced an experimental day’s activities 
and to monitor their time histories. Reynolds and 


Table 1. Mean values for days with convective clouds, July 1974 



Miles City, 
Mont. 

Goodland, 

Kans. 

Big Spring, 
Tex. 

NE 

Colorado* 

Palmer 

Lake 

Divide 

South 

Park 

Upper 
Arkansas 
River Valleyt 

Average cloud cover, % 

4.9 

10.1 

8.3 

4.4 

11.1 

20.3 

27.0 

Average number of clouds per 10 4 km 2 

4 

6 

9 

8 

7 

18 

19 

Average size of cloud, km 2 

144.1 

172.1 

87.3 

54.7 

158.1 

115.2 

140.5 


* National Hail Research Experiment, 
t Leadville and Buena Vista. 
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Fic. 3. SMS-2 2 km visible sector image taken at 2045 GMT on 28 August 1975 showing a line 
of thunderstorms developing between Goodland (GLD) and Hill City (HLC), Kans.. extending 
to the southwest. The 2100 GMT surface observations are plotted. Note the dew point gradient 
across the thunderstorm activity. 


Matthews (1976) and Reynolds and Vonder Haar (1976) 
have discussed these types of techniques in detail. This 
type of research analysis allows stratification of a given 
day’s events into certain categories, which may be filed 
for later reference, regarding the time when a seed/ 
no-seed decision must be made. One of the most im- 
portant features that must be defined by this study is 
determining those days when organized rather than 
isolated convection will be present. With the use of 
conventional meteorological data, those days when or- 
ganized convection was observed with imagery to de- 
velop may be studied and classified. It is felt that 
ultimately, to generate meaningful increases in pre- 
cipitation, generation by seeding of organized, long-lived 
convective systems is needed. Satellite information may 
aid in developing the technology for identifying these 
situations. 

A very important observation that the satellite is 
capable of making, and one that is critical to the evalua- 
tion of any weather modification program, is a com- 
parison of the convective clouds behavior on seed versus 
no-seed days. This specifically addresses the problem of 
relative differences between cloud characteristics on the 
randomized seed/no-seed days. The visible digital satel- 
lite imagery can quantify such cloud characteristics as 
size, number of clouds, reflected brightness, anvil extent, 
cloud duration, movement (both speed and direction 
with respect to the environmental winds), cloud or- 
ganization (line cluster, etc.), cloud separation distances, 
etc. The IR data can monitor CTT (height) variations 


of clouds on the seed/no-seed days to determine both 
temporal and spatial growth patterns on the different 
days. This may be important in assessing whether a 
dynamic seeding approach is modifying the growth pat- 
terns of these clouds. The IR data will also be important 
for they will be the only data available after sunset. 
This technique will require several seasons of data before 
an evaluation can be made, but from present experience 
it should offer an important evaluation mechanism. 

The satellite data may also provide information on 
exactly how random a seeding program might be. This 
would concern itself with whether the "equality of the 
draw” was really equal. A comparison can be made of 
the general surrounding cloudiness (areal visible extent, 
vertical development (IR)) to determine if on the seed 
days more generalized convective development was ap- 
parent than on the no-seed days and vice versa. Thus in 
the evaluation of the seeding results, if a bias did exist, 
it could be eliminated before the final statistics were 
tabulated. This method might also lend itself to the 
determination of covariates for objective decision mak- 
ing of the day’s events. For example, if more well de- 
veloped cloudiness is observed, a look at the synoptic- 
mesoscale weather analysis may show a particular fea- 
ture that could then be classified for future use. 

Cloud parameter determinations. With use of digital 
satellite data, certain quantitative measurements can be 
made of individual clouds and cloud systems, and these 
results can be compared to other ground and aircraft 
observations. Emphasis must be placed on collection of 
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Fig. 4. (Top) The 0115 GMT, 22 June 1976, enhanced 1R 
image showing two large thunderstorms located in the Texas 
panhandle. (Bottom) The 0215 GMT, 22 June 1976, enhanced 
IR image showing that the northern storm has increased in 
height while the southern storm appears to have remained 
constant. 

digital rather than photo image data for these quanti- 
tative studies. There are several reasons for this: 

1) Digital data can be earth located, "navigated,” ac- 
curately for proper location of target sites and co- 
location with other data (Smith and Phillips, 1972). 

2) Data can be normalized for sun angle changes and 
viewing angles using bidirectional reflectance mod- 
els (Sikula and Vonder Haar, 1972). 

3) Data can be directly input into digital processing 
systems (see Section 4) for analysis, and this allows 
computer mixing of several different types of data 
as well as time-sequencing capabilities. 

Probably the most important parameter that must be 
measured in any weather modification program is pre- 
cipitation. Given an area as large as the High Plains 
region of the United States, it is very difficult to have 
enough rain gages available to obtain good estimates of 
rainfall. Radar coverage is very limited in these regions 
but is the next most accurate tool for use in rainfall 
measurements. However, recent studies in H1PLEX 
show a single Z-R relationship is very difficult to obtain. 
Over the past several years, work has been done in at- 
tempting to estimate rainfall from visible and IR satel- 
lite data (Griffith et al., 1976; Follansbee and Oliver, 
1975: and Scofield and Oliver, 1977). In an attempt to 
extend some of this work, we compared digital SMS 


satellite data (available every 7.5 min) to digital radar 
data over the Goodland/Colby HIPLEX sites (Negri 
et al., 1976) for 28 August 1975 during weak convective 
activity. The main conclusions from this study were 
that: 

1) with use of a reflected brightness thresholding tech- 
nique, precipitating versus nonprecipitating clouds 
may be defined; 

2) a lag of as much as 20 min may elapse between the 
time this brightness threshold is reached and a radar 
echo is noted; 

3) brightness does not seem to correlate well with any 
radar echo parameters such as area, volume, or in- 
tensity, but it does seem to be related to cloud size, 
and the center of maximum brightness locates well 
with the maximum radar reflectivity during the 
growth stage of the storm; 

4) the areal growth rate of a cloud may relate to storm 
severity and to determine these growth rates, rapid 
scan data are necessary (5-15 min interval data). 

These studies are in their very formative stages, and 
much work is needed to test the possibility of this pro- 
cedure. 

In relating cloud visible brightness to cloud parameters 
of interest, a separate study was needed to determine the 
relationship between changes in cloud brightness and 
simple changes in the cloud’s vertical or horizontal ex- 
tent (Reynolds et al., 1978). This is a very important 
problem if we are going to use satellite brightness data 
to quantitatively measure changes in intensity of pre- 
cipitation. Preliminary results seem to indicate that 
geometric factors outweigh change in cloud microphysi- 
cal characteristics in influencing cloud brightness 
changes. McKee and Cox (1974, 1976) have discussed 
this problem in detail on a more theoretical basis, and 
the question must be studied further. 

One other aspect of using quantitative satellite data 
is for input to and verification of cumulus cloud models 
(Kreitzberg, 1976). One-, two-, and even three- 
dimensional cloud models are being used in many 
weather modification programs. Ideal input to these 
models requires a closely spaced rawinsonde network 
for determining die atmospheric temperature and mois- 
ture structure and the vertical and horizontal wind 
components. With the advent of vertical temperature 
sounders on NOAA polar-orbiting satellites, high hori- 
zontal resolution (approximately every 70 km) tempera- 
ture and moisture measurements became available for 
use in these types of models. Hillger and Vonder Haar 
(1977) have used VTPR (Vertical Temperature Profile 
Radiometer) data in conjunction with the conventional 
rawinsonde network to obtain mesoscale temperature 
and moisture fields over the High Plains. Figure 5 is an 
example of the use of VTPR data to derive mesoscale 
fields of temperature and precipitable water. (Although 
vertical resolution from satellite data is poor, having 
gradients in these parameters over short horizontal dis- 
tances provides much information at the mesoscale.) 
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Fig. 5. (Left) Combined NWS and VTPR 500 mb temperature (in degrees Celsius) field for 9 
August 1973. The thermal winds for the 700-300 mb layer are shown as derived from the VTPR 
data. (Right) VTPR-derived total precipitable water field (in centimeters) for 9 August 1973 at 
1700 GMT. 


This type of information is now being integrated into 
HIPLEX so that it will be available in real time for 
input into cloud models or for analysis of mesoscale 
features likely to cause convective activity. When 
GOES-D is launched in the 1980s, it will provide the 
first geosynchronous vertical temperature sounder data, 
allowing soundings to be made at rapid intervals over 
small areas where convective activity is likely to occur. 

4) Specialized application areas 

Tropical cumulus modification study. Two general 
areas can be specified when discussing tropical weather 
modification activities. The first includes such activities 
as the Florida Area Cumulus Experiment (FACE) 
(Woodley et al, 1976), which is involved in the seeding 
of tropical cumuli to increase rainfall. The seeding 
hypothesis of the FACE program is to increase the or- 
ganization of the convective activity by seeding certain 
elements of a cloud line increasing the dynamics of the 
system so that it may precipitate for a longer period of 
time. As was mentioned earlier, Griffith et al. (1976) 
have been using satellite data to determine rainfall over 
the target area for a limited data sample. Other uses of 
satellite data are to simply pinpoint areas of existing 
cloud lines or clusters so that seeding aircraft can be 
deployed to these areas in time for seeding to occur. 
With the IR data now avilable from SMS-GOES it is 
possible to monitor these cloud lines and note any 
changes in CTT during seeding or the extent of cirrus 
for the seeded versus non-seeded storms. These IR data 


will also allow comparison to radar data to further ex- 
plore satellite-rainfall relationships. 

The other specific application area that is immedi- 
ately thought of is that of hurricane modification. Satel- 
lite data allow the only consistently available monitoring 
system for determining storm strength, motion, and the 
changes in these parameters. A more direct verification 
method to evaluate hurricane modification can be ac- 
complished by using geosynchronous satellite data to 
determine wind speeds from cloud motions within the 
storm system. With a clear view of the eye, measure- 
ments can be made of cloud motions along the eye wall 
(Gentry et al., 1976). Dvorak (1975) has developed a 
technique for determining die intensity of tropical dis- 
turbances that can be used as a type of climatology when 
seeding occurs to determine if any radical effects have 
taken place. Project Stormfury, which is to begin again 
in 1978, will rely heavily on the SMS-GOES satellites 
for the monitoring of these storms as well as for the 
evaluation of the experiment. 

Extra-area effects (downwind, etc.). The question of 
whether seeding effects are limited to the local area in 
which they are applied has been discussed as a political 
and scientific issue since weather modification activities 
began. There are several ways in which summertime 
cumulus activity and the modification of this activity 
can affect other areas outside the proposed target site. 
One method is that the cirrus outflow from the seeded 
storms could be more extensive than in non-seeded 
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Fic. 6. DMSP 0.6 km visible image showing thunderstorm 
development along the front ranges and mountains of Colo- 
rado and New Mexico. Note the cirrus blow-off and the new 
cell under cirrus. 


cases; this could decrease the local heating under this 
cirrus canopy, thus reducing convection. Satellites pro- 
vide a complete view of all areas both in and out of the 
target areas and can easily monitor the cirrus outflow 
of both the seeded and non-seeded cells to determine 
any differences in cirrus extent and whether convection 
is growing through the cirrus canopy or suppressed in 
areas where the canopies existed. It may also be that the 
cirrus can act as a natural seeding device by injecting 
large amounts of ice crystals into developing cumulus 
clouds to increase cloud growth. Figure 6 is an example 
of a DMSP 0.6 km visible image showing convection 
along the continental divide and front range areas of 
Colorado and New Mexico. Note the extent of cirrus 
outflow' of these naturally forming clouds and also the 
new cell development under the cirrus outflow of an 
earlier cloud. It is this type of monitoring that can aid 
in the analysis of extra-area effects of seeding programs. 

b. Wintertime cloud modification experiments 

I) Cold orographic and Pacific frontal 

RAINBAND CLOUD MODIFICATION 

The potential for increased precipitation from super- 
cooled orographic stratiform and frontal rainband pre- 
cipitation clouds is well established (Grant and Mielke, 
1967; Elliott et al., 1971; Chappell et al., 1971; Grant 


and Elliott, 1974; Gagin and Neuman, 1977). In fact it 
may be that these types of clouds offer the best possi- 
bility for seeding and for increasing precipitation over 
those areas where they exist. Some of the research ex- 
periments that have dealt with this type of cloud seed- 
ing are: 


Group 
EG8.-G Inc. 

Colorado State University 
CSIRO 

Desert Research Institute 
E. Bollay Assoc. 

Fresno State University 
Montana State University 
New Mexico State 
University 

North American Weather 
Consultants 
Utah State University 

University of Washington 

University of Wyoming 


Area for Field Studies 

San Juan Mountains of 
southern Colorado 
Central Colorado (Climax) 
Tasmania 
Sierra Nevada 
Park Range in northwest 
Colorado 
Sierra Nevada 
Bridger Range in Montana 
Jemez Mountains of New 
Mexico 

Mountains- near Santa 
Barbara, Calif. 

Wasatch Mountains of 
northern Utah 
Cascade Mountains in 
Washington 

Elk Mountain and Wind 
River Mountains in 
Wyoming 


Grant and Elliott (1974) have summarized the results 
from many of these studies and others to show that seed- 
ing effectiveness is highly dependent upon CTT (see 
Table 2.) Their studies show that increases in precipi- 
tation can be expected when seeding is performed on 
clouds with CTTs in the range of — 10°C to — 25°C. For 
clouds colder than — 28°C, decreases in precipitation can 
actually occur after seeding. During many of these pro- 


Table 2. Santa Barbara 2 (1967-70) average precipitation as a 
function of estimated CTTs for seeded and non-seeded cases.* 


CTT, 

°C 

Seeded 

Non-Seeded 

Precipi- 

tation 

Ratio 

Seed/ 

No-Seed 

Precipi- 

tation, 

mm/band 

No. of 
Cases 

Precipi- 

tation, 

nun/band 

No. of 
Cases 

-17 

10.7 

13 

5.1 

22 

2.10 

-18 

18.8 

12 

6.4 

23 

2.96 

-19 

16.5 

24 

6.4 

21 

2.60 

-20 

12.5 

22 

6.6 

17 

1.88 

-21 

9.9 

21 

6.4 

14 

1.56 

-22 

10.4 

21 

7.6 

11 

1.36 

-23.5 

8.6 

10 

7.1 

8 

1.21 

-24.5 

10.2 

10 

8.9 

8 

1.14 

-26 

11.4 

10 

9.7 

10 

1.18 

-27.5 

15.8 

9 

12.5 

7 

1.26 

-28.5 

15.0 

7 

8.9 

5 

1.40 

-31 

9.4 

6 

5.6 

4 

t 

-33.5 

7.4 

3 

5.1 

2 

t 


* Values of precipitation for moving means of two temperature 
steps of 2.5°-3°C are centered on the indicated CTT. (Reproduced 
from Grant and Elliott, 1974). 
t Sample too small to be significant. 
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SATELLITE CLOUD -TOP TEMPERATURE. *C 
(MEAN) 


Fig. 7. Relationship between mean satellite-observed CTTs 
and radiosonde-determined CTTs for the entire data set of 
1973-74 for the San Juan Mountain area. 

grams, estimates of CTT were derived from either 
rawinsonde data or interpolated 500 mb temperatures. 
These have frequently proved to be inadequate in that 
only 3—5 rawinsondes can be launched per day and cloud 
tops cannot always be adequately deduced from the 
soundings alone. In addition, the 500 mb temperatures 
do not represent CTTs in many instances. As was stated 
earlier, Dennis el at. (1973) noted that one of the most 
beneficial uses of satellite data would be in monitoring 
CTTs from a geosynchronous satellite to provide half- 
hour information on CTT and its changes. 

Real-time operational support. The term “opportunity 
recognition” has been used in the field of weather 
modification to designate those environmental conditions 
in which it has been determined that seeding will have 
positive effects. As we have pointed out, the CTT is a 
very important parameter when discussing these types of 
cloud systems. Another important question of "oppor- 
tunity recognition” concerns the spatial and temporal 
location of seedable clouds. The inhomogeneities of 
actual clouds and their seedability are so great that con- 
siderable variability generally exists over even moder- 
ately sized weather modification targets. The ability of 
IR satellite data to provide CTT information has been 
demonstrated by Dumont et at. (1974) in a preliminary 
study using NOAA polar orbiter satellite data. Using 
the IR portion (10-12 >*m) of the scanning radiometer 
(SR) on the NOAA satellite, they were able to measure 
CTTs for clouds over the San Juan Mountain area of 
Colorado during the Bureau of Reclamation's Colorado 
River Basin Pilot Project. Figure 7 is a comparison of 
satellite-determined CTT and radiosonde-determined 
CTT for all cases available. The results show that on the 



Ftc. 8. Enhanced GOES-W IR image from a digital display 
for 1745 GMT, 15 March 1977, showing a Pacific storm enter- 
ing the California coastal region. Darker gray represents 
CTTs from —43° to — 47°C; medium gray, from —39° to 
— 43°C; and white, from —31° to — 39°C. 

mean, the satellite observes a warmer cloud top (+1.8°C) 
than does the radiosonde. Although such problems as 
time differences between satellite and radiosonde obser- 
vations, cirrus clouds overlying the orographic cloud, and 
cloud emissivity problems affect these comparisons, this 
preliminary study demonstrated the general feasibility of 
using satellites for this type of measurement. A second 
season (1974-75) of data is now being processed; NOAA 
polar orbiter data, again for the San Juan Mountain re- 
gion, are being used. During this season, the first geosyn- 
chronous IR data became available from SMS-1 posi- 
tioned at 75°W. Digital data from this satellite are avail- 
able, and maps of CTT over the San Juan Mountains 
can be generated on a half-hour basis for a 24 h period. 

To provide the real-time support necessary for field 
operations there are several methods in which the cloud 
top information could be provided. The first is using 
the enhanced IR imagery data where a specific gray 
scale would denote clouds with temperatures in the 
seeding window. These images would be available to the 
forecaster on a half-hour basis so that he could decide 
not only in what locations out of the entire storm com- 
plex seeding should take place but also when to begin 
and end the seeding operations. This is particularly 
important since several studies have strongly indicated 
that the main seeding effect is an increase in the dura- 
tion of precipitation (Chappell et al., 1971; Gagin and 
Neuman, 1977). 

A second method would be to strip out the specific 
area of interest for the seeding activities from the digital 
data taken directly at the satellite ground station and 
either to map this information with a computer printer 
or to display the data on a CRT using a video recording 
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Fig. 9. SMS-2 2 km visible image for 1645 GMT, 9 March 1977, showing a 
Pacific frontal system situated over the Sierra Nevada. 


ORIGINAL PAOR II 
OP POOR QUALITY 


system. These types of data could be transmitted to a 
field site by phone line to give the forecaster a constant 
update of CTTs. More on this video-digital display capa- 
bility will be given in Section 3.b. 

Use of digital satellite data in postexperimental analy- 
sis. During the postexperiment phase of the seeding pro- 
gram the digital SMS-GOES imagery collected can be 
used to stratify each day’s seeding event into CTT cate- 
gories. This is being done for the 1974-75 Colorado 
River Basin Pilot Project experiment using NOAA 
and SMS data and will be done for the Bureau of 
Reclamation’s Sierra Cooperative Pilot Project, which 
began during the 1976-77 winter season over the cen- 
tral Sierra Nevada of California. The Bureau of Recla- 
mation is now involved in a preliminary investigation 
of the characteristics and natural variability of con- 
vective rainbands as they impinge on the Sierra 
Nevada of California. In cooperation with the Bureau 
of Reclamation, Colorado State University has recorded 
digital SMS full resolution visible and IR data con- 
tinuously from the time the bands are just off the coast of 
California to their passage through the Sierra Navada. We 
hope to obtain information on the capability of the 
satellite to determine the accuracy to which the top 
temperatures of these rainbands can be measured, how 
these temperatures vary with time, and the effects 
topography has on these bands (see Fig. 8). As verifica- 
tion, a cloud top aircraft is being flown over these bands 
to estimate their height and temperature. Photographs 
will be taken, and the presence of clouds, such as cirrus, 
that may overcast these bands will be noted. As was done 
for the summertime projects, radar-satellite comparisons 
will be made to determine if any relationships exist and 


whether the satellite data can be used to quantitatively 
determine precipitation amounts. 

Not only can CTTs be determined, but an assessment 
of the types of synoptic situations that give rise to 
seedable clouds can also be determined through the use 
of satellite imagery and conventional meteorological 
data. The satellite allows an overall view of the synoptic 
situation (Fig. 9), which can allow categorization of 
optimum seeding situations and can be used as a fore- 
casting aid for future events. 

The same application of the satellite data can be made 
in the wintertime program as was made in die summer- 
time programs. This relates to the relative difference in 
cloud characteristics of the seed/no-seed events as well 
as determining any biases in the randomized seeding 
events ("equality of draw”). For the wintertime clouds, 
such observable cloud features as stratiform layers versus 
presence of imbedded convection, CTT variation and 
magnitude, and visible cloud appearance (i.e., fuzzy (ice 
present?) versus distinct cloud edges (possibly only super- 
cooled water clouds), etc.) could be analyzed for these 
different days. Chappell et al. (1971) have found that 
for the Climax Experiment, seeding tended to increase 
the duration of precipitation rather then the intensity 
of the precipitation event. It may be that by either 
monitoring the CTT (height) fluctuation through a 
storm period or determining the duration of the oro- 
graphic clouds’ fuzzy appearance, the satellite data might 
be used to decide if the precipitation event is of longer 
duration or more intense in nature. 

For determining whether a true random sample has 
been chosen for a given experimental period, the satel- 
lite information can provide measurements of overall 
cloud extent both in and out of the target region. The 
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mean CTTs for the entire storm complex can be de- 
termined, and the direction and speed of cloud motions 
over the mountain barriers and indications of band 
structures (associated with either fronts or mesoscale 
features) can be observed for all seeding events. If any 
biases are present on the seed or no-seed days, then they 
can be removed before final analysis. Again, these same 
measurements might aid greatly in a covariate determina- 
tion technique. 

2) Specialized wintertime applications 

As was mentioned with the summertime program, there 
has been concern with the effects downwind from the 
wintertime seeding programs. Brier et al. (1974) review 
much of the work carried out to date in determining 
the downwind effects from seeding cold orographic 
clouds. Many of die results to date show varying degrees 
of effect, from large increases to slight decreases. Brown 
et al. (1976) have recently summarized the results of the 
seeding program over Santa Barbara, Calif., where 
Pacific frontal rainbands were seeded over a 7-year pe- 
riod. Their conclusions stated that the primary cause of 
extra-area effects is probably a dynamic intensification 
of organized convective activity that produces increases 
in precipitation — 150 km downwind from the seeding 
source and ~30° to the right of the 700 mb wind flows. 
With satellite data it should be a fairly straightforward 
problem to monitor this area in an effort to determine 
whether the activity is more intense or continues for a 
longer period of time. 

Mulvey and Grant (1976) have proposed a separate 
mechanism for extra-area effects downwind from the 
Climax Experiment in Colorado. Their hypothesis is 
that seeding material carried up into the orographic 
cloud is not entirely used and is transported downwind 
along with unused ice crystals. The seeding materials 
and ice crystals that survive the downwind transport can 
seed lower upslope cloudiness along the Front Range of 
the Rockies, which are warmer and deficient of the 
proper number of ice nuclei for efficient precipitation. 
This inadvertent seeding then should supply the cloud 
with an increased number of nuclei for increased pre- 
cipitation. Scheetz and Grant (1976) have prepared a 
climatic estimate of the frequency with which High 
Plains upslope clouds are seedable. With both visible 
and IR satellite data it is possible to monitor the ice 
blow-off from the orographic cloud. In the area where 
this intersects the upslope cloudiness, a change in both 
the visible appearance of the cloud and in top tempera- 
ture can be monitored to see if this seeding has any 
effect. Without aircraft data this is the only way to 
monitor this type of extra-area effect directly. 

3. Present and developing technology in satellite 

observation systems and data processing 

a. Present and near future capabilities 

We have discussed many of the present capabilities of 
meteorological satellites as they relate to existing weather 
modification activities. However, we have certainly not 



Fig. 10. Proposed geosynchronous satellite coverage for the 
late 1970s showing that all proposed WMO Precipitation En- 
hancement Experiment sites will be in view during all local 
times. 

discussed all the widely varying activides and uses of 
the data. There is one important area that should be 
discussed as it relates to the World Meteorological Or- 
ganization (WMO) Precipitation Enhancement Experi- 
ment (PEP) (List, 1976). Site selection is one of the first 
priorities of the WMO. Although there are many politi- 
cal, social, and economical influences on these decisions, 
one overwhelming factor must be the requirement for 
“seedable” clouds to exist at the site frequently enough 
that seeding can be done. The satellite data can observe 
die presence of clouds and their appearance, i.e., strati- 
fied (stable) versus convective (unstable) over a particular 
site, aiding in the decision of whether to use a static 
(microphysical) seeding approach or a dynamic seeding 
approach. As was mentioned earlier, satellite data can 
provide this type of information on a worldwide basis. 
There exist many years of polar orbiter data with after- 
noon cross-over times that can and have been used to 
monitor convective activity. Also, on board many of the 
satellites are IR sensors allowing the measurement of 
CTTs of clouds large enough to fill the field of view of 
the sensor, so that it can be determined whether the 
clouds are cold enough for seeding materials to be effec- 
tive. 

By 1979, in cooperation with FGGE there will be four, 
possibly five, geosynchronous satellites girdling the globe, 
(Fig. 10), allowing for the first time high spatial and 
temporal frequency satellite data continuously for the 
whole globe. These types of data will prove invaluable in 
assessing future weather modification sites both for sum- 
mer and for winter programs. 

We should mention here that there is one other im- 
portant use of the SMS/GOES satellite that can be a 
real-time source of information besides providing satel- 
lite visible and IR data. This involves the communica- 
tions channels on board the satellite that can interrogate, 
every half-hour, the transmitters of ground-based data 
collection platforms (DCPs). Such observations as tem- 
perature, humidity, and wind can be obtained and 
transmitted down to the ground station and displayed 
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on the aforementioned remote terminal Also, rain gage 
data or stream gage data can be monitored, especially 
for analysis of seeding effects or for flash flood alerts 
and snowpack conditions, to see if additional seeding or 
suspension of seeding is necessary m a given area This 
capability of SMS should be integrated with the digital 
imagery to provide extremely useful and timely informa- 
tion to all aspects of field operations and forecasting 

b Real-time satellite data analysis and 
display systems 

As has been mentioned in the previous discussions, of 
equal importance to the gathering of new data from 
new systems is the development of quantitative digital 
data analysis methods These digital data allow the 
monitoring of cloud brightness, GTTs, and cloud mo- 
tions In addition, accurate navigation of the data can 
be performed through the recording of the line docu- 
mentation being sent by die satellite (Smith and Phillips, 
1972) This includes the housekeeping information of 
the satellite such as the geometry of the earth, sun, and 
satellite and time marks for each scan line This in- 
formation can then be input to an analytic model for 
navigation of each visible image to ±1 (1 km) satellite 
element (relative picture to picture) Knowing the earth, 
sun, and satellite locations also allows brightness nor- 
malization to be made, which can account for changing 
sun angle effects and can try to compensate for cloud 
anisotropy 

Generally, satellite data, to be of maximum usefulness, 
must be available in real time. Over the past several 
years, development of digital data-processing systems 
with man-computer interactive capabilities and CRT 
display have been developed The following is a list of 
several of these systems now in use throughout the 
United States 

1) AOIPS — Atmospheric and Oceanographic Informa- 
tion Processing System (NASA Goodard Space 
Flight Center), 

2) LARS — Laboratory for Applications of Remote 
Sensing (Purdue University), 

3) McIDAS — Man Computer Interactive Data Access 
System (University of Wisconsin, Space Science and 
Engineering Center), 

4) MMIPS — Man-Machine Interactive Processing Sys- 
tem (NOAA/NESS), 

5) ADVISAR — All Digital Video Imaging System for 
Atmospheric Research (Colorado State University), 

6) ESIAC — Electronic Satellite Imaging Analysis Con- 
sole (Stanford Research Institute) 

These systems, when located with the satellite ground 
station, can provide real-time processing and analysis 
of cloud growth rates, cloud motions, change m top 
height, etc, that can be useful to field operation. (Real- 
time navigation is now in use at the Direct Readout 
Ground Station at the Space Science and Engineering 
Center m Wisconsin) These systems operate through 
input of digital satellite data, which can be stored on 


digital systems and displayed sequentially through rapid 
refresh of a CRT monitor Thus, 2-3 h of satellite data 
can be "looped” through, allowing cloud motions to be 
determined and thus providing information on meso- 
scale divergence, convergence, etc Also, the operator can 
interact directly with the data to enhance through color 
or black and white shading methods, he can zoom in on 
small features or rotate or translate the image and can 
determine, through use of a joystick and cursor, CTT 
or brightness of any cloud he so desires Through recent 
and ongoing advances in solid state electronics, these 
types of systems can be remotely placed in the field 
away from, ground stations with just a telephone hook-up 
to a keyboard and CRT necessary for data access and 
manipulation. At the present time, work is under way 
in developing such a system that can meet the needs of 
not only weather modification programs but also severe 
storm forecasters and many other satellite data users, 
including industrial meteorologists (Bristor and Raynore, 
1977) 

4. Future advances in satellite observing systems 

Future advances m satellite technology should he as im- 
pressive and as rapid as those that took place during the 
last 10 years One important new development that 
cotild provide dramatic technological advances will be 
the Space Shuttle The Space Shutde will put much 
larger sensor payloads into orbit either by constructing 
them m space or by launching them on the Shuttle and 
then boosting them into geosynchronous orbit by space 
taxis Thus, large microwave (passive and active) an- 
tenna systems can be launched allowing direct measure- 
ments of precipitation over land and oceans and also 
allowing measurements of water vapor profiles at very 
high horizontal resolution (~10km) and having 2-5 km 
resolution in the vertical Shenk and Kreins (1975) dis- 
cuss some of these future breakthroughs that will come 
about in the 1980s 

Not all of the advances will have to wait for the Space 
Shuttle Planned for launch around 1981 will be the 
GOES-D satellite equipped with VAS (VISSR Atmo- 
spheric Sounder) This will be the first sounder placed 
into geosynchronous orbit and will provide soundings 
at — 30 km spacing every half hour for a 750 km north- 
south latitude band and entire earth width The VAS 
will have nine channels, sensing temperatures (resolu- 
tion to ±2 5°C) below the 100 mb level (thus tempera- 
ture every 100 mb) and two water vapor channels for 
determining precipitable water 

The next generation o£ satellites being considered 
will be Stormsat (Shenk and Krems, 1975), which will 
improve both the imaging and the sounding capabilities 
of the GOES series of satellites This will be a three-axis 
stabilized spacecraft, which means it will always be point- 
ing toward die earth and will be mechanically stepped 
both in the east-west and in the north-south direction 
Thus, very rapid scan information over a 750 knr area 
will be available Also, soundings using microwave 
channels can be made in overcast regions with 2°— 3°C 
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temperature resolution at a 50 km interval every 30 min, 
and these microwave channels will allow moisture 
determination (10-15% accuracy in relative humidity) 
for 30 km intervals every 30 min This type of informa- 
tion can be used directly in two- and three-dimensional 
numerical cloud models for both initialization and 
verification 

There are now even plans to fly radiometers that will 
determine something about a cloud's physical properties, 
l e , cloud top pressure level, the density and phase 
(i e , ice versus water) of condensed water in the clouds, 
a drop-size parameter, and the cloud thickness, both 
optical and geometrical 

5 Summary 

We believe that satellite information has the greatest 
utility for weather modification application when it is 
incorporated into a research operation along with other 
data sets Over the past several years, satellite technology 
and analysis technology have improved greatly The last 
section of this paper was designed to show the projected 
advances that may take place in the next decade, even 
though they may be overly optimistic Satellite sounder 
systems offer one of the highest potentials for applica- 
tion to weather modification programs, and it is hoped 
that, through the use of multispectral instruments, a 
high degree of vertical (albeit still crude compared to 
radiosonde) and horizontal resolution can be obtained, 
both m temperature and in moisture Continued re- 
search is needed on quantifying observed features from 
the satellite image, quantification must be such that the 
numbers can be incorporated into numerical models to 
help predict cloud development, seeding potential, and 
possible extra-area effects 
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1. lNiTtODUCUOH 

Sn tell xtc i-ngcry oroduced b>* Lhe 
Defense ;h tc.orologrral Sartllifo P> ogr im >(D?ISP) 
hu, piuvided very high rL.ool.it ion visual images 
of LhL t.e>LOru United States oner daily during 
d rylipht hours from Gael) satellite. ihis study 
utilizes D’IRP dan available foi t.ui summer of 
lA^A to construe!. o cloud cli.aaf ology of areas 
where wo a her ircdif lent ion or cloud physics 
experiments are currently being conducted or 
pljnneJ. Those areas include the enter High 
Plains hxpcririeni. (MTU X) sites at Lolny, 

Kansas (CBY), Big Springs, Texas (I3GS), and 
Miles CiLy, Montana (IILS); Soutli Par!:, Colorado 
(SPk) , wheie tup South Pat It Area Cu-'ulus 
Experiment (SPACI.) is being conduct ca; and the 
National Hail Research Expo merit (NliKI.) area m 
northeastern Colorado. The Primer hale Divide 
region (I’LD) located bt tween Denver and Colorado 
Springs, Colorado, and the upper Aikansas P.jver 
Valley (ARK) between Buena Vista and I.er.dville, 
Colorado, vt_ic chosen for tiller proximity to 
cuncnt weather modification protects and tneir 
potential as sites fur new weather modification 
studies. 

The satellite used m this *'tudy 
is sun-dynehl onous polar orbiter passing over 
the study jioas at approximately 1200 local time. 
Digital data Lapes were produced from the trans- 
parencies and a computer analysis v»vs accom- 
plished, pLoduLing areas and numbers of clouds 
for each area and day as well as pioviuing data 
for cumulative averages over selected Line 
periods for each studv area. Cloud locations 
and sizes were si tat if red with respect to a 
subjective cloud classification scheme as tell 
as with respect to synoptic weather conditions. 
The relationship between cloudiness and topo- 
graphic lenLincs was also examined. 

2. DELI1JI Ar'IOK 01‘ STUDY AREAS 

rigiice 1 shows til.’ location and 
shape of each of the areas analysed in this 
sludy. The tlnco Illl’t.IX sites arc circles with 
a 112 Ian (60 n.m) radius from Ill.il, CRY, anu 
BCS. 1 hose sites are included teen though a 
similar study has already been performed for 
the suLirn is of 1972-7A using data from tiic ATS- 
III satellite (Reynold > and Yonder Jl.rnr, 1975) 
since the 3)1151’ s iccllite da La has better reso- 
lution. The NRItU area war described for analysis 
purposes by using the 1973 borders of the NURD 
sLudy area (! CAR, 19 /A), lire PI D area was 
chosen ,;itll the western border uL the tdge of 
the I’ront Range of the Colorado Rockies. The 
eastern border is the mosl easterly exLerit of 
the divide itself, about 330 kn (70 n.m) east 
of the mountains. The umth-south extent of the 


study area i’ rqual to Lhe wrdth of the divide, 
from the souLh edge of Denver Lo the norLb edge 
of Colorado Springs, a distance of about 90 k.> (50 
n.m). The wt stein boundary of this area vis 
deter nine. I to have "hot spots" Lr it favoi ed 
thunders icin’ foir.iaLion (Ilonz, 1973). lhe 
South Park st w!\ area (SPIC) r as chosen to provide 
satellite data coverage of tno gi-ogi ankle. I ntei 
of interest in the South Park Aren Curmlus 
lxperii'icnt (53‘ACE) pre ,er.tly being candu.-U « b” 
Colorido State University, lht are r was <K -igned 
Lo include a] 3 of the lopo&'Mp’.ira J features 
known as South Park as well as to have i 
gulav shape for case in compel 'T processing of 
aata. The area for which satellite data !S 
analyzed is bounded by tno cro'-r of the Mosquito 
Range to Lnc wesL, the Kenosha and Tr’rynil 
Ranges to Lhe east, a latitude line through 
liuenn Visln or the south, and the Continent tl 
Divide Lo the north. lire area is jpprorimitcly 
45 Urn (23 n.m) cost— west and 65 km (35 n.m) nanli— 
south. The ..evencn area eucompasse 0 the uppir 
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A' i ii!- is Kivri villoy hounded rm'lilv by 1 1>_ 

Coal menial Hindi no) 111 ol It idviile, on tin 
north, the lautii.le of huen i VIst i cm t ho Ruulii, 
the Mr itch r mge Lo the to .1 .nut the Ilosquj l o 
r.nijo ro tin* o.i. i. Hie are i of studv hot onus 
n.irrowtr to tie iiorLh .is lIio i ivor valley n.ii rovit . 
Tho appioxin 1 L 0 sj.o is 65 l.w (15 n.m) noi Lh- 
south ind .ivoi iges 30 Km (36 11 . 1 . 1 ) t ist-wst. 
ll.ls area was chosen lor study because of 11 s 
pt oxiiuily to SPACE. 

3. SAJuLUir SIKSOli MAI A 

The IlMSP sifcHilo that produced 
the data used In this >tudy is .1 polar oibilcr 
with an evernf e spacecraft height of 853 kri (450 
n.m). Hie orbit i« circular with an angle of 
Incliiuation of 93.7°. This inclination angle 
was selected lo insure Lne satellite orbit would 
be sun-syncln oituuc at this altitude. (Dickinson 
er al, 1974) The nodal penod of Lhis sutt- 
synchionous orbit is 101.56. Since Lhe earth 
revolves unoer the orbit to the cnsL, and the 
oibit processes slightly to the cast, each 
nodal crossing is appi ovinia Lely 25.4 degrees west 
of the previous crossing, lhe sensor scans ?6„6 
degrees of latitude across the subtiack. The 
very high resolution visual rod 1 one tor chat 
provider Lhe data for this sLudy has 1 spectial 
range of 0.4 to 1.3 micrometers. lhe resolution 
of the imagery is nominally 0.6 km (.33 n.m) at 
the satellite subtract.. This degrades slowly to 
about 3.7 km (2 n.m) at Lhe edge of the data due 
to foreshortening. It is interesting to note 
that while each nodal crossing is 25.4 degrees 
west cf the pi evious crossing, the width of Lhe 
scan noth js 26.6 degrees of latitude. Tins 
causes a 1.2 degrto overlap at the equator 
betw«.en successive pa^sc > winch becomes 
prcgi ersively largci as Lhe satellite approaches 
eilher pole. Although the resolution at the 
edge of Lhe scan is degraded by foreshortening, 
if the area of intere-.L falls near Lhe cd ;c of 
two successive passes, data can be obtained on 
cloud novoinenL and cloud growth occurring during 
the 101.56 minute orbital period. The satellite 
used for this study (1)6111* 8531) has a noontime 
sun-synchronous 01 nit. Tor the area of interest 
in this paper, 1 latitude range from 32 degrees 
to 47 degrees norlh, the longitude of the 
ascending node of Lhe satellite is about 7 to 
13 degrees east of the longitude of lhe satellite 
subtiack. Since best resolution is obtained when 
the satellite subLi ack passes ovci or very close 
to the area of interest, the local solar time 
of satellite passovtr varies from about 20 
minuLcs before local noon al JIGS Lo about 35 
minuLcs before local noon at IIL5. 

4. RITE IDEM IFICATIOK 

The method chosen for accurately 
locating the study areas on Lhe satellite imagery 
involved producing an overlay containing several 
properly located geographical reference points 
large enough to lie easily visible on the imagery. 
The study areas lore then accurately positioned 
on the overliy with respect Lo these landmarks. 
Since most days during the summer had large areas 
of clear skies, the meLhod worked well. On days 
when Imdmiks were obscured by cloud, the study 
area overlay was accurately positioned using the. 
grid supplied with the satellite linger}. It 


has been estimited tint laud or tlnud fcatuiis 
can be JoLiled to wrthrn 2.S km (3.5 n.m) across 
Liurh rod 5.6 km (5.0 n.i») on tho center line 
using Lhe grid (liicUieoii cL tl . 19/4). 5 he 

landmark method appealed to work equally Will. 

5* SYNOPTIC l'l AT >11 It PATTI k, 3 

In or ilt r to provide a better 
bncl gi cm nd for the cloud 1 limnology daCi, 1 
biicf MMiiaiy of sigml leant weather p’LUiii. 
for 6u.nu. 1 3974 Js included hei e The wtatbci 
of June 19/4 was ch ir.u Lei iced by a tell 
developi d broad rid-e over the western si .11 is 
which was accor.p’nied by per-i stent cli ir skies, 
warmer than normal nuximum temperatures, and 
genet ally below norm-J precipitation except lot 
portions of western Kansas and eastern Color ido. 
The first week of June vis cool mid wot 'Hit 
snow falling m cite ( olocado nioun Largs June 8 
and 9. lhe firm and di v pattern hec.rc 
established the second week of June and main- 
tained itself for the rest of the month 
(Taubensce, 19/4). The mean 700 mb ridge 
migrated slightly eastward during July tad 
became estiolishe.l over the Central Pints. The 
1I1PI.ET areas were warmer and dnet than uornnl 
with the greatest precipitation dofic h n<-y cuing 
observed in tho CliY and BCR areas. The Colorado 
nouniain and plains areas general 1 v reunited near 
normal tci iporatures in Lite nioumains ("aguer, 
1974). There were walked circulation 1 hinges 
from July to August over the western si ales and 
Great Plains area. Tiouglung it 700 1 . 1 I 1 ropl ued 
Che mean ridge over the central plain 1 resu 1 1 Ing, 
in cooler thin normal It mpi 1 lturer ovoi the 
entire t egion. Precipitation was gtn'i-ljy 
above not iiial in the southein and eenltai pi tins, 
near normal la eastern Montana , wiL’i all th 
areas in Colorado involved in this stud/ 
receiving less than normal precipitation 
(Dickson, 1974). 

In Colorado during the summer 
months, storms typically form ovei the mount ams, 
translate and propagate eastward onto the ihgn 
Plains during the afLetnoon and evening. lhe 
mechanism for storm formation is a combination 
of diurnal boating and local mountain-valley 
breeze circulation. The propagation seems Lo he 
controlled by lo* level fluxes of lioaL and 
moisture, Lite stability, and the atmo-plu. rit wind 
s true tut c (Erbes, 1976). Since this study is a 
climatology 1 eprcsentlng only one point in time 
during the da_, , this effect is reflected lit tho 
size and number of clouds observed in the ai eas. 

6. DATA ACQUISITION 

Veiy high resolution visual data 
for the months of June, July and August, 3974 
were obtained from DiiSP archives maintained by 
Lhe University of tlisconsin, in the fot tn or 
transparencies produced from original •dtcllilc 
analog dila. Eighty-eight d iys of very high 
resolution (0.6 ki,i) d ita were available and four 
days of high rcsolut ion (3.7 Ian) data making a 
complete sol of 92 days of satellite imngeiv 
near local noon available for analysis. The 
scale of Lhe 0.6 In data is 1:7,51)0,000 and 
the scale of the 3.7 km is 1:15,000,000. Two 
data reduction melhods were employed. fhe first 
involved subjective clas sif icarions of clonus in 
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the project arear, by cloud typo and cloud amount.. 
The second method involved digitizing the images 
ot protect areas on the Optical Rata Digitizer 
and Display System (013-*) at Colorado State 
University. The components of the Od 3 system 
are a videcon camera, a television screen, a 
snail corpuLer, and a typewriter for transmitting 
institutions to the :yslem. The computer was 
programmed to sense brightness levels from the 
screen and convert these into digital data. The 
data is stored on magnetic tape for later 
pioccssing. l'oi the digiti-mg process, one 
pixel or one data bit was equated lo Lite resolu- 
tion of the im-igei y such LlnL the area repre- 
sented by one pixel was equal to Lite highest 
resolution of the sensor. 

7. DATA CLASblflCATIOM 

Subjective classification of the 
data with respect to synoptic and itasoscule 
inflt’enccs was per foiled by inspection of the 
satellite rmageiy. Table 1 shows the c iLcgoties 
used and the percentage of days during the 

pci tint each of these cloutt types and 

weather Lypes was presettt. Three major cate- 
gories , cold front, upslopc, and mesoscalc, 
were considered. Mcsosrate systems wetc Cut Liter 
divided and values m paicntbeses show the 
percentage of mesoscalc day. affected by the 
respective cloud types. These divisions were 
patterned after a cloud study of the It 1PM X aicas 
by Reynolds and Vontler lliar (19 /5). The .small 
cumulus category was added because of the 
increased resolving power available from the 


D!ISP sensor, Mesoscnle processes dominated 
the weather piLtern over the region for most of 
the summer. The Lwo most common tnesoscnlc cloud 
types observed on the noontrme nnagei y at all 
areas were the small cumulus, and the isolated 
cumulonimbus. Trom Tabic 1 it is seen that 
imbedded ctir.nl oniiuhl account for 9% of the 
mesoscalc occurrences at '1LS, 2 Z at CI1Y, and 97, 
DCS. If the small cumulus category is removed 
to male the ttula more cor.partblc with the 
Reynolds and Vender tiaar study and the raoso- 
scale percentages are t ccalculaLed, Imbedded 
cusrulommbi then account for 1 17; of the r.esosc ale 
weather occunences a t HIS, 3 % at CRY, and 14% 
at RGS. These perceiitaj.es ate l.nger than those 
determined by the pi cvimtsly mentioned study, 
no doubt clue in part to the net Lor resolution of 
the b'til’ data. I rental aetivny was more 
appaicnt at >iL5 and CUT Lisin ot any of the other 
areas. In both cisos, most of this activity wis 
observed during the month of August. Surface 
frontal systems did not appear to pcnelinLa ns 
far south as ",CS until tlu last few dn\s of 
August. A large number of clays when no clouds 
could be seen occuircd at every iicj C'cept SPK 
and ARK. The fewer number of clear dais at these 
statiors can be considered an indication of the 
effectiveness of the serf ice. be I Ling, mountain- 
valley circulation regime an pioducinj, cumulus 
clouds. These statistics, though in .similar- 
format lo Reynolds md Vender ilna r (19/5), are 
not directly co.upar ible due to the inclusion of 
1972 and 1973 data in the* previous study. 
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DATA AKAIV’.JS 

The digitized tl it a dis.ui.id in 
Section 6 has heel Mibjected ip a mi"hu] of tests. 
Th e first m*’oJ ved delci Niini'j> cloud n *os end 
cloud min, bo is for onch of 1 ho jica'i for llio 
sunmer of 1074. fins is done hy selecting a 
bi j ghtness thieshold lor the < loml border .aid 
summing ail of the pi iJs within that border to 
determine Urn lol.il area of the clouds. 


'ilia cloud brighi ness thieshold was 
chosen by comparison of pruUimU, of clue I 
pictures wiib pi inLoots of at. Lu if digitized 
brightness values for each fiLc, Hie computer 
senses each closed threshold brightness value 
contour as a cload. The total uuancr of clouds 
per sjte, the total cloud ooveiagi foi the day 
studied, and t lie average si/e of the clo..ds 
within the site area were calculated. Table 2 
presents a summary of this output for Cue month 
of July 1574. Three tabulations are presented, 
the average percent of cloud cover, tno average 
number of clouds normalized to an area of 
10,000 hrn^, and the average si/e of the cloud m 
hm^ for each study area. At both MLS and CiJf 
the average cloud size was relatively lai gc but 
the number of clouds per 10,000 km^ was low. 

These data would seem to indicate that when 
convective ictivity was present at noon it was 
already well developed. The data for liGS in 
conjunction with Table 1, showed that large 
numbers of small cumulus clouds predominated at 
noon but that I be few eases when larger 
cuwulou jinbi wet e present increased the average 
individual clou i ai ea. Min'd', PLD, SPK, and Alt!:, 
when considered together, corroborate the 
description of a typical woollier day over the 
mountains and high plains of Colorado. Since 
SP1C anu ARK an both located m the Mountains, 
a noontime climatology should show a relni ively 
high amount of clouds m these areas. This was, 
in fact, the ease with SPK and AI.K having an 
average cloud cover of 20.3% and 27% respectively. 
Thete wete a large number of clouds present with 
individual clouds ha. mg relatively large Si.cs. 
PLD, located adjacent to a "hot spot" for thunder- 
storm development had less average cloud cover 
than the mountain areas and fewer clouds. Since 
all of the activity analyzed at PLD for tins 
month was convtctive, the clouds that weic 
present aL this Lime were alieady well developed 
and they tended to be positioned at the western 
edge of the atea. NIIRE, being further away fron 
the mountains had a low average peiccnl cloud 
cover, few clouds, and a small avci age cloud 
si/e. On most d tys during the month the 
mountain tlmndr rsloin cyi le did not begin early 
enough for Inge cells to icicli NHRI by pass time. 
AlLbough not zipplicthlc to those daLa sets, cloud- 
iness should hive increased in the NIIRE, MLS, C1IY, 
DCS, and PLD areas later in the afternoon due to 
eastward translation and pi opagation of mountain 
induced thunder .tor.as, with « conospondjng 
dccieasu in mountain cloudiness. 


9. SUMMARY 


Mesoscale processes dominated the 
weal her pattern over the entire study region, 
moi c so over the mountains than over the plains 
areas, lhe smallest and least developed clouds 
were observed at NUKE. More developed, larger 


clouds vc i e ub"ei vi d it t lu 'tin to Ilfl'I'a me is, 
more so in the noil I urn and i nntnl pi lies, 
piobably due to a gitaltr freepu m y of synoptic 
sc lie effects, PID, SPK, and ARM were si I 
sjnll ir in cloud *i/e. The least numhei of clouds 
pre*eni were on the* oenLral portion plains 
averaging 4 (o 6 pi i 10,000 ^vtn' , . More cl oik’s 
v,. ic piesent at Bf.S me! Pill, but, hy far, I he 
gi eate"t number of clouds were present over the 
mountain ate, is. A*’c : lga coverage was spirsu 
everywhere east of Ine noun La in, at this tue 
villi Lhe gi eatest pei cent cloud cover ovi_r Luc 
mrumLaiii sites. This comp ices well with the 
typical diurnal summertime cloud regime ovei 
Colorado. These satellite d Jta are a valuable 
tool for determining cloud popolitioi* and wuh 
Lhe advent of goostne’n onous sitellite dull. 

Cloud systems can K followed f i oin develop non! 
to dissipation at half hour Llire increments. 
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EFFECTS OF CLOUD SIZE AND CLOUD PARTICLES 
ON SATELLITE OBSERVED REFLECTED BRIGHTNESS 


David W. Reynolds, T B. McKee and K. S. Danielson 

Department of Atmospheric Science 
Colorado State University 
Fort Collins , Colorado 


1. INTRODUCTION 

Modeling of tne visible radiation 
scattered from semi-mf inite to finite size clouds 
has been attempted for several years. Hansen 
(1969) and Twomey, Jacobowitz and Howell (1967) 
resolved that for a plane parallel atmosphere 
with semi-mf mite clouds having optical depths 
near 100 (cloud 1 - 1.5 km deep and having liquid 
water contents of 2 gm/nr^) that these clouds 
would be at their maximum brightness. These re- 
sults showed that remote sensing of clouds m the 
mid-visible wavelengths would be of little value 
for clouds of any real vertical extent, and it 
would be impossible to determine anything about 
the clouds microphysical properties. However, 
observational measurements from satellites of 
convective clouds of significant vertical extent, 
(z>2 km), have shown increasing brightness with 
increasing height, Griffith & Woodley (1973) , 
Reynolds & Vonder Haar, (1973), Griffith et. al, , 
(1976). Recent work by Busygin et. al., (1973) and 
McKee and Cox (1973, 1975) have shown that the 
fin’te cloud poses a particular problem for moni- 
toring its brightness due to energy passing through 
the vertical sides of the cloud. Thus their rea- 
soning for this, observed brightness-height change 
is related to the fact that as the cloud grows it 
becomes wider as well as tnicker making side ef- 
fects less important and allowing more light to 
be reflected off the top. However, even this 
theory does not account for some of the observed 
change in brightness for semi-infinite clouds 
whicn will be reported in this article. 

During this past summer, a unique 
data set was obtained during the South Park Cumulus 
Experiement (SPACE) . Data on cloud dimensions and 
cloud position were collected through the use of 
two cameras which simultaneously photographed the 
e<pennental area. Since these stereo photographs 
were taken in a time lapse mode, the growth charac- 
teristics of the clouds could also be deduced 
Cloud microphysical data for the experiment were 
collected by aircraft, radar, and surface obser- 
vations. Three aircraft were used to penetrate 
the cumulus clouds The CSU Areocommander and the 
University of Wyoming Queenaire both are instru- 
mented to study cloud microphysical processes 
Emphasis is placed on collecting data pertamirg 
to cloud particle distributions, IN and CCN con- 
centrations, and liquid water contents as well as 
the state parameters The KOAA/NCAR Explorer 
sailplane, because of its ability to remain m the 
clouds for extended periods of time while collect- 
ing microphysical data, was a particularly unique 


data gathering system which was utilized exten- 
sively in SPACE. The Limon WSR-57 10 cm radar 
and the CSU K-33 10 cm radar provides data on 
position, intensity and motions of clouds with pre- 
cipitation si 2 ed particles. The CHILL dual 10 cm 
and 3 cm wavelength radar with doppler caoaoilities 
was used to determine position and intensity of 
natural clouds and was also used in radar chaff 
tracer exneriments which have yielded substantial 
data on cumulus cloud microphysical processes and 
on cloud morphology. A surface meso-net collected 
data on low level fluxes of heat and moisture A 
surface chase vehicle provided data on precipita- 
tion type and size distribution. Several daily 
air mass soundings were taken to provide data on 
the changes experienced by the non-cloud air mass 
and provided a reference for satellite derived 
cloud-top temperatures. All these data complement, 
support, ana strengthen observations made by the 
satellites of the South Park region 

Along with this data, digital SMS-2 
(Synchronous Meteorological Satellite) visible 
(.5 - .7 pm) (.9 km resolution at SSP) and infrared 
(10 5 - 12.5 pn) (9 km resolution at SSP) radiance 
data were received for this area through the DircKC 
Readout Ground Station at White Sands MiesiL 1 in i 
This allowed us to observe reflect inc" mf-riix >f 
the clouds as well as estimate ti of top hciriits. 

Using the Jlonte Carlo cloud model 
for finite clouds develoood !•'. McKee and Cox, (1973) 
we used different distributions of drop sizes and 
numbers (Deirmendj ran, 1969) along with vnryi ng 
cloud depths and widths to determine how theory 
would predict what the satellite would view from 
its given location in space Results of these runs 
along with satellite observed reflectance will be 
presented m the following sections 

2. MODEL CALCULATIONS 

Fig. 1 shows the two drop size dis- 
tributions used in the model (C.l, C 3), along with 
the average distribution observed during the last 
3 days of SPACF. Note that the two distributions 
used bounded that observed during this period so 
that C.2 results would fall between the curves 
shown m Fig. 2 From Deirnendj lan (1969), we ob- 
tained the 2 phase functions and volume scattering 
coefficients for a water cloud at a wavelength of 
7 urn The phase function labeled C.l is charac- 
terized by a very strong forward scattering peak 
while C.3 is not nearly as strongly peaked. Accu- 
racy of the computations depends on the number of 
photons processed through the computer program 
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Fig. 1. Particle size and number distributions 

available for model input (Deirmendj lan, 
1969) South Park observed distributions 
are shown as the dashed line Note 
curve continues upward showing large 
number of small particles. 


1 2r 



u u i . i i i i i i — i. i 

0 20 90 60 80 

OPTICAL THICKNESS 

Fig 2. Model results showing small variation in 
relative radiance for C.l, C.3 orop size 
aistributions over the optical depths 
shown Note also that a cloud slab 
having a width to depth ratio of 5-1 by 
optical depth 80 is very nearly equal in 
brightness to an infinite layer. 


The accuracy of the relative radiance for a semi- 
irfmire cToud and the top of a cubic cloud is 
indicated by error bars. 

Fig. 2 shows the results from this 
computational study. The first test was to deter- 
mine what effect the tw T o drop size distributions 
would have on the reflected brightness for the 
top of a cubic cloud. To orient the reader, for 


a cloud having a C.2 particle distribution and a 
t = SO, its’ vertical depth would be 1 S km and 
have a liquid water content of 15 gm/m-*. The two 
curves show that over a wide range of optical 
depths, there appears to be i sr.-ill effect due to 
phase functions which were determined from the 
drop size distributions This confirms the early 
work by Twomay, Jacobowitz and Powell that remote- 
ly sensing clouds in the mid-visible wavelengths 
for determining their micropnysieal structure 
seems unlikely. 


The second study carried out was 
to determine the differences m reflectance be- 
tween infinite and cubic clouds Fig. 2 indicates 
the relative radiance for the semi- infinite cloud 
is about 50Z greater than for the top of the cube 
at optical depth of 80 The relative tnfference 
increases for smaller optical depths. The semi- 
mf mice cloud is near a theoretical limit for 
optically thick clouds as the slope of the curve 
is nearing horizontal. The cubic clouds hive a 
theoretical limit identical to the sem-mf mite 
cloud but require a much larger ootical thickness 
to approach this limit. Consequently, the cubic 
cloud would continue to get brighter for increasing 
optical thickness. 

A third feature is indicated by 
calculating the relative radiance for a cloua with 
a width to depth ratio of 5 to 1 while maintaining 
a square top. Theory indicates the radiances are' 
closer to the semi-infinite laver tnan co the cube. 
Observations illustrated in Fig 4 indicate a ratio 
of 10 to 1 needed to approach a iraviraurt brightness. 
Real clouds do not hove flat tops as is the condi- 
tion of the cube. Irregular structure as likely 
to slow the transition to a semi-infinite laver, 
since small irregularities will retain character- 
istics of smaller clouds. 

3. SATLLI.HF OBSERVED REFLI CTED 

BRIGHTNESS 

Digital SMS-2 VISSIR (Visible and 
Infrared Spin Scan Radiometer) data was obtained 
for August 6, 7, and 8th during the initial stages 
of convection over South Park, Colorado. The SMS-2 
satellite was positioned at 115 “V during this per- 
iod giving it a nadir angle of the South Park area 
of 46° and 17° west of a due south view This 
gives S MS a ground resolution in this area of 1.2 
km in the visible. The infrared sensor on board 
allowed cloua top tenperature/heights to be deter- 
mined for clouds with horizontal dimensions greater 
than 12 km. Thus only large clouds could bo viewed 
For those clouds where heights could be determined, 
the maximum heights ranged to 1 5 to 2 hn Fig. 3 
is an SMS-2 view of Che South Park region showing 
the range of cloud sizes observed on August 6 which 
was fairly typical of the clouds observed on the 
following two days although tins day may have been 
a little more active. A comparison was first made 
of satellite measured visible radiance versus satel- 
lite derived cloud top temperature for co-located 
SMS visible - IN digital sectors. Nine separate 
sectors were compared for the three days, all within 
1 hour of local noon. Approximately 250 data points 
per picture were correlated From this, correlation 
coefficients ranged from -.3 to -.65 These were 
fairly low* correlations ard seemed to show that 
visible radiance was not well related to top temper- 
ature/height . The relationship between visible 
radiance and cloud horizontal dimensions was nl'-t' 
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stage until much further in its life cycle 
when the top begins to flatten and glaciate. 
We should add here that many of the larger 
clouds viewed were conglomerates of clouds 
spaced closer together chan 1 2 Km so were 
averaged by the sensor into a large cloud. 
There may De some interaction between the 
cloud sides that we see but this is assumed 
to oe small compared to cloud top reflectance 
properties. Ranges of visible radiance 
measured for these clouds for the wavelength 
interval used agrees well with what theory 
predicts. The discrepancy occurs in the 
sizes of clouds over which this range should 
exist. 


4. CONCLUSIONS 

Several results have been obtained 
through looking at theoretical and observa- 
tional cloud reflectance properties 


Fig. 3. SMS-2 .9 km visible image for 1845Z 
6 August 1975. Area of interest is 
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1) Cannot remotely sense the micro- 
phvsicnl structure and change taking 
place m cumulus clouds m the mid- 
visible spectral region for clouds 
greater than 1.5 km in depth. 

2) Can predict that geometrical factors 
will strongly affect the cloud 
brightness and far outweigh micro- 
physical changes but have not pre- 
dicted these changes for clouds of 
width to depth ratios as large as 
that observed from satellites (10.1). 

3) The theory used can predict the 
range of brightness that should be 
observed, but does specify the width 
to depth ratios in accord with ob- 
servations. 

4) Discrepencies from theor> and obser- 
vations may be due to theory not 
adequately representing the non- 
uniformitv of the cloud top which 

we feel stronglv controls the bright- 
ness observed from satellites. 


Fig. 4 Satellite derived maximum cloud bright- 
ness versus horizontal width for 65 
clouds investigated during this three 
day period. Numoers under points repre- 
sent multiple data points at this loca- 
tion 

investigated Clouds were chosen from the three 
days at three different times around local noon 
and their horizontal dimension and maximum bright- 
ness were determined. Fig. 4 shows the results 
of this study which indicate that the brightness 
does not level off until the width approaches 
20 km For clouds approaching the 2 km neight 
the ratio of width to depth is 10 to 1. This is 
well past' the 5 to 1 ratio suggested by McKee and 
Cox where the brightness changes should level off 
This curve may be demonstrating a lack in the 
model in its handling of the cloud top surface 
features. As was mentioned earlier, the model 
hanales only flat tops of cloucs, while we know- 
that growing cumulus are distorted, rough turrets 
It is felt that due to this cloud top configura- 
tion, that a cloud may not approach the infinite 


Work is progressing on mouifjmg this 
program to handle more of the shape factors involved 
with clouds as well as streamlining die computa- 
tional scheme to reduce the noise problem It ap- 
pears that some reflection if needed on previous 
work which has compared satellite brightness to 
cloud heights and rainfall rates and liquid water 
contents. U*e may better a]>prcciece what these 
apparent brightness changes in observed clouds are 
actually telling us through both Ll.eoretic.il model- 
ing and comparison from ground and aircraft observed 
cloue structure, as well as satellite observations. 
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ABSTRACT 

The relationship between a cumulus clouds’ brightness, horizontal dimension and internal microphy sical 
structure are investigated Cumulus clouds located over the South Park region of Colorado are observed 
by the SMS-2 satellite and their brightness and size are determined Aircraft observations were made m- 
cloud to obtain the drop size distributions and liquid water content (LWC) of the cloud A Monte Carlo 
cloud model is used to imitate the sun-satellite-cloud geometry m an effort to understand the role of cloud 
size and microphvsical structure m affecting cloud brightness 

Results show that for clouds of optical thickness betw een 20 and 60 (i e., LYV C of 0 03 7 gm -3 and 0 1 1 gin' 1 
for a 2 km deep cloud), information about a cloud’s LWC may be obtained through monitoring cloud bright- 
ness for ebuds of umiorm depth and variable width Theoretical results using this Monte Carlo method 
approximate very closely the relative brightness changes of clouds of the size and depth monitored bv the 
SMS-2 satellite for these lew days Theory and observation both conclude that a cloud having a width to 
depth ratio of approximately 10 1 (and constant optical thickness) is nearly reaching it= maximum bright- 
ness Theory predicts that geometric factors affect cloud brightness more than microphy sical changes. 

It is also discussed that the previously reported work on the cloud height-doud brightness relationship 
may indeed be seeing increasing brightness with increasing horizontal size changes (size being related to 
height) with finite small perturbations on top of the growing cloud slowing its approach to maximum 
brightness 


1. Introduction 

Modeling of the visible radiation scattered from 
semi-infinite to finite size clouds has been attempted 
for several years Hansen (1969) and Twomey el al 
(1967) found that for a plane-parallel atmosphere 
with semi-infinite clouds having optical thicknesses 
near 100 (cloud 1-1 5 km deep and having liquid 
water contents of 0 2 g m -3 ) that these clouds are 
within 20% of their maximum brightness These 
results showed that remote sensing o! clouds in the 
mid- visible wavelengths would be of little value for 
clouds of any real vertical extent, and it would be 
impossible to determine anything about the micro- 
physical properties of the clouds However, observa- 
tional measurements from satellites of convective 
clouds of significant vertical extent ()>2 km) haxe 
shown increasing brightness with increasing height 
(Griffith and Woodley, 1973, Reynolds and Vonder 
Haar, 1973) Recent work by Busygin el al (1973) 
and McKee and Cox (1974, 1976) has shown that 
the finite cloud poses a particular problem for moni- 
toring its brightness due to energy passing through 
the vertical sides of the cloud Thus their reasoning 


for this observed brightness-height change is related 
to the fact that as the cloud grows it becomes wider 
as well as thicker, making side effects less important 
and allowing more light to be reflected off the top. 
To further understand cloud brightness changes, 
a small study w*as performed using measured satellite 
brightness data and comparing these to expected 
theoretical results 

The South Park Area Cumulus Experiment (SPACE) 
collected microphysical data on cumulus clouds m 
South Park, Colo , during the summer of 1975 The 
data set included cloud drop size distributions made 
with an electrostatic disdrometer (Kelly and Millen, 
1960) mounted on the NCAR/XOAA Explorer sail- 
plane This instrument measures cloud droplets from 
4 to 16 juni in 1 5 jam increments and has a ninth 
channel for cloud droplets >19 gm Liquid water 
contents were simultaneously measured with a Johnson- 
Witliams liquid water content meter mounted on the 
Explorer 

Along wuth these data, digital SMS-2 (Synchronous 
Meteorological Satellite) visible (0 5-0 7 jam) (0 9 km 
resolution at SSP) data were received for this area 
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through, the Direct Readout Ground Station at White 
Sands Missile Range. This allowed us to observe 
visible brightness of the clouds as well as estimate 
their top heights from the infrared readout using an 
appropriate temperature sounding. 

A' Monte Carlo cloud model for finite clouds de- 
veloped' by McKee and Cox (1974) was run using 
different distributions of drop sizes and numbers 
(Deirmendjian. 1969), while varying the cloud depth 
and width to determine how theory would predict 
what the satellite would view from its given location 
in space Comparison of these results to the satellite 
observed reflectances will be presented m the fol- 
lowing sections 

2. Model calculations 

Fig. 1 graphically compares the two theoretical drop 
size distributions (C.l and C 3) used in the numerical 
simulation of satellite measured visible brightness, 
with the average drop size distribution observed on 
6-8 August of the 1975 SPACE field program The 
specific values but not the shape of the depicted 
theoretical distributions are dependent on the liquid 
water content (LWC) of the clouds. C 1 and C 3 m 
Fig 1 are adjusted to a LWC of 0 055 g m~ 3 which 
was the average observed LWC of the three days 
under study. The daily average LWC values ranged 



Fig 1. Particle size and number distributions used for model 
input (Deirmendjian, 1969) Average South Park observed dis- 
tr.buuons for three case study days are also shown 



OPTICAL THICKNESS 


Fig 2 Model results showing small variation in relative 
radiance for C 1 and C 3 drop size distributions over the optical 
thicknesses shown Mote also that a cloud slab having a width 
to depth ratio of 5 1 by optical thickness 80 is very nearly 
equal in brightness to an infinite layer 

from 0 047 g m -3 on 8 August to 0 07 g m -3 on 7 August 
This was lower then one might expect for cumulus 
clouds but we feel these numbers are representative 
The observed drop size distribution is also plotted 
on Fig. 1, primarily to allow a comparison of the 
shape of the observed distributions with the theoretical 
distributions utilized in the numerical simulation The 
drop size distribution observed with the electrostatic 
disdrometer has several problems which are discussed 
by Dye (1976) One prominent difficulty is that the 
LWC integrated from the electrostatic disdrometer 
is typically half the LWC observed with a Johnson- 
Wiliiams (J-W) liquid water content meter The 
average LWC as observed by the J-W on 6-8 August 
was ~0 18 g m -3 Second, the distributions are skew r ed 
toward smaller drops by the electrostatic disdrometer. 

The observed distributions were measured along the 
entire vertical extent of the clouds For the three days 
studied, the distributions were generally monomodal 
and showed good colloidal stability The sailplane did 
observe ice particles in the clouds but ice particles 
are not considered m the calculations presented 
Ex aluation of both theoretical distributions C 1 
and C 3 were processed through the model because 
neither distribution exactly approximates the observed 
distribution and the sensitivity of the simulation may 
be sensitive to the shape of the distribution The 
theoretical droplet distribution C 1 and C 3 and cor- 
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responding phase functions and volume scattering 
coefficients (at 0.7 «im wavelength) for a water cloud 
were obtained from Deirmendjian (1969). The phase 
function labeled C.l is characterized by a very strong 
forwardscattering peak, while C.3 is not nearly as 
strongly peaked. 

The relative radiance from the cloud top for the 
satellite geometry is shown in Fig. 2 for a cloud layer 
of specified thickness and of semi-infinite horizontal 
extent (C.3), a cloud with width 5 times the depth 
(C.3), and a cubic cloud with two different particle 
size distributions. Relative radiance is presented as a 
function of optical thickness defined by 


-r 


0ds, 


( 1 ) 


in which r is optical thickness, s the distance and J 
the volume scattering coefficient. The volume scat- 
tering coefficient is determined by the drop size dis- 
tribution and the liquid water content. For a given 
drop size distribution the scattering coefficient varies 
linearly with liquid water content. Consequently, the 
optical thickness varies with the three parameters of 
the drop size distribution, with liquid water content 
and with geometric distance. Accuracy of the com- 
putations depends on the number of photons processed 
through the computer program. The accuracy of the 
relative radiance for the semi-infinite cloud and the 
top of the cubic cloud are shown by error bars for 
one standard deviation. For reference, a cloud with 
a C.l particle distribution, an optical thickness of 80 
and a liquid water content of 0.15 g m -5 has a vertical 
depth of 2.0 km. The two curves for the cubic cloud 
indicate only a small change in relative radiance due 
to particle size distributions for a large range of 
optical thickness. 

In contrast, the changes in relative radiance for 
the semi-infinite cloud is about 50% greater than for 
the top of the cube at an optical thickness of 80. 
The relative difference increases for smaller optical 
thickness. The semi-infinite cloud is approaching a 
theoretical limit for optically thick clouds as the slope 
of the curve is approaching zero. The cubic clouds 
have a theoretical limit identical to the semi-infinite 
cloud but require a much larger optical thickness to 
approach this limit. Consequently, the cubic cloud 
would continue to get brighter for increasing optical 
thickness. 

Another feature of shape is indicated by calculating 
the relative radiance for a cloud with a width to 
depth ratio of 5: 1 while maintaining a square top. 
Th- ory indicates the radiances are closer to the semi- 
infinite layer than to the cube. 

3. Satellite observed reflected brightness 

Digital SMS-2 YISSR (Visible and Infrared Spin 
Scan Radiometer) data were obtained for o 8 August 


1975 during the initial stages of convection over 
South Park. Colo. The SMS-2 satellite was positioned 
at 115°W during this period resulting in a satellite 
zenith angle for the South Park area of 46° and a 
relative azimuth angle of 17° west of a due south 
view. Since all observations were made within one 
hour of local noon a fixed solar zenith angle of 23° 
was used for model calculations along with the satel- 
lite viewing angles. The ground resolution of the SMS 
visible data in the South Park area is 1.2 km. The 
infrared sensor on board allowed cloud top tempera- 
ture and heights to be determined for clouds with 
horizontal dimensions greater than 11 km. Thus only- 
large clouds could be viewed in the IR. For those 
clouds where heights could be determined, the maxi- 
mum heights ranged from 1.5 to 2 km above cloud 
base (cloud base determined by surface observations). 
Fig. 3 is an SMS-2 view of the South Park region 
showing the range of cloud sizes observed on 6 August 
which was fairly typical of the clouds observed on 
the following two days although this day may have 
been a little more active. A comparison was first 
made of satellite measured visible brightness versus 
satellite derived cloud top temperature for co-located 
SMS visible-IR digital sectors. Nine separate sectors 
were compared for the three days. Approximately 
250 data points per section were correlated including 
some noncloudy regions. From these correlations the 
nine coefficients ranged from —0.3 to —0.65. These 
were fairly low correlations but for this sample size, 
still significant at the 1% level. Thus, although noisy, 
some relationship exists between visible brightness 
and cloud height. In an attempt to further explain 
this, we investigated the relationship between visible 
brightness and cloud horizontal dimensions. It is felt 
that as a cloud increases in height, it also increases 
in width with the width increase dominating the 


t 18:45 06AU75 







Fig. 3. SMS-2 1.2 km resolution visible image for 1845 GMT 
6 August 1975. The South Park area of Colorado and sur- 
rounding areas is shown in the outlined box. Scale : 1 mm = 6.2 km. 
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Fig. 4. Digital SMS-2 1.2 km resolution visible display of clouds over the South Park region 
of Colorado on 8 August at 1817 GMT. Using a thresholding technique only clouds are shown. The size 
and maximum brightness of the clouds chosen at random for this study were obtained off printouts such 
as this. Latitude in degrees, minutes, seconds is given on the right along with the character representation 
per count interval used. 


brightness change. We hope to show this in the fol- 
lowing results. 

Clouds were chosen from the three days at three 
different times around local noon and their horizontal 
dimension and maximum brightness [here maximum 
brightness refers to the brightest pixel measured for 
each individual cloud as taken from the digital display 
(see Fig. 4)3 were determined. Fig. 5 shows the results 
of this study which indicate that the brightness does 
not level off until the width approaches 20 km. For 
clouds approaching a height of 2 km, the ratio of 
width to depth is 10: 1. 

The solid lines in Fig. 5 are the theoretical relative 
radiances (brightnesses) from clouds with fixed thick- 
ness of 2 km and a variable width which ranges from 
a cube at width 2 km to a slab 40 km wide. Maximum 
radiance from theory has been scaled to coincide 
with the highest satellite maximum (count of 57) 
brightness so that relative changes could be compared 
with theory’. (Note the convergence of the two lines 
at optical depth 40.) Two different liquid water con- 
tents are represented to illustrate the effect on cloud 
brightness. Optical thickness 60 is derived from a 
liquid water content of 0.11 g m -3 and optical thickness 



Fig. 5. Satellite derived maximum brightness versus horizontal 
width for 05 clouds investigated during the three-dav period. 
Numbers under points represent multiple data points. Also 
shown arc the theoretical results dor satellite geometry i using 
the C l distribution for a cloud of fixed depth 2 km) and vari- 
able width. The r of 20 and 60 correspond to LWC’s of 0.037 
and 0.1 1 g m -1 , respectively. 
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20 from 0 037 g m“ 3 The theory predicts both the 
magnitude of the brightness change and the shape 
of the change rather well and indicates a sensitivity 
to the liquid water content The sensitivity to liquid 
water content is actually restricted to clouds with 
optical thickness <60 Results of Fig. 5 indicate a 
distinct possibility of inferring liquid water content 
of cumulus clouds from satellite measurements. More 
work is needed to separate eSects of particle size 
distribution and liquid water content on the scat- 
tered radiation. 

4. Conclusions 

Several results have been obtained through looking 
at theoretical and observational cloud reflectance 
properties • 

1) For clouds of optical thickness between 20 and 60, 
information about a cloud’s liquid water content may 
be obtained through monitoring cloud brightness 
changes for clouds of uniform depth and variable 
width. Beyond this point, geometrical factors dominate 

2) Theoretical results using this Monte Carlo 
method approximate very closely the relative bright- 
ness changes of clouds of the size and depth monitored 
by the satellite for these few days. 

3) Theory and observations both conclude that a 
cloud having a width to depth ratio of approximately 
10 1 is nearly reaching its maximum brightness for 
a specified optical thickness 

4) Theory predicts that geometrical factors will 
strongly affect the cloud brightness and far outweigh 
microphysical changes Thus, finite perturbations on 
top of a large cloud (»2 km in depth) may account for 
these types of clouds increasing in brightness past the 
theoretical limits show n here 

Some initial work has been done m adding finite 
shapes onto semi-infinite clouds in the Monte Carlo 
program These show that the finite perturbations do 
decrease the rate at which a cloud will reach maximum 
brightness Thus theory lends evidence to our final 
conclusion that a cloud may continue to increase in 
brightness past previously indicated size limits if the 
top of the cloud has distorted finite shape factors, i e , 


growing cumuius clouds. As satellite resolution con- 
tinues to increase, this finite cloud problem must be 
dealt with to a greater degree if satellite reflected bright- 
nesses are to be interpreted correctly with respect to 
cloud heights, rainfall rates and liquid water contents 
A better appreciation of these apparent brightness 
changes in observed clouds will be revealed through 
both theoretical modeling and comparisons to ground, 
aircraft and satellite observations 
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INITIAL CONVECTION DURING SPACE 1 , 19?5 

by 

Daniel Breed 


INTRODUCTION 

The reason for locating a cloud physics experiment m South 
Park was that the early formation of small, isolated cumuli 
over the Mosquito Range and their subsequent growth and move- 
ment eastward over South Park created a natural laboratory for 
surface observations and aircraft operations. This study at- 
tempts to determine the time and location of initial convec- 
tion ("first puff") m the South Park area using SMS-2 satel- 
lite photos for the 19?5 season ( July?-August ?) . 'ihe formation 

i -* 

of these cumuli are attributed to the differential heating of 
the mountains and to the adequate summer moisture supply m the 
park. Location of primary heat source areas are attempted, 
and one case study of thunderstorm organization and movement 
is presented as an introduction to various studies that are 
feasible using this data. 

The data base for this study was SMS-2 visible photograph: 

2 

scanning a particular sector every half hour. One mile rtboj w - 
tion photos (one data point per square mile) for the KB8 sector 
(central North America) were available for the entire period 
and the visible pictures usually began by i4:45 GMT, For part 
of the period, one half mile resolution photos were used for 

■^South Park Area Cumulus Experiment 

2 

There were a few SMS-1 pictures also. 
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more accurate location of first puffs when the scanned sector 
included South Park (e.g. 32N102W, printed as the last row of 
numbers on the upper right portion of the photo) . Table 1 lists 
the days, times, and resolution of the available visible satel- 
lite photos that include the South Park area. A 'log book of 
field observations and first puff forecasts from the morning 
rawinsondes were used for verification and comparison of the 
satellite-determined initial times. The 07:00 MDT sounding was 
usually followed by an afternoon sounding between 12:00 and 
l4: 00 MDT which would be useful in following the growth of 
these cumuli. Other data available for use with the satellite 
photos include cloud stereo-photography (terretrial photogram- 
metry) , a doppler radar (with hail signal) , a mesoscale network 
measuring pressure, temperature, relative humidity, and winds, 
and cloud physics aircraft penetrations. Development of precip- 
itating cumuli, distinction between hail and non-hail convective 
systems, wind shear and flow fields are a few satellite projects 
that could use the surface observations as ground truth. 

PROCEDURE 
A. Set-up 

The 0 d 3 was used to magnify the pictures so that visual 
analysis could be made more easily. The digitizer and compu- 
ter were not used in this study, because the abundance of vis- 
ual landmarks in South Park and the complications of setting 
up orientation overlays to account for the wobble of the SMS 
satellite justified the procedure of visual locations of ' first 
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DAYS 

TIDIES (GMT) 

RESOLUTION 

DAYS 

TIMES (GMT) 

RESOLUTION 

7/7 

14:45-23 :45 

1 mile 

7/24 

l4: 46-23: 45 

1 mile 




n 

14:46-00:15 

5 mile 

7/8 

14:45-23:45 

1 mile 







7/25 

14:00-23:45 

1 mile 

7/9 

14:47-23:45 

1 mile 

it 

14:46-00:15 

| mile 

11 

20:45-21:15 

|r mile 


* 





7/26 

14:45-23:45 

1 mile 

7/10 

15:15-23:45 

1 mile 

U 

14:45-00:15 

mile 

* 

7/1 1 

14:45-23:45 

1 mile 

7/27 

15:15-23:45 

1 mile 




tl 

15 : 1 5 - 00: 15 

\ mile 

7/12 

14:45-23:45 

1 mile 

) 




\ 


7/28 

14:45-23:45 

1 mile 

7/13 

. 14:45-23:45 

1 mile 

>1 

14:46-00:15 

ts mile 

7/1 4 

13:30-23:45 

1 mile 

7/29 

13:30-23:45 

1 mile 

ir 

13 : 30 - 00:00 

mile 

It 

14:45-00:15 

mile 

V/15 

15:45-23:45 

1 mile 

7/30 

14:46-23:15 

1 mile 

II 

13 : 00 - 23:30 

4 mile 







7/31 

14:45-23:15 

1 mile 

7/1 6 

14:45-23:45 

1 mile 




1 ! 

17:30-2301 

\ mile 

8/1 

— 

1 mile 

7/17 

12:45,14:45-23: 

45 1 mile 

8/2 

— , — 

1 mile 

7/18 

14:45-23:45 

1 mile 

8/3 

14:45-23:15 

1 mile 

7/19 

15:15-23:45 

1 mile 

8/4 

1^:45-23:45 

1 mile 

7/20 

14:45-23; 45 

1 mile 

8/5 

14:46-23:45 

1 mile 

1 ! 

l6:45-23;45 

| mile 

tl 

14 : 45 - 00: 15 

4 mile 

7/21 

14:45-23:45 

1 mile 

8/6 

15:15-21:45 

1 mile 

tl 

15:15-00:15 

2 mile 

II 

15 :15-22:16 

5 mile 

7/22 

14:45-23:45 

1 mile 

8/7 

14:45-21:46 

1 mile 

K 

12:45,14:45-00: 

15 i mile 

It 

1 * 1 : 45-22:15 

h mile 

7/23 

14:45-23:45 

1 mile 




tl 

14:45-23:45 

2 - mile 





TABLE 1 

Available visible satellite photos that include South Park (1975). 
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puffs, "he equipment could be used in determining cloud top 
brightness or other characteristics in a hail study. Exten- 
sion tubes were tested to increase the magnification of the 
picture onto the TV screen. Eut, the contrast became so poor 
and the quality of the picture became so coarse that the trade- 
off between magnif ication and contrast resulted m the use of 
the camera lens with no extension tubes. The picture to be 
studied was put on the light table under the camera and flat- 
tened by the glass cover. It was then adjusted so that South 
Park was centered on the TV screen. The f-stop was adjusted 
for the best contrast and the focus was adjusted for the sharp- 
est picture. This was the initial set-up f.or examining the 
satellite photos. 

B. Initial Time 

The time of the first puff was determined by examining the 

-a 

early morning photographs by eye or with a magnifying glass, 
and then choosing two pictures, one with the first visible con- 
vective element and one a half hour earlier. These two photos 
were examined on the 0 D- , first checking to see if there was 
any trace of the cloud In the earlier one, and then using the 
photograph that contained the first visible cumulus to deter- 
mine location. Accuracy of the Initial cloud time is £ 15 min- 
utes, using the average time of the two photos (no clouds and 
first puff). Also, since the resolution is 7 mile at best, 
this time would tend to be an overestimate because of the growth 
time necessary to develop a cumulus large enough to see on the 
satellite picture. There were a few days when high clouds and 



raesoscale systems obscured th?- Louth Park area, and rough esti- 
mates were made when the first convective elements could be 
detected on the pictures. 

C. Location 

The first puff on any particular day was located as accur- 
ately as possible from surrounding landmarks. Mountain valleys, 
peaks, reservoirs, and ridges gave a gray scale contrast on the 
SMS-2 photos that was readily comparable to an ERTS composite 
photograph of Colorado. In order to locate' the South Park 
area and orient the shapes of the landmarks as they change 
with the wobble of the satellite orbit, the ERTS photo and tne 
SMS-2 photo with no clouds (and a clear shot at South Park) 
were kept handy for reference. Although detailed topographic 
maps of the area were available, the size of the National Forest 
maps worked best for the accuracy of this study. 

The Arkansas River valley to the west of the Mosquito Range 
and the westernmost edge of South Park are separated by 10 
miles at the narrowest part. About one fif/th of that separ- „ 
at ion was detectable over another fifth. In low contrast areas, 
this detection limit would decrease to about one fourth. So, 
although the resolution was at least one mile, visual detection 
was limited to 2 to 2.5 miles. Allowing for other errors, a 3 
mile grid (9 square miles) was set-up on a Pike National Forest 
map of the South Park area (Figure l). A San Isabel National 
Forest map would have been useful as it includes the area west 
of the Mosquito Range, but it was not immediately available. 
Figure 2 shows the gridded area and some landmarks on a trace 
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from the EETS photo. 

The grids are numbered with four digits, the first two 
identifying the \ 5 rows (01 the northernmost-, 15 the south- 
ernmost) and the last two identifying the 12 columns {01 the 
westernmost, 12 the easternmost). So, m Fig. 1, the bold 
number 5 would be in 0809, the bold 3 m 0409, and the bold 
8 in 1112. (The bold numbers happen to be the locations of the 
mesoscale weather stations.) The locations of the initial 
cumulus clouds were identified using this grid. 

RESULTS 

Initial times and locations are assembled m Table 2 and 
the X's in Fig. 1 show the grid frequencies of the locations. 

The general comments in Table 2 describe other convective areas, 
areas of further or rapid development, the quality of the photos, 
and the days when the park was obscured or partially obscured 
by early clouds. If a line of clouds had formed, the brightest 
cloud was considered to be the first. Two or more grid loca- 
tions on a single day indicate that there were tv/o or more 
clouds of equal brightness on the examined photo. Although the 
west side of the park was of the most interest, there were times 
when convection was clearly dominant m the north or east areas. 
The comments should clarify the differing situations. 

The morning sounding was used to predict the time of the 
first puff. A knowledge of the normal morni] g heating rate in 
South Park, a judgement of moisture adequacy, and a synoptic 
forecast were used to determine when the surface inversion 
would be broken and how fast and deep the clouds would gro w. 




TABLE 

? lime and 

OF foot 1 PAGJS 18 

W0Ti QUAirm 

location of first puTfs (1975)* 

DATE 

TIME (GMT) 

GRID 

LOCATION 

COMMENTS 

7/7 

l 6 : 00 

0405 

Cu line along Nos. Range, also Cu north 
Kenosha Range. 

7/8 

15:30 

0804 

Many landmarks obscured by early clouds 

7/9 

l 6:00 

0303,0703 

At 15:45 a Cu shoos over 02Q1, cloudy 
over Leadville valley 

7/10 

15:00 

0607 

Partially obscured by early clouds so 
first convective element reported 

7/11 

l4: 30 

not 

possible 

Obscured by early clouds, convection 
starts on 14:45 photo 

7/12 

16:30 

(0701,0801 ) 
(0504,0405) 

Two convective elements reported, mostly 
obscured by early clouds, photo is fuzzy 

7/13 

15:30 

0902 

One cloud at 1504 on 14:45 photo that 
dissipates 

7/14 

15:15-15: 

: 30 0604 

Missing 15:15 KB 8 photo 

7/15 

15:15 

0704 


7/1 6 

16:00 

0704 

Some Cu north of Kenoshas at 15:45 

7/17 

14 : 30 

0803 

Partially obscured by early clouds, Cu 
form over Kenoshas at 15:15 

7/18 

15:00 

0306 

Possibly Cu on l4:45 photo but very light 

7/19 

15:30 

0105,0903, 

0703 

Cu also over south Kenoshas and Tarryall 
Reservoir (0612) 

7/20 

case study day early convection obscured by clouds 

7/21 

0 0 
0 0 ^ 

Oavpv 

01123 

1104 

(of marginal interest) 
Over Nos. Range by 16 : 15 

7/22 

16: 00 

0506 

Many early clouds, Cu also over lios. 
Range by 16:15 

7/23 

14: 30 

0903 

14:45 photo poor quality, 15:15 shows 
rapid development 

7/24 

15:00 

1103,0811 

A line develops just south of Kenoshas 

7/25 

15:30 

1104 

A couple less bright Cu near Sheep Kt. 
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TABLE 2 

(con. ) 


DATE 

TIME (GMT) 

GRID 

LOCATION 

COMMENTS 

7/26 

15:30 

07l2 

l 6 :l 5 shows rapid development m this 
area, and development beginning north 
of the park 

7/27 

16:30 

1103 

This is a lone cloud, none over park 
until 21:15 

7/28 

15:30 

0504,0403 


7/29 

15:00 

0404 

Cloudy over Leadville area, colic bo — 
r.m~ ir.~ in-gay- sout-hw-e-e-t 

7/30 

15 : 00 

0403 

Cells beginning m far southwest 

7/31 

15:00 

0605 

Convection reigns by 16:45 

8/1 

not done 



8/2 

not done 



8/3 

20:30 

1101 ( very 

faint) Best convection is south of Buena 

Vista and over Eenoshas at 22:45 

8/4 

16:00 

0502 

Cu over Thirty-Nine Mile fit. at 15:45, 
major cells develop far south end 

8/5 

17:30 

0311,0110 

Cu north of Kenosha Pass at 17:15, con- 
vection started at north end, over Nos. 
Range by 18:15 

8/6 

16 : 00 

0407,0505 

16:15 photo is fuzzy 

8/7 

16: 00 

0704 

Two Cu just north of Copper Kt. 
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This was operationally important in determining 1 ) whether 
or not the day was an experimental day, 2) the best time for 
NCAH's sailplane to take-off, 3 ) the possibility , of an early 

Arcp 

morning chaffy for the doppler radar to follow, and the time 
to start the photogramme try cameras. Table 3 contains the 
first puff forecasts from the morning soundings, observations 
at the time of launch, and any additional comments from the 
daily log that were pertinent to this study. Satellite-deter- 
mined times in NDT are included for comparison. 

CONCLUSIONS AND DISCUSSION 

At the beginning of this paper, differential heating was 
sited as a factor influencing convection during SPACE. These 
conclusions continue alonj that line. All but five days have 
initial convection along the eastern slop>es of the hosquito 
Range or southeast of the Continental Divide (northcentral to 
northwest border of South Park) . It is obvious that the early 
morning sun strikes these areas first, causing an initial 
temperature difference to form between these slopes and the 
environment. This is also supported by the fact that convec- 
tion was prevalent over the Kenosha Range on most days. This 
range has a northeast facing slope with no ridges to the east, 
therefore enhancing its ability to differentially heat since 
the sunrise was a little north of east during the period stu- 
died . 

Grid 0?03 is reported twice and grid 0704 is reported 
three times. This hot spot Is characterized by the proximity 



TABLE-3 Initial times and comments from radiosonde data (in hDT) 


SAT. 


DATE 

RAOB FORECAST 

OBSERVATIONS 

LOG BOOK COMhENTS 

TIhE 

7/7 

Not Available (N/A) 

• 

N/A 

10:00 

7/8 

9' 30, good moisture and 
instability 

7:45, hazy, stratus 
over north mts. 

N/A 

9:30 

7/9 

10:00, very good instab. 

7:15, clear, stratus 

N/A 

10:00 

7/10 

9*15, good moisture 

7:05, fog on foothills 
hazy, low clouds 

10:00, clouds popping 
all over park, still 
gen'l cloudiness 

9:00 

7/11 

10:00, upper inversion 

8:20, fair, cirrus 
and middle clouds 

Cu' deck til 10: 30 
l4:30 small isolated Cu 

8:30 

7/12 

9:30, low CCL 

7:00, fair, scattered 
fog 

8:00 low clouds 
12:00 scattered Cu 

10:30 

7/13 

9:00, two upper level 
inversions 

7:10, fair and cold, some 
fog,Cu pancakus on mts. 

Rain at 15:30 

9:30 

7/14 

9:00 

7:15, fair, clear 

9:00 first small Cu 9: 

15-9:30 

7/15 

9: 15 

6:45, fair, some cirrus 
and cirrocumulus 

10:50 Cu starting to 
close in on park 

9:15 

7/1 6 

9:30, pretty unstable 

7:55, some haze, cirrus, 
Cu pancakus 

N/A 

10:00 

7/17 

8:45, good instab, til 
inversion at l 6K 1 

7:20, fair, cirrus, some 
Cu over west mt s . 

Convection not material- 
izing thru 12:00 

8:30 

7/18 

9:15 

6:50, fair , clear , some 
Cu over mts. 

15:00 fair w:c Cu 

9:00 

O 

7/19 

8:45, upper inversion at 
16.5K 1 (possibly unreported 
high clouds) 

N/A 

Clouds are small 

9:30 



TABLE 3 (con.) 


SAT. 


DATE 

RAOB FORECAST 

OBSERVATIONS 

LOG BOOK COMMENTS 

TINE 

7/ 20 

9: 45, not good instab, til 
l4. 5 K' 

7: 10, fair, middle deck 
of Cu, cirrus 

9:00 cameras on 
11:00 small showers 

— 

7/21 

10:00 

6:35» fair, clear 

9:00 cameras on 


9:30 

7/22 

9:30 

N/A 

12:00 clouds of 
dium size 

me- 

10:00 

7/23 

9:00, two upper inversions 

6:35, fair, clear 

N/A 


8:30 

7/24 

9:15 

6: 35, fair, clear 

11; 00 sailplane 

up 

9:00 

7/25 

9:45 

6:50, fair, stratus 

9:4l cameras on 
1 0:30 good iso. 

Cu 

9:30 

7/26 

9:45 

N/A 

no comments 

© © 

9:30 

7/27 

10:15 

6:45, fair, clear 

11 


10:30 

7/28 

10:00, good instab. 

7 : 00 | fair , stratoCu ,Ci 

11 

& > 

9:30 

7/29 

9:^5 

7:00, stratoCu 

no comments 

JO 

9:00 

7/30 

radiosonde was late and puffs 

had already formed by 9: 

00 

& £ 

*' 9:00 

7/31 

9:30 

6:45, fair 

sonde very stable 

9:00 

8/l 

9:15 

N/A 

7 : 30 chaff drop 
picked up by Cu 


— 

8/2 

9:00 

N/A 

sonde stable and 
13:40 fair wx Cu 

ary 

-- — 

8/3 

9:00 

6:55, sunny , clear 

N/A 


14:30 

8/4 

9:15 

8:30, balmy 

time (?)-fair wx Cu 

10:00 





TABLE 3 (con.) 



DATE 

RAOB FORECASTS 

OBSERVATIONS 

LOG BOOK COMMENTS 

SAT. 

TIME 

8/5 

9:15, good instability 

6 : 50 , fair, clear 

AH very dry, few 
small Cu 

11:30 

8/6 

9:15, two upper level 
inversions 

7:00, fair, some 
cirrcCu 

12:00 sailplane up 

10:00 

8/7 

9:00, canning: inversion 
at 19 K'*‘ 

N/A 

8:20 first puff 

10:00 
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of Sheep Mountain and the valley containing Four Mile Creek 
(site of the Mountain Met. course two summers ago). Another 
frequently sited location is the 1103, 1104 area (twice each). 
Buffalo Peaks are to the southwest, and Buffalo Meadows and 
Long Park surround this area. The area north of Fairplay is 
rather active also. Although there is no favored grid location, 
this part of the park is a broad sloping area up to Pit. Silver- 
heels, Palmer Peak, and Little Baldy Pit. Also, there are many 
gulches that could act as pockets for moisture. The mechanism 
most probat'le for initial convection is the heating of the ele- 
vated surfaces to a higher temperature than the surrounding air 
causing convective currents. This heating breaks through the 
inversion more quickly being at a higher altitude. Also, some 
mixing must occur in the lower levels prior to the break, because 
the CCL of the morning sounding was consistently lower but al- 


ways close to the true cloud base. The mixing does not have to 
be very intense, as lower levels with higher moisture content 
are close to these hot spots. Hence, convection occurs rapidly 
at these locations, transporting moisture effectively to create 
the first cumulus clouds in the South Park area. 

Some personal observations follow that support these find- 
ings. The area southwest of V/eston Pass- appeared as a valley 
between the Buffalo Peaks and the Mosquito Range from the radar 
site (see Fig. 2). This was often the site of the first radar 
returns with the showers often raying northeast over the radar 
site. The three or so occasions when hail occured directly 
above the radar appeared to come from this region also. It is 



interesting to note that the four locations where data was 
taken by the chase van (equipped for raindrop distribution 
measurements and electric field changes) were situated just 
downwind from these hot spots (e.g. l80°-270°). The chase van 
usually locate^cis close to the strongest echo as possible. 

In summary, the first puffs appeared above slopes with the 
longest time of intense sunlight, and followed fairly closely 
the predicted times from the morning sounding. Also, it is 
evident that the strongest echo-producing cells may grow from 
these first puffs, or at least take advantage of these hot 
spots. While these results are not earth-shattering, many 
interesting studies could be done with the data presented. 

As an introduction along these lines, a- descriptive case study 
of July 20th is presented. 

CASE STUDY: JULY 20, 19 75 

This day was chosen for many reasons. The first was the 
fact that initial convection could not be looked at on this 
day because of cloud cover, but the day looked too interesting 
to pass up. Secondly, half-mile resolution pictures that in- 
cluded South Park were available and of good quality. Another 
important reason was that initial cursory examination revealed 
some real intense storm cells developing during the day, and 
upon reviewing the radar log, the 20th turned out to be a day 
•of many hail events. 

The sequence of pictures to be discussed begins at 16:45 GMT 
and ends at 21:49 GMT. It appears that the Colorado-New 'Mexico 
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area is under some general dynamic lifting because of the fre- 
quency of thunderstorm activity. This is not immediately ap- 
parent from the 500 mb pattern. However the map does show 
light and variable winds in the South Park area, accounting 
for the variability of direction and of speed wath Which the 
reported storm cells move. Early storms are evident along the 
southern Colorado border and in southeastern Colorado, and con- 
vective cells have formed over the Breckenridge-Copper fit . 
area. These cells organize^ move eastward as a storm center. 

It loses definition during 18:15 and 18:45, but becomes well 
defined as an intense center in the 19^15 photo { 19*15 is miss- 
ing) . It eventually organizes with a larger system to the 
north. 

During this time, another strong cell develops m the south- 
ern Arkansas River valley near Salida. This one grows and 
moves northeasterly, organizing with a line of convective cells 
probably originating over the Collegiate Range. This line 
moves northeasterly over the Kosquito Range, and possibly organ- 
izes with the first storm between 18:45 and 19:45, making it so 
intense and well defined. General cloudiness and some small 
cumuli move in behind this system. However, the area around 
Buena Vista has been relatively clear for an hour or so, in- 
creasing the instability in that region. The small cumuli, 
upon moving over this area ( 19 : 45-20 : 1 5) grow rapidly as sep- 
arate entities. By 20:45, some rather intense cells have de- 

/ 

veloped, and continue to grow and organize to some extent by 
21:49. Although South Park is obscured after this time, the 
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growth and extent of some of the cells over the park become 
very impressive. The radar data available on this day could 
serve as groung triith for investigating storm cells with a 
possible hail signature. (Note: Around 16 : 00-1.7: 00 NDT, 
pea and marble sized hail accumulated in large amounts 
(drifts) at the radar site.) 
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1. , INTRODUCTION 

The importance of extra -storsr scale 
mesoscale systems to the genesis and propagation 
of ru.r.ulimous systems is only recently becoming 
fully appreciated. A number of workers (Sktlkcs 
and Riehl, 1964, ’latsanoto et al, 1967; Lavoie, 
197?; and linurral’tar, 1973) have found tnat 
cumulus activity is strongly correlated vith 
masosccJe and synoptic scale convergence rones. 
Chang and Orville (1973) found that a two- 
dimensional cunt' las rodel responded much more 
vigorously and penetrated to a consioerable 
depth ’ -hen large-scale convergence was imposed 
on trie boundary of the cumulus model. Further- 
more, Cotton, Pielke, and Gannon (1976) showed 
that e one-dinersa onal time-dependent cumulus 
node! developed a significantly deeper, Iciger 
lasting cloud when the initial sounding was 
replaced by theoretical soundings predicted with 
the Pielke (1974) sesoscale model in the vicinity 
of observed cumulus models. 

In addition, Cotton at al (1976) 
found the sea breeze mesoscale system modified 
the cumulus scale environment by- 


penetrating do.mdrafts, that the extra-storm 
scale mesoscale convergence fields are no longer 
of prime importance. Recently, Pielke (1974) 
developed a fully three-dimensional tesoscale 
model of the Florida sea freeze circulation. Tt>e 
sea breeze was initiated and evolved in resoorse 
to tne boundary layer ilu'es of heat and 'nomes '-up- 
end interacted with the large scale flow The 
effects of deep, precipitating cumuli < nd 
cumulonimbi on the mesoscale circulation were not 
considered m the simulation. In spite of tnis 
neglect, Pielke found that on synopticalj-y 
undisturbed days with generally weak flow t (trough 
the deptn of tne tropospnere, the p-edicLcu 
patterns of convergence agreed favorably with the 
observed locations of clovis and showers througn- 
out an extensive portion of the day. Contrary to 
the generally accepted view point mentioned above, 
the agreGm°nt between the locations of predicted 
convergence zones and locations of echo patterns 
of extensive tnunderstorm complexes improved as 
the day progressed. Figure 1 illustrates such a 
correspondence. 


1. Increasing the depth of the 
planets; y boundary layer. 

2. Indues 'g larger surface fluxes 
of momentum, heat, and moisture. 

3. Changing the vertical shear of 
the horizontal wind in lower- 
levels of the atmosphere. 

4 . Developing intense, horizontal 
convergence regions of heat, 
moisture, momentum, and cloud 
material. 

2. INTERPRET VTION OF FLORIDA 

MODELING AMD OBSERVATIONAL 
EXPERIMENT 

Still, the primal y ro] c of 
e ;Lra-stortn scale mesoicale systems is 
often viewed simoly as an initiator of 
cum'ilonimous convection. Once such 
systems form, they so perturb the plane- 
tary boundary layer by raean3 of intense 
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Fig. 1, The model prediotec vertical motion field at 1.22 
km and the radar ecno map at equivalent times for 29 June, 
1971. 




The interpretation oi this result 
is as follows. Individual cumulus congestus and 
cumulonimbus clouds can be initialed by small- 
scale moisture anomalies nol necessarily well- 
defined by a mcsoscale model. Once initiated, 
a cell may survive anu piopagute by virtue of 
penetrating downdrafts for periods of 45 minutes 
to 1 1/2 hours, A complex of cumulonimbus, 
systems, on the other hand, proceoses so much 
moisture that they cannot survive for an exLen- 
sive period (one to three hours) without the 
cnticiiment of moisture convergence by extra- 
storm scale mcsoscale systems. This inter- 
pretation has been corroborated by the analysis 
of several case studies, one of which has 
recently been reported by Cotton and Pielke 
(1976a, b). In this study, they found that the 
general complexes of cumulonimbi evolved and 
migrated across the Florid,, peninsula in response 
to the lower tropospheric winds similar to the 
predictions of the Pielke model. Individual 
cumulonimbus cells, on the other hand, were 
observed to propagate in response to deep trop- 
ospneric winds 50 to 75 km away from the main 
centroid of cumulonimbus complexes. In addition, 
the magnitude of the ratio of eddy transport of 
moisture at cloud base to tne 11 km mean moisture 
transport was found to be on the order of 2 to 
6%. At the same time, it was found Lhat the 
ratio of eddy kinetic energy Lo mean kinetic 
energy ranged from a high of 1700 in the vicinity 
of cumulus congestus and isolated cumulonimbi to 
a low of one below extensive cumulonimbus 
complexes. It is thus seen that thermal scale 
and cumulonimbus scale eddies make substantial 
contributions to the kinetic energy budget of 
the subcloud layer, while the cJoug base moisture 
fluxes are dominated by mesoscale eddies. 

3. OBSERVATIONAL AND THEORETICAL 

STUDIES OF THE COLORADO RIDGE- 

VALLEY CIRCULATION 

We are now in the process of deter- 
mining whether or now our conclusions from the 
study of the Florida sea breeze can be general- 
ized to other thermal) y-dnven mesoscale circu- 
lations such as the nage-valley circulation 
which forms over the Colorado Rockies. As m 
the case of the sea breeze circulation, the 
idea)ized two-dimensional solenoidal field set 
up by elevated heat sources is fairly well 
understood. Using a two-dimensional vorticity 
model, Dirks (1969) showed that a prevailing 
flow with vertical shear produced a pronounced 
tilt to the circulation cell above the slope so 
that part of the cell appeared to break off and 
form a much larger circulation cell out over the 
plajns. Figure 2 illustrates such a circulation 
cell. In addition, a simulation of the onset of 
the noctural regime demonstrated a strong 
enhancement of vertical motion over the plains 
as the downslope flow carried the upper level 
westerly moment-urn downward and outward into the 
plains. This is shown m Fig. 3. 

A.s we have found m the simulation 
of the Florida sea breeze, the actual three- 
dimensional flow is considerably mot e complex. 
This is particularly true of the mountain- 
induced mesoscale flow where the complex 
topography illustrated in F,lg. 4 can initiate 
intense dynamic interactions with the prevailing 


leol 


~s / ( 

■'-‘ s \ I J 


*\ , 

* ^ X/ 

7S \ / 

V "o v is -*i "" T 


I 

IS 



Fig. 2. Evolution of the deviation stieam 
function (solid lines in 10^ m2 scc _ l) and poten- 
tial Lemperature deviation (dashed lines in C deg) 
fields for Case D. Top of superaaiabatic region 
is shown by heavy dashed line (A) . 



Fig, 3. Evolution of the deviation stream 
function (solid lines in 10^ m2 sec - -*-) and poten- 
tial~Temperature deviation (dashed lines in C deg) 
fields for Case E. Time is in hours past start 
of diurnal sine wave, 

flow aloft as well as complicated patterns of 
elevated heating which can drive the upslope flow. 
To study the spatially and temporally varying 
pattern of surface temperature, cross sections 
of surface temperature were measured from air- 
craft with a Barnes PRT-5 IR radiometer. Figure 
5 illustrates an east to west cross section 
observed at 0900 LST on August 2, 1975, from a 
point nearly half way between Denver and Colorado 
Springs, Colorado, over the plains, westward to 
Leadville, Colorado. The data have been averaged 
over 11 km intervals and the variance with respect 
to such an average is shown m Fig. 5c. Not 
surprisingly, the regions of contrast between the 
plains and foothills, and mountains and plateau 
are the regions of maximum temperature variance. 
What is most surprising, however, is chat the 
South Park elevated plateau heats up considerably 
faster than any of the surrounding eastward 
facing slopes or the plains Lo the east. The 
magnitude of the difference m temperature 
between the plateau and the plains to the east 
reaches it.s maximum between 0900 and 1100 MST 
with values greater than 10°C. While the exact 
reason for the differences in response to solar 
heating of the relatively horizontal surfaces is 
noL fully understood, it appears Lo be largely a 
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F3S- 4. Smoothed topogiaphv (2 and 4Ax trends 
removed for a 1] km gild) used as input m Lnc 
Pielke mesoscale model, (contour interval 250 ra) 

function of differences in thermal inertia of 
the soil types. The importance of such a 
temperature difference is even more striking 
vhen one realizes that the plateua represents 
a large-area elevated heat source some 2600 to 
4000 feet above the surrounaing plains. Thus, 
it would appear to be a major driving force m 
the mountain upslope flow. 

Another characteristic of the 
pattern of elevated heating is that the plains 
to the east rapidly catch up to the elevated 
plateau exhibiting a surface temperature 
companble to the plateau. Shortly afLer local 
noon the plateau begins to cool, partly m 
response to extensive developing cloud cover 
and precipitation along the higher mountains to 
the west and partly due to the lower tnermal 
inertia of the plateau soil. During the same 
period, the plains to tr.e east continue to 
warm, finally reaching their ma—imum surface 
temperature around 1400 LST. Associated with 
the cooling of the upper level plateau is a 
reversal m flow from a generally easterly - 
ups] ope to a prevailing westerly current. The 
impact of this pattern of elevated heating and 
associated flow reversal on the genesis of 
Cumulonimbus systems over the western plains 
will be analyzed by combined numerical experiment 
wj tb the Pielke model and satellate and radar 
data analysis. 


The SMS-2 satellite photograph tit 
2345 CMT (1445 1ST) illustrated in Fig. 6 shows 
the extensive cloud cover which shrouds the 
higher elevations, as well as two preferred regions 
of convective activity over the western great 
plains. The regions of cumulonimbus activity are 
generally located over the Cheyenne ridge to the 
north and the Palmer ridge to the souLh. Also of 
possible importance is the fact that both ridges 
lie to Lhe east of major plateaus, namely South 
Park and North Park. 

SUMMARY AND CONCbUSIONS 

It is clear that thermally driven 
mesoscale systems swell as the Florida sea breeze 
and Colorado ridge/vallcy cjrculations play an 
important role in the genesis, organization and 
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Fig. 5. East-west cross section of IP, surrace 
temoeiature observed at 0900 LST on August 2, 
1975. Top represents smoothed elevation of 
topography. Middle represents 11 km average 
surface temperature. Bottom rcorosonts variance 
with respect to 11 km smoothed surface temper- 
ature. 
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Fig, 6, SMS-2 satellite photograph at 2145 GMT 
(1445 LST) with 1,85 kui resolut ion . 



intensity of cumulonimbus systems. lu the case 
of the Colorado ridge/valley circulation, how- 
ever, the importance of the circulation to the 
genesis and prooagation of severe storm systems 
over the western gre.it plains still remains to 
lie demonstrated- The emphasis in the oral 
presentation will be on the evidence indicating 
the role of such a circulation in the genesis 
and propagation of hail producing storms. 
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THE CHARACTERISTICS OF EVOLVING MESOSCALE SYSTEMS 
OVER MOUNTAINOUS TERRAIN AS REVEALED 3l RADAR AND FAM 
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Fort Collins, Coloraao 80523 


1.0 INTRODUCTION 

The 1977 Colorado State University (CSV) 
Sotztn Park Area Cumulus Experiment (SPACE) was 
a comorehensive summer t ire field orojecc 
designed to investigate the structure and 
evolution of cumulus clouds and convective 
mesoscale systems over mountainous terrain 
Most field measurements were taken m Soucn 
Park, Coloraao, a flat mountain park measuring 
about 100 km north-south and 50 km east-west, 
uirn a mean elevation of 2.3 km MSL Located 
about 150 kn southwest of Denver, South Pars is 
a natural genes ls area for cumulus clouus, with 

tne 4000 , MSL Mosouico Range d^recciv to the 

west and smaller ranges on tne southwest, nortn- 
west, and northeast edges. Within the park, 
intensive measurements were taken fron 10 July 
1977 to 13 August 1977, using powered aircraft, 
t a NCAR sailolane, rawmsonaes, boundary layer 
ro'iLers, acoustic sounders, micrometeorologi- 
cai towers, lidar, seatcn radar, tnple- 
DoDoler radar, and a square grid of surface 
meceoro logical stations, olacea at 10 km 
intervals In cooperation ■ ith the Bureau of 
Reclamation High Plains Experiment (HIPLEX) , 
am e tended area fron Soutn Park east to 
Oooclann, Kansas was also covered using search 
radar, satellite, ravinsonaes, and surEace 
nataorologiCal stations. Figure 1 shows the 
locations of search radars, raumsonde units, 
ird surface meteorological stations that were 
used m tne coraDinea SFACE/HIPLEN experiment, 
along with averaged 1000 fooc (304.3 n) 
terrain contours. 

The presence of the Rocky Mountains, 
ioothills and South Park immediately west of 
tne Colorado riign Plains creates an interesting 
sec of nesoscale -dynamic ana thermodynamic 
influences on the convective precipitation 
process. The initiation of small cumulus 
occurs earliest m the day over the Mosquito 
Range and ocner high peaks and ridges This 
is cue to the nigntcime drainage of cool, 
staolc air off the slocos, so that the 
morning solar heating causes the planjcary 
bounaar La>er (PBL) to deepen to cloud nose 
nut 1 - -ore ouickl, (ofeen n; 1000 IDT in 
Soul rk) thin over the pla.rs to the east 
I'hiv :eivc tru cool ,'rn-jy. mr 
Stror, aolar heating on east-facing slooes 
generates an upslope flow in tne lowest 
lavers, creating convergence of mass and 
noiscure abo.-e easc-facirg mountains and 


ridges. Dirks' (1969) two-dimensional vorticifcy 
model over a heated east sloee indicated that 
uooer return flow fron this upsloDe circula- 
tion would create a zone of subsidence 50 - 
100 km east of Che mountains, over an area 
where ooserved summertime rainfall is in fact 
less tnar. m aimer the nouncuns to the west 
or the plains to the east 

Thunnerscor-n activity is also observed 
to occur earliest near east-facing slopes. 

Henz (1974) u gea data from the Vgnonai 
Weather Service (NWS) 10 cm search raaar 
(WSR-5’’) at Liron, Colorado, to snow mat 
oreferraa regions of echo genesis are found 
near the east-facing slopes of tne front 
range Cunulonibus cells from these ' not 
spots" often propagate eastward, providing 
iruen of tne summertime precipitation vmcn 
falls on the nearby Plains Hourly rainfall 
data show an eastward shife of the time of 
maximum precipitacion, beginning lr the after- 
noon m tne Front Range, midnignt m western 
Kansas, and early cor-’ing in eastern Kansas 
This supports the hvpothcsis that early "hot 
soot" convection is responsible for Che 
formation of or participation in a dynamic 
meso-uave wmen propagates eastward 



Figure L Mao of the SPAC E/H[PtE\ extended 
coverage area, ntn ^onoutar 
averaged terrain contours every 
1000 feet (304 8 m) The large 
circles represent radar effective 
coverage areas 




July 19, 1977 was chosen as Lae first 
case tor inters rve analysis. On cnis day, 
strong cumulonimbus cells formed n South 
Park as well as otner parts of tne Coloraao 
Rock_es and m tne plains just east of the 
footnills, while ail convection was suppressed 
east of Lnon, Coloraoo Cells tenaea to ->ove 
norc-mr.-d witn t K e rean flow, r»3 rot 
initiating or participating m any one of 
eastward-moving -esoscaie dLsturoamce. A 
contrasting case stuoy day is 4 August 1977 , 
in whicn a strong squall line with high 
precimtacion forced m tne eastern part of 
Sootn Par< anc prooagateG eastward througn 
Kansas and Missouri We will concentrate on 
the analysis of 19 July 197 7 m the written 
presentation given here. Oar oral ora- 
santation „rll include a cornarstiva anal, sis 
of these two oays. 

2. RADAS AND SURFACE EQUIPMENT 


The National Center for Acmospneric 
Researcn Field’ Observing Facility i\C.AR/FOt) 
pro /idea tne Portable Automatic Mesonet (PAi ^ 
and tne' CP-3 5 5 cm Doppler raoar. PAM 
provided one minute averaged winds, wet and 
dry oulb tenperacures, oressure, ana rainfall 
at 20 locations, shown m Figure 2. '’ote 
chat three of the locations were at the cresc 
of the Mosquito Range at elevations of over 
3.8 km MSL. Data were teleneterea to a centra., 
site, where they were immediate*’/ availaole for 
d,splav and were also storea on tape (Cotton one 
George, 1978). 
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Figure 2 Map of Conch Park, snowing the 
PA '/CP-3 local coverage urei 
S- iwures 3-L5 hive ecactlv tne same 
scale md coverage as this nap 
Actual terrain contours are sho.-', 
with solid contours being even 
1000‘s of feet (304-3 ml, dasned 
contours 500’s of feec (.152 A o) 


Alchougn originally installed for triple 
Doppler use, t^e CP-3 radar aoubied as a 
searen radar oefore 1 August 1977 , wnen t'-“ 

13U FPS-L8, 10 cm searen radar, came into 
service. CP-3 woula periodically interrupt 
its Doaoler scans to perform tall volume scans. 
CP-3 data for 19 Jury 1977 are extremely 
detailed ard complete, and or /erv hign 
qualicy . 

3.0 LARGE SCALE ANALYSIS FOR THE MORNING 

OF 19 JULY 1977 

At I2C0 GMT (0600 .!DT) on 19 July 1977, 
a vast subtropical hign dominated the weather 
from tne Eastern U S west to Colorado. At 
500 mo, a branch of tms nigh was cantered i" 
eastern Kansas, witn resultant southerly flow- 
bringing not, dry desert air over eastern 
Colorado. In kyornmg and norenwesteen 
Colorado, a baroclimc zone provided souen- 
■ esterly winds wmch advectad deep Pacific 
moisture into the Rockies aurmg the day. 

.A wea< snortwave m this zone turned the 
Grana Junction, Colorado 500 -b winds to 260°. 
Good 500 no moisture extended from Grand 
Junction through Wyoming to Soucn Dakota On 
tne surface, a row ™as centered m Souch 
Dakota, with a surface trough e:. tending south- 
ward to cover eastern Colorado Figure 3 is 
a surface analysis for 0600 MDT, showing tne 
low oressure trough and associated light wires 
m eastern Colorado, with stronger winds at ere 
soutn ana east edges of the trongn 



Figure 3. Surface analysis for 0600 MDT 
(1200 CMT) , 19 July 1977, from 
NWS nourlv reports Pressure 
contour interval »s 2 mb. 

An east-west cross-section of 1200 CMT 
(0600 MDf) rawinsonucs is shown in figure A 
The cross-section extended for 680 km and 
includes rawinsondes from Grand Junction, 
Sout 1 ' Par-, sad - Linon, Colorado, and Gooulund, 
Kansas Potential cemperacure 9 (°K) and 
specific Tj-.-d it' q v (gm kg -1 ) are concoured, 
„rd «.Rd ure shown uc A0 nb intervals, dotn 
Limon ur.u Gooaland n id deep dry-ad lubat ,c 
layers, with low values of q v (5gm kg -1 -). 

This is tne result of several days of strong 
heating with little horizontal moisdru 
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auvection or moist convection. Tne Grand 
Junction sounding showed the oresence of deep 
Pacific moisture and the effects of deep 
cumuloniciDus activity which occured the 
Dteyious day and night (as revealed by 
satellite pnotographs) 'finds .n tne lowest 
SO mb were north or north-easterly at pocn 
Scum Pert, and Llnon, snowing tne effects of 
tne surface trouun over eastern Colorado 
Moiscure a; tne 350 tab level over South Para 
defined the location of a cirrus deck which 
covered the eastern slopes of the Rockies 
early in tne morning and subsequently fissi- 
on ted . 





barb *25 nt- Fi ^ , 

Figure 4 East-west cross-section of 0600 MOT 
(1200 MDT) 19 July 1977 soundings 
Lowest line shows aoproximace surface 
elevation. 


4 EVOLUTION OF THE PLANETARY BOUNDARY LAYER 

Five-minute average surface winds, 8, and 
q v observed witn PAM are shown m Figure 5. 

At 0600 MDT (before sunrise) , drainage winds 
pooled the coldest air into the center of South 
Park. At 0800 MDT and 1000 MDT, lignt winds 
were present in tne park, with stronger 
westerly winas on the mountain crest bringing 
a somewhat mgher q v into the park The 1000 
MDT South Park sounding (Figure 6) showed that 
the planetary boundary layer already extended 
up to 600 -ib, with 8 = 321°K and q v = 6.5 gm 
kg — ^ No subsidence inversion was present 
The cloud condensation level (CCL) las 540 mb, 
giving a wcc-oulb potential temperature 6 W of 
22 l°f This indicates a convective instability 
of aoout 1°C at 500 mo, with clouu tops of onl\ 
7-8 km If, however, tile mixed-layer q v was 

5 gm kg - ^ on this sounding (as was available 
on the 0600 MDT Grand Junction sounding), a 
8 W equal to 23 2°C at tne CCL of 580 mo, a 
500 nib instability of 2.5°C, and cloud tops of 
at least 12 5 '-an MSL would be obtained. A 
time section of the four South Park soundings 
(Figure 7) shows that the PEL hno deepened to 
500 no bi 1300 MDT, WLth 0 ^ 324 5 u rl arc a 

mi ed-l 'yer q v = 5 5 gm kg -L The decrease in 
moisture was due to the entramrent of drier 
air into the P!1L. Li'-on radar and satellite 
phoco.raphs show that eop convection nad 
begu efire noon to t a west ana nurenwest 
of South Park 
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Figure 5 PAM surface data for 0600 0800, and 

1000 MDT. 19 July 1977 DM mea-s 
data missing tCe 1 on page o 



Figure 6. South Park sounding at LOGO MDT, 
19 July 1977, skew-T, log-? 







Figure 7. Plot of the four South Park 

soundings vs time, showing cne 
probaoie evolution of the planetar-, 
boundary iaver over South Par< 
during Che day. Note the addition 
of hLgner values of q v after 1500 
MDT in the PBL Key on Figure 4. 

5. AFTERNOON RADAR A V D SURFACE ANALYSIS 

At 1230 -3T, the CP-3 radar began 
recording data Figures S-I5 are composite 
CP-3/PAM maps Echoes are contourea for a 
scan elevation of 5.5 = projected onto a 
horizontal plane, with only echoes of 30 dBz 
or greater intensity shown. PAM data and 
derived parameters are five minute averages 
aoout the time of tne radar scan. Locations 
of convergence lines or gust fronts are 
estimates iron 1 single-station time plots by 
determining the time of tuna shifts, temperature 
crops, etc., ana advectir.g the gust front aiorg 
witn tne speed and direction of the winds 
benmd tne gust front. 

Around 1200 MDT the first Dree ipita ting 
ceils formed m South Park along the west am 
northwest slopes. Figure 8 snows PAM and raoar 
presentations for 1239 :S)T ComDanng this 

V 



CoToired PAM/CP-3 racar pLots for 
123 n JDT, 19 July L977. Key to 
sy-cois on page 6. 


figure with Figure 7, we tan see chat Cne 
easterl_ winds Lnaicatea at low levels on Che 
1246 MDT Soucn Parti rawinsonae were confined co 
a small area of the park that was sheltered 
from the orevailirg westerly and southvasteriv 
winas av the mountains The oalloon was 
ooser- ed to ascend "ear a cumulonimbus cloua, 
probaolv the nook snaped echo north of the base 
shown n Figure 8. Several convective cells 
were active in the north end of the park. -with 
one strong precipitation area near CP-3 (at 
this close range, the scan is onlv a few 
hundred meters aoove ground level). All of 
these cells were moving towara the north at 
about 5 ns - ^, aooarenciy in response to upper 
winds 


Small racar ecnoes continued to form 
northwest of the base and west of CP-3, in 
the climatologically most favorea echo 
locations (Huggins, 1975) One cell was 
ODserved on radar to intensify greatly as it 
approached Hoosier Pass from the south (see 
Figure 2). By 1333 MDT ("Figure 9), a line 
of strong thunderstorms nad formed in the Blue 
River valley north of Hoosier Pass. Range 
Feignt Indicator (RHI) scans, which were 
synthesized from a series of constant elevation, 
full volume scans (Knupp, et al, 1978) show 
that these cells had tops of over 12 5 km, 
wnereas the smaller cells within South Park had , 
tops of 7-8 krt This indicates thac the 
increased oounaary layer moisture (8 gn kg - -} 
needec to supDort deep thurcerstorns, as 
■"ndicaced by ravmsonde data, was available 
just to the northwest of tne park PAM data 
in Figure 9 indicate that this moisture was 
being advected across the mountains The 
remainder of t^e nark, ecciuairg the north ana 
west edges, continued unaar strong heating, 
witn surface q v decreasing to 6-7 gm kg - - as 
drier air continued to be entrained front above. 
Litnon radar data for 1330 MDT snowed tnat the 
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Figure 9. Same as Fig 8, f 
July 1977. 


Figure 8 
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Hoosier Pass storms were pace of a larger 
convective oana mien ewtenaed nortneastward 
for 150 km. 

Ac aoout 1430 1ST, surface winds m tne 
southeastern oarc of me park pagan co snift 
coward the ease t -;s f_av gradually aftecrcc 
Che encire eastern enga of the park. The 
resultant convergence pattern generated 
nu-ero-s snail corvsccive ceils from 143G- 
1600 MDT. Surface te-oeracures renamed nigh, 
inuicatmg strong neating ana increasing 
instability in the center of Soutn Park. 

Around 1600 MDT the easterly flow mvadea cne 
center of the park, and together with the 
stronger westerlies m me western portion of 
tne park, forreo a convergence line wmch 
generated a line of thuna era corns. At this 
time, two incense cells formed at the noren 
end of the park near wenosna Pass. 

In Figures 10-15, the convective cells 
are numbered sequentially relative to their 
pcooaole initiation times. Ac 1623 MDT 
(Figure 10), four cells (C3-C6) were visible 
over the convergence line, wirn the second of 
the strong cells tC2) farmer „orch. Mote that 
C2 moved well co the west (left) of the upper 
winds, while the smaller cell C4 moved from 
160°, 3Doro--i-'3tely along tne line of con- 
vergence PAM ram gauge uata snow that 
station 23 received 1-2 mm or precipitation 
per rr.nure for at least 15 minutes wmle cell 
C2 was overnead. A visual ooservation of a 
hail snaft from this cell was reported at 
about 1600 "SDT A mesoscale cold front nad 
formed soutn of cell 02, pernaps consisting of 
outflow frov C2 and ceils farther north, Wina 
velocities in tne southwest seccor had in- 
creased co 7-10 ms - -, and considerable 
cvclonic wind shear existed near tne frontal 
boundarv (cell C3) . 

By 1702 MDT (Figure 11) , the meso-cold 
front had roved abouc 20 km south on tne west 
siae of tne convergence line (traveling at 
3 ms }. In succession, cells C4, C5, C6, 
and C7 haa intensified to over 40 dBz 
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Figure-. L0-I5. Combined PAM/Cp-3 ndir plots 
for 1623 '!DT througn 1939 MDT 
at 20 -mute internals key 
co s'. —do Is above. Arrows are 
used to sjggesc possible non- 
quancitative streamlines 


refleccivitv as tne frontal winds approached 
Cells C4 and C5 haG oarcmllv aeua. ea o> 1702 
MDT These cells -ovea from aoout 170° at 
5 ms - ^, with new cells naviog continuous! 1 , 
formed on the south end of the ime. Pre- 
cipitation at rates of up to 100 — n nr -as 
measures at the west ease ar uo kote the 
weak, radar recurn directly over tne Continental 
Divide, wit 1 ' remtens^tied ecnoes over the 
31aa '"alley norenwese of tne marie. 

At 1742 MDT (Figure 12), surface winds 
were more or less northerly througnout the P\M 
network. The latest cell (C9) was also the 
most intense (over 4S d3z). 'vest of the 
convergence line, wines were stronger, and 
values of 8 e (equivalent potential temperature) 
higher, chan that under ana east of the line 
of cells. 

By 1822 12)T (Figure 13) the character of 
the lme had begun to change Cell C8 had 
split into two pares, witn the more intense 
part (C8) moving at 150° while the smaller 
cell (C8a) movea almost due roren. A gust 
front, possibly causea oy outflow from CS, C9, 
and CIO, was -ovine into cne oases of new 
cells Cli and C12. Cell Cll grew very quiekl- 
mto the aonmant storm of tne oay. By 1856 
MDT (Figure 14) a noo^-snanec ecno protruded 
to che southwest of Cl'l A surface scation 
under this hoot, e 'oenenced precioitacio"’ rccas 
of aoout 200 mn nr -1 , prooaoii including nail 
Cell Cll movea at 4 ms - ^ from 125°, ~ell to 
cne left of the mean winds Moist outflow frc-i 
other cells concmuea to feed this cell from 
cne nortn Cell C3 moved fron 150°, while 
other ceils in the lme movea even farther to ch 
ngne. Cell Cll continued co move slowl; and 
— intensified to at least 56 d3z at 1922 MDT, 
wmle the otner cells m the lire decreasea m 
intensity 3y 1939 MDT (Figure 15), cell Cll 
had oegun to decrease in intensity wmle 
increasing its speed to 10 ms - All the other 
cells m the line had dissipacea, while a new 
line of cells, oriented east-west but still 
moving toward che north, had oegun in tne 
northeast part of Che nark. Winds in South 
Park became ligneer and more vanaDle after 
the passage of Cll, witn cold, moist air 
occupying the area. Althougn no furcher radar 
data are availaole, PAM rainfall data indicates 
that Cll may have moved over tne crest of tne 
mountains at about 2015 MDT 

Total precipitation oaca for tne penoa 
1600-2400 MDT are shown m Figure 16. Vote 
the ma-cima at the e<cre-e nortn (.from C2) , 
ce-cer (from CS) , a-d excre-e south (fron Cll) 
RamrulL at station ‘t 28 raa\ have been higher 
than snown because it occurred at 1600 "DT, 
when there were voids m cne PAM aaca. 

6. LARGF SCALC AVALISIS FOR U00 MDT, 19, 

JL'U 1977 

A very conple' picture of lower crooospher. 
f low patterns emerged -itn che analysis ot 1800 
MDT (0000 Gill, 20 juii 1977) svnnpcic and lame 
mesoscale data. Ac 500 -o, southerly flow hac. 
resumed enrol gh rust of Colorado, but the Soucm 
Park sounding tor 1710 MDT still had a 235° 
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wind a c 560 rb (Figure 7), indicating cor- 
tin Led advection of Pacific moisture across tne 
mountains. Surface aata (Figure J.7) shoved a 
strong lev pressure center in the plains 200 -on 
soucncasc of South Park, parr of rme pars is rend 
eastern Colorado surface trough This low* 
pressure cencer naa deepened curing tne day. 
probably in response to ootn surface neatirg 
ana a weak short wave vnxen may nave passed 
througn cne region. 

Satellite nnotograpns na<en at 1330 'IDT 
show the entire Rocky Mountain region covered 
with clouds Tne deepest storrs were located 
in nree north-soutn bands, eacn extending t'rc-a 
central Coloraao into souenern k'yonirg. The 
easternmost of these bards extended aorc i-ror ch- 
eastward fron Soutr Part to a point norm or 
Chevenne, Wyonng This convective sand 
probaolv marked the easternmost advance of 
moist Pacific air in the planetary boundary 
layer Surface data iron 1500 HDT - 2200 .T5T 
show chat the first occurrence of strong 
thunderstorms east of the Rockies was at 
Chevenne at 1500 'IDT. At this time, Cneyenne 
winds saifteu to northerly, and tne tempers I are 
drooped about 10°C. This phenomenon then 
occurred sequentially at each station farther 
south along the front range. At 1300 MDT, 
the Boundary of tnis wind smft had passeu 
soutn of Denver, T 'nile tne nortrerl* wire! 
stuccs oDServea m South Park naa also passed 
througn In Figure 17, a cold front symbol 
is used to show the possible location of tn*s 
air n* n ss boundary 

Tne Denver sounding for 1S00 MDT showed 
a layer >i Xignt nortn-norcheasterly wirgs 
e<tanding from the surface (at 640 mo) to count 
700 mb This laver nad a potential temperature 
of about 32I°k and an average q v of aoout 7 5 gm 
kg -1 In Soi ch Park at 1710 'IDT (see ? igur » 7) 
there were norc u erlv and rorcneasceri r *irca 
from the surraco (719 mb) to about 690 mb, 
easterlies above tnat to o00 id (perhaps an 
eascu'-lv :low regime driven bv tne soueneast 
Colorado surface low), soucnwescerlies from 
600-520 mo, and soutnerl/ winds aloft This 
sound mg was laun*_ned only a few -‘mutes after 
the northerly surface winds covered the Soutn 
Park rawinsonde site 


Figure 17. Surface analysis for 1800 MDT, 

19 July 1977 (0000 CMT 23 Juiv 
1977) from WS Hourly reports. 
Pressure contour mter/al is 2 
mb 

Earlier in the day, at 1333 >Df, Soutn 
Park radar and PAM data (Figure 9) nac shown 
that tne mountain barrier extending to aoout 
Che 650 ro le/el (norc«nest of the Par-.) nac 
been sufficient to block /er/ moist Pacific a_r 
from advecting into the park. This moistvte 
hac caused deep chunaerscorm activity just 
outside the earn, vet the onset of aeep scors 
in South Parti was delayed for at least three 
more hours This when conomea with evidence 
from the sounding wine fields at Denver and 
South Park, indicates chat advection of Pacific 
moisture across the Rocky Mountains took place 
at very low levels, Droocoiy below 700 id. In 
souenern tT yommg, a gap m tne Rockv Mountain 
barrier allowed such aavection co proceed "'ote 
rapidly into tne Cheyenne area. 

After 1500 MDT, deep tmunaers corns, wnich 
had beer occurring previously in tne mountains 
west of Denver, invaded the plains alon* the 
front range. Persistent northerly flow in tne 
lowest levels was observed at to is time m a 
bana about 100 km v’ide This was an area of 
east- west oaroclinicity in tne lowest levels, 
while at 500 mb, little heignt gradient was 
apparent Low level air to tne west of this 
zone vd s cooled by a como matron of tharma.^. 
advection, evaaor^tmg tnunderstorn downdrar^s, 
and lack of surtace heating due co cloudiness 
To c e east, extreme heating continued in tne 
doud-frea zone. The resultant surface 
pressure gradient, shown in Figure 17, 
prooabLy provided the forcing for Che nortaariA 
v mds 

i. OVCiL'SXOS 

Jnlt L9, 1977 was selected as our first 
t ** siii'.v d/i*' for storm-scale and me c os<_ale 
a 1 v*ls beciuse of its ipparent relative 
simplicity. Extensive aial*sis ot this d >v 
r vc lied a very complex fie turn of tHund* r-* 
scorn scale, large-mesoscale and synopcic- 
seale interactions. 



Eariv convecfnn "is initiated over fe 
’!osqu-to Range in Che western and norcnwestern 
portion of Souch Park, m the roist, Pacific, 
westerly-flow airmass As cne cav progressed, 
easterly flow, driven m part o> tee surface 
'iow southeast of Souch Park and in oart by cne 
o-urnaiiy driven nage-valley circulation, 
invaded the Park and sec up a north-soucn 
oriented convergence line. Soutn Par*. was 
apparently the location of the convergence 
line Decause it remained cloud free for a 
mz-or portion of the day and thus became a 
significant elevated heat source. 

Late m the afternoon a cooler airmass 
began invading the northern end of tne Park 
At first it appeared to be outflow from a 
single cell As time proceeded, however, the 
norcherlv flow became quice persistent and 
coherent across the p ark. By 1742 MDT 
tne cool, norther! ; -flow eirmass nad mvaaed 
the entire PAM network. Because of the 
temporal coherencv and spatial consistency of 
tne northerlv-flow airmass, it resample- a 
synODtic-scale frontal oassage. The 
occurrence of such a frortai passage was 
conceivaole since a weak surface front • as 
analyzed on the morning surface maps along Che 
southern Wyoming border. By 2000 MDT, nowever, 
Che nortnerly flow over the PAM network had 
weakened in intensity and become less consistent, 
suggesting a diurnal influence. 

Tne difficulty m isolating mesoscale 
versus svnoptic scale influences is that tne 
synoptic-scale baroclinicicy was so weak that 
it cannot be opjectively uefined At tne same 
t^ne, the Colorado Rocky Mountain barrier is so 
large (in its vertical and horizontal extent) 
that its associated dxurnally-dnven raesoscale 
disturbances develop a distinct synoptic- 
scale cnaracter Thus considerable further 
analysis and perhaps numerical experimentation 
is required to isolate the scales naving the 
major controlling influence on tms day. 

It ts interesting to note tnac the north- 
soucn orienceu line of convective cells wnich 
was preceded by a north-south oriented surface 
convergence line, retained its general north- 
south orientattve during and subsequent to the 
passage of a well-defined, larger-scale, 
generally east-west oriented convergence line 
The reraarkiole strength of the earlier observed 
mesosrale convergence line is attested to by 
the fact that this feature was sustained 
during the oassage of a vigorous larger-scale 
convergence line. 
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develooea tne ?\_1 system and CP-3 radar. Tne 
excellence of tne PAM and CP-3 data is large!, 
due to the hard core performed by FOF field 
personnel 
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A RAD AX CASE STUDY ANALYSIS OF A HEAVILY PRECIPITATING 
QUASI-STAXIONARY CONVECTIVE STORM SYSTEM 


K.R. Knupp, iC S. Danielson, W.R. Cotton 
Department of Atmospheric Science, Colorado Stete University 

Fort Collins, CO 80523 


1.0 INTRODUCTION 

i 

i An integral component of the 1977 South 
Park Area Cumulus Experiment (SPACE) , was tne 
NCAR Field Observing Facility (FOR) CP-3 C-band 
radar. This Doppler radar system was used to 
obtain reflectivity and racial velocity data 
from s-mmercime convecti"e storm systems occurring 
daily within Souch Park and adjacent areas. The 
NCAR/CP-3 radar was nart of a trxrle-Doppler 
racar array which included the two NCAA Wave 
Propagation Laboratory (WPL) X-bard radars. 

Pig. 1 illustrates che locations of these radars 
witnin Souch Park. 

On 19 July 1977 an intense and long-lived 
convective storm system was initiated and main- 
tained by phenomenon organized on the mesoscale. 
The general synopcxc-reateras accompanying tne 
mesoscale system were low vertical wind shear 
and moisture advection carriea by relatively 
weak southerly flow. George and Cotton (1978), 
m a companion paper. Have investigated the 
general features and evolution of this mesoscale 
system. The morphological characteristics of 
the July 19 radar echos, wmen is the subject of 
this paoer, will form a basis in selecting cells 
for subsequent triple-Doppler radar analyses. 

Several distinct and intense convective 
precipitating systems occurred on 19 July 1977 
These convective entities had unusually long 
lifetimes, high reflectivities, and high echo 
cops when compared with climatological radar 
data from Soutn Park reported on b> Huggins 
(1975) Two specific intense cells, which differed 
greatly m echo structures ware probable hail 
producers Ono of the convective systems was 
composed of multiple nigh reflectivit) cores, 
etibte. 1 non-steady-stace features, and preceded 
Che fc , non of convection organized on the 
mesoscale. The other convective system formed 
afeer initial mesoscale organization haa taken 
place. Contrastingly, this system displayed 
quasi-steady-state echo characteristics. 


2.0 GENERAL ECHO CHARACTERISTICS 

Initial convective activity on July 19 
occurred shortly before noon in the northwestern 
quadrant of South Park along the eastern slopes 
of che Mosquito Range By 1530 (all tires are 
MDT) the major area of convection had shifted to 
the northeastern portions of South Park, where 
rapidly growing cumulus merged and formed an 
intense convective system composed of several 
hign reflectivity regions. 



Figure 1- The location of the CP- 3 C-oz'-d r~icr 
and the HOin 1-ba.na radars within Pc it a Pc. . 

The S-bana C3S-4 radar a team operatvO'zzl a; ,'r 
3 August 1977 and was used for iurvsilla'ioe co- 
poses. i denotes radar locations , • aerates ? 
station locations. Labeled elevation contours 
are thousands of feet hSL 



t — 



Figure 2. Constant 2. S'" elevation radar seems from 1537 'IDT (2a), 1803 !’DT (2b), and 1907 !1DT (2a). 

Contours ere Cretan at 5-dB intervals beginning at 20 c3Z. The dotted areas represent reflec- 
tivity factors greater than or equal to 35 d3Z. 


As this system intensified, a broken north- 
south line of isolated convective ecnos over 
eastern South park extended southward from the 
intense convective system in the northeast. The 
develoomg line at 1637 is shown m rig. 2a. 

This ime remained m a quasi-stationary position 
over eastern South Park for the next two hours. 

By 1304 I1DT the echo line had attained a length 
of M.00 km and a mean width of 10 —15 km. Fig. 

2b clearly shows that by 1804 MDT the most 
incense convective cells were located at tne 
soutnern end of the line as opposed to the 
nortnern end at earlier times. A typical defin- 
able convective cell in this line would form on 
the Southern end of the quasi-stationary echo 
line and propagate northward through the line 

with a mean speed of 'Mj ms . Two notable 
exceptions to this sceneno formed and intensified 
within the line’s interior. 

Between 1830 and 1900 I33T the convective 
echo line underwent another transition as a 
larga organised, intense convective system on 
the line’s soutnern end began to nominate. As 
tms system intensified, cells to the north 
weakened and the organization of the line, vmch 
had been a persxstent quasi-stationary feature 
for a,2 hrs. rapidly dissipated. The 1906 5.5° 
elevation scan is shown in Fig 2c. Noteworthy 
features are the organized appearance of the 
convective system at the line’s southern end and 
the lack of intense echos (only one appears) to 
the north. A closer examination of Figs. 2a - 
2c reveals a clockwise turning of tne line’s 
main avis This effect may be oua to local 
changes in convergence of the associated neso- 
scale system. The northwestward proDngation of 
the dominanc convective region m Fig. 2 c may be 
a factor .n the line’s axis cnange during the 
latte ' t.-e periods. 



Time (MOTl 

Figure 3. The total echo areas , enclosed by 35, 
40, 45, and SO d3Z contours, as a f.aiction of 
time The ecno areas include only the segment of 
the eel jo lice south of the CP- 3 radar. 

A summary of the major convective cells on 
19 Jul> is Dresented m Table 1. These data, 
when comoared with Huggins (1975) data, indicate 
that cell cnaraccenscics, including echo life- 
time, maximum echo top, and ecno intensity are 
well above normal for South Park thunderstorms. 

In his analysis of 1973 South Park radar data, 
Huggins found chat convective echos haa mean 
echo too heights of 9.0 km MSL, with extremes 
of 13 0 kn MSL. As shown in Table 1, echo cops 
on July 19 often exceeded 12 and 13 kra. According 
to Huggms, -lean ecno lifetimes of South Park 
cells were 30-40 minutes, with cells infrequently 
lasting 30-90 minutes. Table 1 inaicaces " that 
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TABLE 1. Summary of the major TRW subsystems occuring after 1400 MDT, 19 July 1977. 


Subsystem ■ 

Tima of Ecno i 
Aooaarance (MOT) J 

Time of Echo i 
Absence (MOT) 1 

t 

Li retire (Min) > 

Maximum Tod of 1 

20-d5Z Ecno 'km MSL) 1 

Maxi mum 

Intensitv (DBZ) 

1 i 

1 

7 

7 

7 1 

\ 

? 

50- 

2 

1500-1510 

| 

? j 

>120 | 

>12.0 | 

So 

3 

1500-1510 

7 

\ 

5 1 

7 

50 

4 

1530 

1740 

140 

13 5 

51 

5 

1600 

1900 

180 

12 0 

51 

6 

1610 

| 1830 

140 j 

n 5 

«7 

7 

1630 

| 1900 

150 j 

13.3 

50+ 

7 a 

| 1730 

! 1850 

: 

[ 80 

12 2 

! 50 

8 1 

| 1645 

1945 

180 

14 1 

51 

9 

! 1710 

i 

1304 1 

- 

13.3 

51 

9a 

j 1720 

1740 2 


i 12 8 

47 

10 

1740 

1910 

90 

13.5 

51 

11 ! 

1740 

well after 
2010 

; >150 

13 7 

t 

1 54 
i 

12 

I 1750 

1 

1900 

80 

12.5 

45 

Average 

i 

1 

1 

130 

] 12.8 

51 


1 Merged with cell 10 

2 

Merged i/itn cell 9 


July 19 echo lifetimes often exceeded 140 minutes. 
Altnougn ettensrve climatological studies of 
echo reflectivities arc limited for Soutn Park 
thunderstorms, Che maximum intensities appearing 
in Table 1 appear to be well above normal and 
aporoach the echo reflectivities measured within 
moderately intense bE Colorado thunderstorms. 
(Foote, « al , 1976.) 

Fig 3 presents the time evolution of the 
total area enclosed by 35, 40, 45 and 50 dBZ 
echoes after the echo line had attained an 
organized st 2 te. The values plotted in Fig. 3 
were oocained from 2.5° constant elevation scans 
fron Cae portion of the echo line south of the 
CP-3 radar. One of the noteworthy features 
appearing in Fig 3 is the relatively constant 
area of echo intensities for the 1700-1830 time 
period This constancy reflects the quasi- 
steady-stace properties of the echo lire as a 
whole. However, for the sane time period, 
individual convective cells within the echo line 
exhibited widely varying echo shapes and in- 
tensities The nigh intensity echo areas reached 
a mati~--a near 1930 1ST resuitirg primarily fron 
the presence of the intense quasi-steady-state 
convective subsystem described earlier. 

The following sections discuss in further 
detail the individual convective subsystems 
wnich exmbited differing echo characteristics 


(The tern convective subsystem is defined here 
as a convective entity composed of one or more 
cells of fcign reflectivity regions Most of the 
subsystems listed in Table 1 appeared to be 
single-celled at a given instant in time.) The 
three suosystems to be described were the only 
ones on the 19th vnose trajectories were signifi- 
cantly to the left of the mean environmental 
cloud-level winds ('VL70“). These subsystems 
were also quire intense when compared to "average" 
July 19 convective subsystems. 

The first subsystem developed prior to the 
formation of tne aesoscale echo line and dis- 
played multi-cellular and non-steady-state echo 
feacures. The second subsystem studied was the 
most dominant one occurring on 19 July. It 
attained quasi-s tcady-state echo characteristics 
and aopearea to communicate with and modify 
certain mesoscaie features. The third subsystem 
was long-lived, and appeared to exhibit non- 
steady-state echo characteristics as the result 
of its interactions with adjacent cells and 
interactions with the mesoscaie 



3.0 ANALYSIS OF SPECIFIC ECHOES 


3 1 Production or a Multi-Cellular Storm Dy a 

Merger Process 

Intense reflectivity regions formed in a 
variecy of ways on July 19. One moce of develop- 
ment vis cnrougn successive mergers of norp-.olog- 
ica±i> mdepencent reflectivity regions. The 
development of a convective subsystem (C2) shows 
en excellent example of Intensification via cell 
merger. The individual cells whicn eventually 
merged to become C2 first became prominent at 
1510 before the echo line described in the 
previous section became organized. PPI data 
suggest that this subsystem resulted from a two- 
steo merger Drocess Approximately 3—10 separate 
reflectivity regions existed initially and 
underwent mergers in groups of 2 to 3, which 
resulted in three larger and more intense cell 
groups. RHl's of the two most intense cells in 
two of the cell groups -.t 1532 incicate tnat tne 
maximum reflectivity core of each lies 3-4 Km 
above tne surface. The elevations of these re- 
flectivity maxima at 6-7 Km MSL are slightly 
below ouc consistent with the median values 
observed for che South Park region by Huggins 
(1975) and for NE Coloraao by Dye (1976) . 

The 1542 elevation scan shown in Fig. 4 
portrays the cell complex C2 shortly after 
merger of the three cell groups. These cell 
groups are now apparent as three maximum reflec- 
tivity regions labelled a, b , and c. Also shown 
in Fig. 4 is tne location of Cl, a less intense 
convective subsystem which was prominent before 
C2. It was noted that anotner cell group to tne 
south weakened and dissipated as the merger and 
subsequent intensification of C2 occurred. By 
1552 C2 contained one large reflectivity region 
with a reflectivity factor maximum of 58 dBZ at 
3.75 Km MSL. It is likely that this reflectivity 
core contained small hail (or graupel) at this 
time The 5.5° elevation scans at 1553 and 1602 
indicate that the 40 dBZ echos of C2 hao completely 
merged to form one large echo. Other noteworthy 
features were the tight echo gradients appearing 
along che southeastern edge of C2 during its 
most intense stages (1652-1602). This is • 
contrary to the usual pattern of nign reflectivity 
gradients appearing along tne leading edge or 
inflow area (the northwest quadrant in this 
case) of intense convective storms. (See Foote 
and Fankhauser, 1973, for an example.) 

During C2's mature stages, the general echo 
structure was characterized by a 45-55 dBZ 
reflectivity maximum situated near the surface 
along the system's southeastern edge, and a 
seconder, suspended reflectivity maximum located 
aloft (6 to 8 kn MSL) along the western or 
norchwestern edge At times a third local 
reflectivity maximum appeared at an intermediate 
altituue Detween the ocher two maxima (see Fig. 

5) Tne average 10-mm. time interval between 
successive PPI scans was too great to track che 
evoiut.nn of these maxima The location of tha 
major nf low/uodraft region along the northwestern 
quadr c of C2 may explain, by propagation 
efrects, why its trajectory was significantly to 
tne left ('150°) of tne mean cloua-ievel environ- 
mental winds during its incense stages. 



Figure 4 The 1542 'OT S. 5° el3Vc.zz.cn scan of 
convectzve subsysteri » oncer 2 (C2) shortly after 
the merger of three pre^tc^slu separate ceil 
groups The celt groaps, rcu euzdsnt as three 
reflectivity Knit, ere labeled a, o, and a The 
reflecr.tvzt i factor consolers are at 25 , 35, eO, 
end 50 dBZ, noth dotted regions greater than 40 
dBZ. Subsysteri 1 is located ,B of C2. 



Figure S. Core Pi" of C2 recorstri cted fret the 
16~2S-50 to the ISZC Cl "IT sector scan Contours 
are drami for 20, 30, 35, 40, 45, arc 50 dBZ. Tne 
dotted areas represent reflectivity factors greater 
than 40 dBZ. The slcoted lire depicts the mccaitum 
elevation (15.5°) scanned ig C?-5. 

3.2 Single-Cell Growth in an Environment Modified 

by Previous Clones 

The subsystem described m this section 
appeared in a region previously occupied by 
several relatively passive echos at the southern 
end of the well— developed echo line desenoed 
earlier This passive ec.no group had moderate 
radar reflectivity returns for approximately 20- 
30 minutes before in; significant intensification 
of tne ecno occurred Tne RHI profile at 1814 
PET indicates tnat this cell (Cli) was multl- 
turreted, with an overhanging echo along che 
northwestern and sout-jra quadrants. A merger 
phenomenon is not evident m this case. Using 
surface mesonet data, George and Cotton (1978) 
.suggest. that intensification of Cll coincided 
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with che arrxval of outflow air from cells 
located 10-20 ka due rorc- E animation of che 
morphology of me ceil (CIO, Fig 2o) from wmcn 
che outflow is nyoochesized, indicates chat ac 
I 8 O 4 , this intense subsystem had an extensive 
overhanging -43-50 dBZ echo along its eastern and 
southern portions. By 1815 mien of the overhang- 
ing ecro had collapsed and maximum intensities 
appeared near the surface along the subsystem's 
southern quadrant. The outflow resulting from 
this collapse, as inferred iron the above descrip- 
tive echo evolution, coincides with the gust 
front detected by the surface nesonet and analyzed 
by George and Cotton. 

The 1837 and 1848 MDT core KHI’s show 
several interesting features of Oil's echo 
morphology in tne vertical Unlike the 1814 RHI 
of this cell, these later data indicate a single, 
intense, overhanging echo core- The secondary 
individual far-sioe echo appearing on tne 1848 
RHI in Fig. 7a is associated with che hook-like 
echo on the soutnern quadrant of Cll. This hook, 
echo, uhlcn xs shown in che 1348 5 5“ elevation 
scan of Fig. 6 a, as pears to oe a flanking line 

of ro-arirg cux- an-amariag from Cii's ec.no 

core region 



Figure 6. Constant 5. 5 ° elevation scans of con- 
vective suts*stem number 11 (Cll) at 1343 MDT (6a) 
and at 1916 IDT (5b). Contours are dram at 5 d3 
intervals beginning at Z0 aBZ The dotted areas 
represent reflectivity factors greater than or 
eaval to 40 d3Z. The disced Xic.es are the hori- 
zontal locations of the RHI's appearing m Fla 7. 

Examination of subsequent PPI's and RHI’s 
reveals chat this subsystem raa the characteris- 
tics of a large ouasi-sceady-stata, intense, 
single-cell storm Steady-state features of 
this cell include a fairly persistent FPI hook- 
ecio pattern (Fig. 6 ), a single, intense reflec- 
ci\ic> core often associated with a weak-echo 
region (WER) to the north or northwest (Fig. 7), 
and a nearly consume maximum ecno top (Fig. 8 ) 
The 1928 RHI and 5.5“ elevation scan revealed 
the oniy exception to the otherwise constant 
single reflectivity core These data show a 
double reflectivic core. The secondary reflec- 
tivity core (50 dBZ) appearing at this time was 
locacw- in an echo overhang in tne cell's north- 
west • .i^ranc. Ic should be noted chat tne 
overha gitig echo structures appearing in Fig. 7 
are consistent with other spatial and temooral 
core structures of Cll from 1840 to 1950 



Fiaure 7. pccorstractea. PRI cross sections of Cll 
from the 184c JD? (7a) a"d 1915 - ID? (7b) sector 
sccws, Che vefZeezsiciszy factor qosi~om*s eve 
for 20, 30, 55, 40, 45, and 50 dSZ The corre- 
sponding FPZ t s atoear m Fig. S. Ti should os 
rioted that a 4-mn. time period uas taker to ob- 
tain this reflectivity data. 



Tims <>!DT, 

Figure 8 Tne tire evolution of Cll 's r -aziim 
ecro intensity in dBZ (top) and maximum 20 aBZ 
echo top (batten). 

The echo-inferred inflow ireas of -Cell 11 
were initially confined Co the southern sectors. 
The inflow region transferred Co the northern 
and northwestern portions as the cell intensified 
At 1837 FfflT the major inflow appeared to oe 
divided almost equally between che southern and 
■northwestern portions. - - — — 


t\ox Echo Top (Km MSL) 





At 1837 a ridge of 13-km 20-aBZ echo tops 
extended from the south-central portion of Cell 
11 to the north-central oorcion. By 1906 three 
local 13-xn maxi-a appeared as toe major inflow 
appeared to shift to the northwest quadrant. 

The 1926 echo-top contours reverted to a single 
area of 13-kn tops. These relatively minor 
fluctuations (both spatially ana temporally) in 
Cell 11' s echo cop maxima seem to indicate a 
slignt non-uniformity in cell’s updraft and 
inflow patterns. Despite the relatively minor 
spatial non-umforaities m echo cores, echo 
structures, and echo tops, the persistent hoos. 
echo and RHI core ecno features indicate that 
Cell 11 was in a steadv-state condition for the 
period 1840 to 1950 >21. 

3.3 A Non-Steady-State, Long-Lived, Interactive 
i Subsystem (CS) 

, Most convective cells interact to differing 
degrees with surrounding convective systems. 

This example provides some clues on tne important 
and complex problems of cell-cell interactions. 
Many of the physical processes wnich conscience 
this interaction nave not oeen well defined. 

The morpnological description in this section 
presents a hignly interactive subsystem as 
inferred from reflectivity data. 

Throughout nuen of its lifetime convective 
subsystem number 8 (C8) displayed both vertically 
and horizontally non-steady-state echo features. 
The initial echo of C8 first aopeared at the 
soutnern end of tne ecno line and coved slowly 
northward through tne line as furtner ecno 
development occurred to the soutn. As the 
subsystem orogressed through the ecno line, 
moderate cnanges xsi reflectivity maxima, in 
reflectivity structure, ana m ecno areas were 
noted 

Sequential KHI cross sections from 1730 - 
1837 MDT exemplify the variability m acho 
structure and inferred flow patterns of this 
cell. Between. 1730 and 1745, the major inflow 
Co Ceil 8, as inferred from an overhanging acho 
structure, shifted from the cell’s eastern 
boundary to its western boundary. C8, which had 
a width equal to that of the echo line, was 
located in the south-central portion of the echo 
line at this time. For the remainder of C3's 
lifeci-’e, major inflow appeared to stay along 
the western and northwestern ecno boundaries 
During this same time, the subsystem (C8) itself 
remained proximal to the western Doundary of the 
echo line. The transition of this inferred 
inflow region bears a striking resemolence to 
Che inflow transition exhibited by CU. The 
availability of low-level moisture was probably 
important m CS’s longevity C^3 hrs.), slow 

average speed of movement (3-4 ms , and 20° 
trajectory deviation to the left of the mean 
winds. 

No significant aaarun reflectivity or ecno 
struct re cnanges accompanied C8's inflow transi- 
tion '.aver, an intense strong cell (C9) did 
develop a v ir CS during the 1730-1745 transition 
period. Tne 1740 RHI vertical sections of 
C9 indicate the presence of extensive echo 
overhang along tms echo's eastern edge. Thus, 
it, appears that C8's inflow change was con-current 


with the appearence of an active cell (C9) which 
received its moisture from the same general area 
in wmch C8 received its low-level molst_re. 

The temporal relations.-.!? of these events suggests 
a causal relationship, cha specifics of wmen 
cannot be discerrsa from the analysis of reflec- 
tivity data alone 

The relative positions of CS and C9 are 
shown in tne 1740 MDT RHI m Fig 9, By 1751 CS 
and the northern portion of C9 appear to have 
merged Shortly after the C8-C9 merger the 
cells could not be distinguished from one another. 

C8 interacted with two otner subsystems, 

C8a and CIO (see Fig. 2b for locations). C3a 
first became prominent at 1757 MDT wnen it 
appeared as a secondary cell attached to the 
nor tneas term quadrant of C3. C3a intansifit-C 
rapidly to become a separate 47-dBZ ecno entity 
at 1804 !JDT (see Fig. 2b) and, for the next 50 
minutes, travelled due north with a mean speed 

of 7.5 ms 1 . The difference in tne rean prooaga- 
tional velocities of C3 and C8a (150' , /3-4 ms 1 

and 175V7.5 as \ respectively) is quite 
striking The soeed differences, at least, are 
prooabiy due to east— west gradients in the low- 
level momentum. 

Specifically, the northerly component of 
low-level momentum west of the echo line was 
approximately twice as large as mat on the 
eastern side. (See George and Cotton (1978) 

Fig. 13.) C8, whicn lay along the line’s 

westarn border, aopeared to be feeding uoon the 
relatively strong low-level moist nortnerlies, 
thus acquiring northerly momentum. C8a, on tne 
other hand, aopeared to be receiving inflow from 
the region east of the ecno line where tne low- 
level northerly momentum was consideraoly weaker. 
This cell, therefore, would have acquired less 
northerly momentum and consequently exhiDited a 
faster southerly speed. Factors concrioucmg to 
the different direction of propagation of CS and 
C8a have not been resolved. 



Dvuw ( ** M |< J 0-1 


Figure 9. As in Fig. 7, except for C3 and C9 
at 1740 'IDT. C8 is or. t he left, centered az 
44 far. 
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4.0 SUMMARY A HD DISCUSSION 

4.1 Comparison of Results With Other Case 

Studies 

C2’s general echo structure, once establish- 
ed, appeared non-steacy-stace, with the nature 
stage containing 2-3 reflectivity sa'uta. One 
maximum was located along Che storm’ s south- 
eastern region walls the ot-er ref leccivity 
maximum was present vitnin tne ecno overhang 
along the storm's northwestern segment. 

Some of the echo characteristics that C2 
exhibited are similar to the multi-cell hail- 
storms examined by ilarwit: (2972b) ana Cnisboim 
(13o7). In ms study, Cnisnoin founa that storm 
"families" (a family is synomonous with the term 
convective subsystem used here) wera composed of 
several smaller hign-refleccivity cellular 
ecnos. Each cell formed in a preferred area of 
the storm family and moved through the storm 
family in a direction different from tne overall 
movement of the family. Althougn the individual 
cells of C2 could not be accurate!.’ cracked, 
they were multiple and appeared to form m a 
preferred location along the storm's western or 
northwestern side. The appearance of multiple 
cores in C2 implies a storm movement to the left 
of the mean cloud layer winds by discrete propa- 
gatu- of individual cells. 

The echo structure of CU, in contrast to 
that of 02, aopearad to be ouasi-steady-state 
and contained, for the most part, a single 
reflectivity maximum, an associated weax-ecsio 
region, and a nook-like ecno. These features 
are similar to those of suoercell storms described 
by Browning (1964) ana Marwicz (1972a). However, 
several lnporcant differences exist Detween, the 
classical suoercell ana Cll. Most case studies 
of suoercells nave found that tne storm usually 
propagates contmously to tha right of the mean 
environmental winds. These winds usually exnlbi t 
appreciable vertical shear and veering with 
height, especially in the Sub-cloud layer. Cell 
11 prooagaced continuously 30-40° to the left of 
mean environmental (in— cloud) wanes wmch con- 
tained small vertical wind shear and whicn 
backea with height The relatively shallow 
northerly winds did induce appreciaole wind 
shear (nearly a 180° change in wind direction) 
in the suocioud layers. Cll had -ay or inflow 
from formerly low-level winds along its north- 
western quadrant and only minor inflow in the 
vicinity of tne hook-like echo. Typical super- 
cell storms generally have suDstantial inflow 
from low-level southerly winds within the weak 
echo region along their southern cuadranc. 
Apparantly, interactions between Cll anc tne 
mesoscaie were more important m determining 
Cll's echo characteristics chan the dynamics of 
wind snear present in the usual supercell environ- 
ment. 

4 2 Interactions Among Convective Subsystems 

It appears that interactions a m ong convective 
cells of -men the line was composed, and inter- 
actions between convecci.e cells and the local 
mesoscaie circulations, are important m determin- 
ing cell Characteristics. Marvin (1972a), 
using his data and case studies oy other 


investigators' concluded that a”ligut wind 
condition in the suocioud layer is a prominent 
characteristic of the nuiti-ceil storm environment 
Light wines in the sudcIouo layer were part of 
C2's environment . During C2's growth and 
development period, surface winds were generally 

5-10 ms 1 m the subcloud layer. Conversely, a 
change in the mesoscaie and synoptic-scale 
features of the subcloud environment surrounding 
Cll (.a auasi-supercell) produced strorger 

10 ns northerly winds at the’ surface and 
greater subcloud wina shear- 7-' 

Observed interactions among cells included 
echo mergers ana outflow kinematics Outflow 
from convective suosystems appeared to be im- 
portant in new ceil production- Initiation of 
new cells apparantly resulted from the outflow 
of older cells, setting up zones of surface 
convergence on a core localizec scale chan the 
mesoscaie. A good example is Cll. It appeared 
that intensification of tms suosystem was 
enhanced by outflow from an adjacent subsystem 
(CIO) Cell mergers appeared to ole; an important 
role in tne intensification and sustenance of 
the July 19 convective subsystems For example, 
the intensification of C2 appeared to be the 
result of several cell mergers. Florica studies 
(Uoodley and Sax, 1976; Cunning, et al. , 1977) 
have indicated chat cumulus merger everts may 
be important in increasing convergence and cell 
inflow at lower levels slitn the sceeo vertical 
moisture gradients oresent on 19 July, the 
effect of low-level inflow organization and 
intensification would appear to be important in 
drawing a significant increase of water vapor 
into the convective suosystsm. This increase of 
— water vapor influx may then be sufficient for 

the liouid-water-deoeadent process or precioication 
formation ana hailstone growth. 

On several occasions the dissipation of 
weak to moderate cells adjacent to intensifying 
or intense convective suosystems was noted. For 
example, between 153S and 1552 MDT, a moderately 
weak con -eccive cell 10 km south of C2 rapidly 
dissipated while C2 rapidly intensified. It was 
also noted that cells 12 and 13, whicn formed 
adjacent to Cll, exhibited relatively short ecno 
lifetimes and echo intensities compared to the 
other subsystems observed on 19 July Radar 
observations in central Iowa of small convective 
showers m tne vicinity of more dominating 
thunderstorms nave Indicated a similar behavior 
(H C Vaughan, personal communication) . Such 
observations suggest that cne circulations 
induced by a vigorous convective s' T ste’-’ apparently 
are sufficient ac times to mnibic significant 
convection in adjacent regions 

5 0 CONCLUDING REMARKS 

The cells and convective suosystems occurring 
on 19 July 1977 contained unusuallv nigh echo 
intensities and echo tops. The ecno longevity 
and nigh echo intensities are attributed to tne 
presence of che organized mesoscaie low-level 
convergence zone discussed oy Ceorge and Cocton 
(1976) Anothar factor in cell longevity and 
intensity may have been the bortn— south orienta- 
tion of the convective mesoscaie line which was 
parallel to the mean m-clouc southerly winds 


This affect "steered" most cells so that they 
travelled along the echo line m a near- saturated 
environment modified by previous convection. 
Erosion and entrainment effects would nave been 
minimized by this type of geometry Also air at 
ridge coo levels and aoove was quite moist, 
again minimizing entrainment effects. There 
appears to be a feedback mechanism between the 
mesoscale circulations ans the inferred circula- 
tions of the convective subsystems. The structure 
of individual cells or subsystems appeared to be 
dependent upon the stage of development of Che 
mesoscale system, and the location of the convec- 
tive entity relative Co the mesoscale convective 
system. Cll's echo structure and inferred 
circulations appeared to be more organized than 
those characteristics from the other convective 
subsystems. This organization was apparently 
closely related to the organized features of tne 
assoicaced low-level mesoscale convergence field 
that was observed at that time. 

It is nopea tnat suosequent analysis of 
storm-scale motions derived from criple-Dopoler 
radar will shed further lignt on the nature of 
storr — scale and mesoscale interactions. 
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1.0 INTRODUCTION 

Since the first launch of meteorological satellites, researchers have 
been utilizing the resultant unique data set in various types of convec- 
tive cloud studies. For example, Griffith and Woodley (1973) and Reynolds 
and Vonder Haar (1973) have related cloud visible brightness to radar 
cloud top heights to develop relationships for remote rainfall measure- 
ment. Grant and Vonder Haar (1975) have used satellite data to determine 
when favorable seeding conditions exist (using IR sensed cloud top tempera- 
tures) for orographic precipitation enhancement in the western United 
States. The detection, study, and prediction of severe convective 
storms (thunderstorms which produce large hail and/or damaging surface 
winds and/or tornadoes) is an area of high interest since these types of 
damaging storms occur frequently over much of the eastern two thirds of 
the country. The importance of real time detection of large convective 
storms which are producing abnormally heavy rainfalls was tragically 
emphasized by the Big Thompson flash flood in Colorado this past summer. 

Efforts have been made to use satellite data in a qualitative manner 
for severe storm identification, mesoscale analysis, and short range 
"now casting" (see Purdom, 1974; Weiss and Purdom, 1974; and Fujita . 
and Forbes, 1974). Arn (1975) related satellite observed cloud bright- 
ness to reported severe storm occurrences, and Shenk and Curran (1973) 
correlated brightness to observed cloud top heights. Attempts have 
been made to quantitatively relate satellite measured anvil growth rates 
to severe storm occurrences (Arn, 1975) and to radar measured echo charac- 
teristics (Reynolds and Vonder Haar, 1975 and Negri, et.al . 1976). 



3 


Efforts of this type have been limited by the long time period which 
elapses between images (usually 30 minutes at best). Recently Negri, 
et.al . (1976) have emphasized the importance of using short interval 
(5 to 7.5 minute) SMS data in studies of convective clouds and storms. 

A multi-phased study of severe storms and their environments, emphasizing 

i 

the use and interpretation of satellite data, is underway at Colorado 
State University and this paper presents preliminary results which 
compare satellite sensed cloud top radiative temperature fields to 
corresponding (in time and space) radar echo characteristics. 


2.0 DATA 

The SMS-1 satellite collected visible (.5 - .7 pm ) and infrared 
(10.5 - 12.5 pm ) data iat 5 to 15 minute intervals from 1300 GMT, 24 
April to 0200 GMT, 25 April, 1975. Numerous severe storms occurred over 
a region stretching from southwest Oklahoma to eastern Tennessee during 
the afternoon and evening hours this day. A destructive, killer tornado 
struck at Neosho, Missouri and hailstones larger than baseballs pounded 
the small town of Wewoka, Oklahoma. NASA, Marshall Space Flight Center, 
was conducting an atmospheric variability experiment on this date and 
provided CSU with upper air data for the eastern two thirds of the 
country which was taken at three hour intervals. The area of interest 
of this study included eastern Oklahoma, northern Arkansas, and much of 
Missouri. Fig. 1 shows this region, some of the storms considered, 
along with radar indicated heights ( x 10 2 ft. MSI ), and the location 
of the upper air station at Monett, Missouri. Distinctive overshooting 
towers, or tops, are visible on the storms in northeastern Oklahoma. 
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Fujita (1972) and Pearl, et.al . (1975) have related the rapid collapse 
of similar overshooting tops to the occurrence of tornadoes and large 
hail at the surface. The large storm complex in Missouri is intersecting 
a WNW to ESE orientated frontal surface and a line of cumulus clouds is 
visible along this boundary. Qualitative studies by Purdom (1975) have 
shown this type of intersection to be a favored location for the occurrence 
of intense severe storms. However, in this particular case the storms 
in Oklahoma are more intense and dangerous, emphasizing the need for 
quantitative means of analyzing and interpreting satellite data. 

NASA, Goddard Space Flight Center, provided CSU with the SMS-1 VISSR 
(Visible and Infrared Spin Scan Radiometer) digital data, on magnetic 
tape, that were used in this study. Radar data used were the hourly 
reports taken at NWS WSR-57 radar sites at Kansas City, Witchita, Oklahoma 
City, Little Rock, Monett, and St. Louis. Intensity contoured, PPI scope 
photographs (at zero degree antenna elevation) from Oklahoma City were 
compared with the satellite data for the storms in Oklahoma. 


3.0 STORM TOP COMPARISONS 

The semi- isolated nature of the storms of interest made a direct 
comparison with the radar reports fairly simple. Radar indicated echo 
tops at 1935, 0035, and 0135 GMT were compared with nearly simultaneous 
(maximum difference of nine minutes) IR satellite data. If different 
radar stations reported slightly different top heights for the same storm 
the values were averaged. A total of 25 radar heights have been compared 
with the satellite measured radiative cloud top temperatures. The sample 
number will be increased when the digital IR data for the remainder of 
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the afternoon are received and processed - IR data were only available 
for the three times listed above. Smith and Reynolds (1976) reported 
on several different methods of estimating cloud top heights for a 
variety of larger cloud types. They compared top heights obtained using 
cloud shadows, IR-T (Z) relations and radar. Their results showed 

V, 

a fairly large scatter with a root mean difference for an IR/Radar com- 
parison of 2.53 km. Only large thunderstorms which reached near to, or 
penetrated, the tropopause were considered in the current study whereas 
Smith and Reynolds studied clouds which ranged from 6 to 15 km in height. 

A reliable technique to remotely monitor significant storm tops from 
satellite would provide a valuable complement to conventional radar data. 

To determine cloud height using IR data it is usually assumed that 

the cloud emits radiation as a blackbody and that the cloud top effective 

temperature is the same as that of its near environment. Negri, et.al . 

(1976) have detailed some particular problems involved in applying 

these assumptions to obtain top heights for small, isolated clouds in a 

hot, dry environment. However, it is felt that for large storms reaching 

to the tropopause, these assumptions are acceptable. The primary problems 

in monitoring this type of storm result from the coarse resolution of the 

IR data (one IR data pixel at 35N 95W represents the mean blackbody 

2 

temperature for an area approximately 22 km ) and also from the nearly 
isothermal lapse rate above the tropopause. 

Fig. 2 is a plot of the 2100 GMT special upper air sounding taken 
at Monett, Missouri. This temperature profile was considered to be 
representative of that of the environment of all 25 storms used in the 
study. The sounding is potentially very unstable with a lifted index of 
-6 and is similar to what is considered the classic central U.S. tornado 
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proximity sounding (see Miller, 1972). The coldest environmental temp- 
erature present was 209° K (-64°C). 

An overshooting tower rapidly becomes much colder than its environ- 
ment as it overshoots its level of zero buoyancy and penetrates above the 
tropopause height (note the rapid divergence of the moist adiabat and 
the environmental temperature above 200 mb in Fig. 2). However, as thunder- 
storm tops boil above the anvil level they are entraining environmental 
air whose temperature changes little with height. This effect masks 
or smoothes out the presence of very cold temperatures within the nega- 
tively buoyant top. Fig. 3 is a scatter diagram of the radar indicated 
cloud top versus the minimum IR temperature measured for that cloud top. 

The dashed line is an estimated smooth fit of the points and the solid 
line is the Monett sounding (plotted using the radar height scale). 

The two curves agree quite well up to the tropopause (approximately 42,000 
ft.) but with the radar tops consistently higher than IR tops derived 
using the environmental sounding. The higher bias of the radar tops is 
likely due to the coarse resolution of the IR data which doesn't "see" 
small, high rising turrets. The IR temperature measured for the large, 
severe storm anvils is indicated on the figure. The anvils were at or 
very near the tropopause. 

The scatter for tops above the tropopause illustrates the problems 
inherent in inferring cloud top heights for severe storms using only 
satellite IR data. More quantitative techniques are needed to relate 
IR top characteristics to storm severity. Adler (1976) has studied' 
areal rates of change for cold cloud top temperatures during the life of 
the Omaha tornadic storm. His technique will be considered as this study 
continues at CSU, as will other severe storm characteristics. Some 
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possibilities include: examination of the time rate of change of the 

minimum cloud top temperature, development of relationships between extent 
of penetration above the tropopause and the difference between the IR 
temperature and that of the tropopause, and development of relationships 
between both IR and visual characteristics of the storm top to..* its 
severity. 


4.0 PPI DISPLAYS AND IR TEMPERATURE FIELDS 

Figs. 4 and 5 present the Oklahoma City WSR-57 intensity contoured 
PPI display and the nearly simultaneous cloud top temperature field, 
(contoured in 4°C intervals with -70°C the coldest region). The very 
cold cloud top temperatures correlate well with the radar echo configura- 
tions of the large storm complexes at 070°/100 n.m. and 130°/75 n.m. 

2 

The intense echo cells (cross-sectional areas of approximately 50-100 km ) 

detected by the radar within these complexes do not show up within the 

general IR cold cloud area; however, good agreement was noted between 

intense echo cells and overshooting tops apparent in the visible data. 

The problem is again most likely one of resolution and jiornts out the need 

for detailed studies of simultaneous radar and satellite data. Pearl, 

et.al . (1975) studied overshooting tops on a number of severe hailstorms. 

The large overshooting domes had diameters of 10-12 km, or cross-sectional 

areas of approximately 80 km . Only a small part of the area is at the 

top of the dome and an IR data pixel would have to be centered exactly 

2 

on the dome to sense the blackbody temperature of the highest 22 km , 

An interesting feature of this situation is the intense echo SSW 
of the radar. This storm was dropping large hail at the time of the 
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photograph, but it does not stand out in the corresponding satellite IR 
data and would pose a serious problem if it were being monitored only 
by satellite. The radar and satellite evolution of this “hidden cell 11 
is being studied in detail. 


5.0 SUMMARY 

This paper has presented some preliminary results and interesting 
features noted during CSU studies of severe convective storms emphasizing 
the use of concurrent high space and time resolution radar and satellite 
data. Comprehensive, detailed studies of the satellite characteristics 
of severe storms are required if the best use is to be made of the tremen- 
dous amount of satellite sensed data which are now available to both the 
operational meteorologist and the researcher. Methods for quantitative, 
optimum analysis of these complementary data sets should be developed. 
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Figure 1. SMS 1 visible photograph of severe storms studied. Radar 

indicated cloud top heights are shown, as is location of radio- 
sonde station at Monett, Missouri. Time of photograph is 2242 
GMT, 24 April 1975. 




Figure 2. Monett, Missouri sounding at 2100 GMT, 24 April 1975. 

Of interest are the strong potential instability and the 
rapid rate at which a parcel following a moist adiabat would 
become colder than the environment above 200 mb. 



Figure 3. Scatter diagram of radar indicated top height versus satellite 
sensed coldest IR cloud top temperature. The 2100 GMT Monett 
sounding is plotted showing environmental temperature as a func- 
tion of height (for example at 30 Q X 10^ ft. MSL the environmental 
temperature was approximately -40°C). 



Figure 4. Oklahoma City (OKC) WSR-57 intensity contoured PPI scope picture 
at 0122 GMT, 25 April 1975. 
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Figure 5. 


SMS 1 contoured 
25, April 1975. 


IR black body temperature field at 0120 GMT, 
Location of OKC WSR-57 radar is indicated. 
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1 . INTRODUCTION 

High-resolution GOES imagery has been 
increasingly incorporated in studies of severe 
convective storms. Most uses of this satellite 
data have been of a qualitative nature involving 
image interpretation. (For examples refer to 
Maddox, 1977; Purdom, 1976; and Weiss and Purdom, 
1974.) However, some researchers have begun to 
use quantitative, digital GOES data in their 
work. Adler and Fenn (1976) computed growth 
rates of IR measured cold cloud top area for 
several severe storms on 6 May 1975 (one of 
which spawned the Omaha tornado) and compared 
the changes in growth rate with the timing of 
severe weather events at the surface. Negri 
et al. (1976) contrasted visible cloud bright- 
ness changes with the evolution of radar echoes 
for several Great Plains' thunderstorms. They 
also showed that computed parameters, such as 
cloud top divergence, may be affected by the 
time interval separating data used in the compu- 
tations. Maddox et al. (1976) compared IR in- 
ferred storm top heights with radar measured top 
heights and showed that large scatter existed 
between radar and satellite heights for severe 
storms which penetrated the tropopause. 

Houghton and Wilson (1975) presented veloc ity 
divergence and vertical motion fields computed 
for a severe convective situation. The winds 
were derived from tracking clouds on satellite 
photographs. The resultant fields exhibited 
horizontal and vertical coherency suggesting 
that such applications might eventually be of 
operational utility. During the current study 
cumulus clouds were tracked to develop low-level 
wind fields prior to an outbreak of severe thun- 
derstorms. Divergence, moisture convergence, and 
vorticity fields were then objectively computed. 
These analyses were compared to similar fields 
computed using only surface data. Subjective 
modifications of the data set were made to dem- 
onstrate the sensitivity of the results to vary- 
ing spatial distributions of trackable clouds. 

The day chosen for study was 24 April 1975 
when several severe thunderstorms developed 
during late afternoon from north Texas to 
Missouri. GOES 1 imagery was taken at 5-min 
intervals beginning at 1200 CST (all times are 

‘Also at NOAA-ERL, Boulder, Colorado 


given in Central Standard Time) and continuing 
for 8 h. A NASA Atmospheric Variability Experi- 
ment (AVE IV) was conducted this day and upper- 
air soundings were taken at 3-h intervals over a 
large region east of the Rocky Mountains. (The 
NASA AVE program has been described in detail by 
Hill and Turner (1977).) An unusually comprehen- 
sive set of data was therefore available for use 
in the study of the severe storms. Significant 
storm reports included: hailstones with a 

diameter of 2.75 in. (7.0 cm) at Wewoka, Oklahoma 
at 1715; a tornado near Miami, Oklahoma, at 1800; 
and a destructive tornado at Neosho, Missouri, at 
1840. 

2. SYNOPTIC SITUATION 

At 1500 (Fig. 1) a slow moving cold front 
stretched from southern Illinois westward to 
eastern Kansas and then southwestward to the 
Texas panhandle. South of the front a narrow band 



Fig. 1. Surface analysis at 1500 CST. Note the 
dry line which extends from northern Oklahoma to 
southwestern Texas. Winds are in kt with full 
barb = 10 kt. 
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of hot, dry air had advanced eastward through 
north— central Oklahoma into a low pressure area 
east of Ponca City. The region of lowest pres- 
sures was located on the dry line rather than on 
the frontal wave. The surface pattern resembles 
that of the sub-synoptic scale surface low docu- 
mented by Tegtmeier (1974) . He found that tor- 
nadic storms were most likely to occur within the 
northeast quadrant of this type low, and this 
particular case was not an exception. The rem- 
nants of a surface boundary, which separated hot, 
moist air from a region of cool air produced by 
nocturnal thunderstorms, is depicted as a dis- 
sipating warm front. Note that in the narrow 
warm sector temperatures had climbed into the 
upper 80 's with dewpoint temperatures holding in 
the uppers 60' s. 

The 500 mb analysis (Fig. 2) showed that 
a weak short-wave trough extended from south 
Kansas into central Texas. The flow field split 
into two branches to the west over the Great 
Basin, and this short wave was imbedded within 
the southerly stream. A stronger short wave, 
within the northern stream, was moving east- 
southeastward across the northern Plains. By 
0000 on the 25th these two features had moved 
into phase and merged just west of the central 
Mississippi Valley. Although severe thunder- 
storms continued through the night, the most 
intense storms were those that occurred along 
the southern short wave during the late after- 
noon and evening. 
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Skew T/Log P plots of upper-air soundings taken 
at Monett, Missouri, and Amarillo, Texas, are 
shown in Fig. 3. (Release times were 1500 and 
1415 respectively.) The Monett plot resembles a 
Miller (1972) Type I tornado sounding. The low- 
est 100 mb layer L.I. was -4, and the mean vapor 
mixing ratio for the same layer was almost 12 g 
kg - l. However, note that the LFC was almost 
200 mb above the LCL with a large negative buoy- 
ancy area indicated through the separating layer. 
Most veering in wind direction occurred below 
the LCL with west-southwesterly flow extending 
from 800 through 100 mb. The temperature lapse 
rate was strongly superadiabatic in the lowest 
30 mb. 


The Amarillo sonde was released within the 
hot, dry air mass prior to cold frontal passage. 
A well mixed, nearly adiabatic layer extended 
upward from the surface to 600 mb. Although 
strong westerly flow was present at middle and 
upper levels, nearlv uniform westerly winds of 
only 10 kt (Sms ”1) existed in the mixed layer. 
Because of this weak flow the dry line moved 
only slightly eastward during the afternoon, 
except over central Oklahoma, where it moved 
eastward toward the low pressure region. 



Fig. 2. 1500 CST 500 mb analysis. Triangles 
show locations of short-wave troughs. Arrows 
depict axes of maximum wind speeds, full barb = 
10 kt, flag = 50 kt. Moist regions with T-Td 
6°C are shaded. 


Fig. 3a. A Skew T/Log P plot of 1500 CST Monett, 
Missouri, upper-air sounding. 


55 


PRESSURE 



Fig. 3b. Skew T/Log P plot of 1500 CST Amarillo, 
Texas, upper-air sounding. 

Wilson (1976) has computed large scale 
vertical motion fields (700 mb u computed with an 
adjusted kinematic method) using the AVE IV 
sounding data. His 1500 700 mb w field is shown 
in Fig. 4, along with the concurrent surface 
analysis. Regions where severe thunderstorms 
developed during the following 2 h are 
indicated. The large scale vertical notion 
fields apparently did not play a dominant role 
in organizing convection on this particular day 
since the storms developed in regions of 
diagnosed subsidence, or of only slight upward 
motion. 

The 1702 satellite photograph is shown in 
Fig. 5 with the 1700 surface analysis super- 
imposed. The three large storms in northern 
Missouri, northeastern Oklahoma and southern 
Oklahoma were all severe. The Missouri storm- 
complex was oriented north to south and inter- 
sected the surface front and an associated east- 
to-west line of convective clouds. Empirical 
rules (Purdom, 1976) suggest that this storm 
would be a likely tornado produces Similarly, the 
large storm in southern Oklahoma would be con- 
sidered highly suspect because of its extremely 
rapid growth rate (refer to Negri et al. (1976) 
for details), and the obvious flanking lines 
which extended west-southwestward from the storm 
complex. However, it was the less spectacular 
storm (as viewed from satellite) located east- 
northeast of the sub-synoptic surface low that 
became tornadic. 
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Fig. 4. 1500 CST 700 mb vertical velocity in 

tibars s _1 (from Wilson, 1976). Regions where 
severe thunderstorms developed during the 
following 2 h are cross-hatched. 




Fig. 5. 1702 CST COES photograph with 1700 

suface analysis. 


3. SATELLITE DATA SET 


A low-level wind field was developed by 
tracking cumulus clouds through a four-picture 
loop which began at 1458 and continued through 
1513. The tracking was accomplished on the NASA 
Atmospheric and Oceanographic Information Pro- 
cessing System (AIOPS) . Manual cursor position- 
ing methods were employed. (The AOIPS and cloud 
tracking techniques are described by Billingsley 
et al. (1976).) Acceptance of a derived vector 
was subjective and was primarily based- upon 
obtaining consistent cloud motion through the 
loop. 
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The set of 214 wind vectors obtained 
is shown in Fig. 6. These data were oDjectivelv 
interpolated onto a 0.4 degree grid (approx- 
imately 9 degrees square and centered on south- 
western Missouri) using a 2-D cubic spline 
technique which was developed by Fritsch (1971) . 
The resultant wind field is shown in Fig. 7. 

The spatial distribution of the data points 
indicated that the validity of the gridded winds 
would be most questionable over northwestern and 
southeastern corners of the domain, regions not 
affected by storms. 

Results of previous studies have indicated 
that cumulus motion represents the cloud-base 
environmental wind to within 5 degrees and 
1 m s -i (Fujita et al. , 1975 and Hasler et al. 
1975). Hasler et al. also noted that there 
were no discernable biases to the differences 
indicating that analysis of a large, closely 
spaced data sample might minimize the error 
present. 

Several dynamic fields were computed using 
the wind field shown in Fig. 7. The represen- 
tativeness of these fields is partially a 
function of the spatial distribution of track- 
able cumulus clouds, and this important 



Fig. 6. Set of 214 wind vectors obtained by 
tracking low-level clouds on the NASA AOIPS. 
One degree of latitute = 22.2 m s - * and one 
degree of longitude = 17.8 m s“l. 
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Fig. 7. Objectively gridded field of wind vectors. 
Scaling is same as in Fig. 6. 


limitation is considered in a subsequent sectio- 
Location of the wind vectors in z-space is also 
a significant problem. Schaefer (1973) noted 
that usually neglected terms which appear when 
divergence is computed on a nonhorizontal surface 
can significantly affect the resultant fields, 
.lowever, in this case it was assumed that the 
satellite derived winds were on a horizontal 
surface, even though rawinsonde data indicated 
that the LCL ranged from 910 to 800 r across 
the region where clouds were trackeu. Al...ougi: 
simplifying assumptions have been made, the 
results do seem promising. 

4. MOISTURE CONVERGENCE FIELDS 

Hudson (1971), Newman (1972), and Sasaki 
(1973) have all shown that severe thunderstorms 
and squall lines tend to develop in regions 
where values of moisture convergence are large. 
Identifiable maxima in the surface moisture 
convergence field may precede development of 
intense convection by as much as 2 to 3 h. 

Studies by Pearson et al. (1967), Steyaert 
and Darkow (1973), and Schaefer (1975) have 
presented results indicating that surface mixing 
ratios might be combined with the satellite de- 
rived low-level wind field to approximate mean 
mixed layer moisture convergence. Surface mixing 
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ratios were computed for the 1500 observations 
using Tetens' empirical formula for vapor 
pressure and the definition of mixing ratio. 

The mixing ratios were objectively analyzed onto 
the same grid as shown m the wind field of 
Fig. 7. These two fields were used to compute 
the moisture divergence (V.wV) , hence moisture 
convergence is of negative sign. 

The resultant field is presented m Fig. 8. 
Comparison with Fig. 5 shows chat during the 
following 2 h storms developed in most of the 
regions characterized by large values of mois- 
ture convergence at 1500. The mam exception 
was the lack of activity within the moisture 
convergence zone in central Kansas. The mois- 
ture convergence field calculated using only 
surface data is shown in Fig. 9. The overall 
pattern of the two fields is similar; however, 
significant differences do exist. The most 
pronounced of these was the band of high mois- 
ture convergence values that stretched across 
south-central Oklahoma in the satellite field. 
The surface analysis indicated moisture 
divergence over much of this area. This was an 
important difference since several severe 
storms developed across this region. 


-oo -08 -Q7 -96 -95 -94 -95 -92 -91 -90 


Fig. 9. 1500 CST moisture convergence field 

computed using observed surface winds Details 
are same as in Fig. 8. 

Maximum values for satellite and surface 
moisture convergence (~8.5 g kg” 1 h -1 ) were 
similar to values found by Newman (1972) to be 
associated with severe thunderstorm occurrences. 
The Neosho tornadic storm began to develop within 
an hour of the satellite measurements and within 
the area of highest moisture convergence values. 
An improved correspondence between maxima and 
subsequent severe storm genesis were character- 
istic of the satellite field. This may have been 
a result of both stronger flow fields at cloud 
base and an elimination of local influences which 
affect measurements of surface winds. 

5. DATA SENSITIVITY 

A modification of the original 214- wind 
vector set was made to demonstrate that the 
spatial distribution of trackable cumulus clouds 
strongly affects the results of calculations of 
the type considered above. The 24 April tornado 
case was unusual because strong synoptic scale 
weather systems and attendant large cloud masses 
were not present. Eleven of the wind vectors 
were obtained by tracking individual cumulus 
within a narrow cloud band that developed to the 
rear of the surface cold front. These vectors 
were arbitrarily deleted from the data set. 
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Fig 8. 1500 CST moisture convergence field 

computed using satellite derived winds. Units 
are g kg- 1 s^ 1 x 10~ 5 and contour interval is 
5 x 10 - ^ g kg" 1 s -1 . 
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Calculated relative vorticity fields are 
shown in Fig. 10 for both the complete and 
modified data sets. The results differ both in 
general pattern and also by as much as two 
orders of magnitude. The data modification was 
chosen to maximize differences; however, cloud 
patterns_ associated with real weather systems are 
also likely to systematically affect the discn- 
bution of trackable clouds, For example, it 
would be most unusual to find a cumulus cloud 
field m the air mass behind a dry line, and 
low-level clouds would often be obscured north 
of a surface' cold front by layers of middle and 
high clouds. 

These limitations emphasize the importance 
of integrating satellite determined low-level 
wind fields into a boundary layer model (both 
diagnostic and prognostic) which incorporates 
all available conventionally and remotely 
sensed data. This is especially true when 
operational applications are considered. 



Fig. 10a. 1500 CST relative vorticity computed 
using satellite derived wind field. Units _are 
s~l x 10“4 an( j contour interval is 5 x 10 


~ 3 a"*. 



Fig. 10b. 1500 CST relative vorticity computed 

using the modified satellite derived wind set. 

6. SUMMARY 

Cumulus clouds were tracked, using short- 
interval high-resolution GOES satellite imagery, 
to derive a low-level wind field for a severe 
thunderstorm situation. The cloud tracking was 
accomplished during a period prior to storm devel- 
opment. Moisture convergence fields were approx- 
imated using these satellite derived winds and 
surface mixing ratios. Subsequent severe storm 
development occurred within regions of diagnosed 
high moisture convergence. Satellite moisture 
convergence patterns specified convective patterns 
that developed much better than did the same field 
computed using only surface observations. These 
data suggest that useful applications, using 
satellite wind fields, may be possible m opera- 
tional severe thunderstorm forecasting. 

It was also shown that computed dynamic 
fields are highly affected by the spatial dis- 
tribution of trackable clouds. Satellite winds 
must be accurately located in the vertical before 
more precise computations can be obtained. Work 
is underway on development of an optimum objec- 
tive analysis scheme which would simultaneously 
analyze satellite winds, surface data, rawmsonde 
data, and other types of remotely sensed data. 
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1. INTRODUCTION 

The launch of geosynchronous satellites has 
provided meteorologists with high spatial and 
temporal resolution data essential in evaluating 
mesoscale features which may force the develop- 
ment of severe local storms (SELS) . At the 
present time, only radiances from two channels, 
the 11-12 pm infrared and the 5-7 pm visible are 
available, so that vertical temperature profiles 
cannot be obtained. However, Vertical Tempera- 
ture Profile Radiometer (VTPR) data are avail- 
able from the polar orbiting NOAA-4 satellite, 
which uses the 15 pm CO2 channels for tempera- 
ture retrieval and the rotational water vapor 
absorption band to infer ptecipitaole water 
values (Hillger and Vonder Haar, 1976) , however 
temporal resolution is poor. This paper explores 
the possibility of combining these two data sets 
to simulate data that will become available with 
the launch of the VISSR-*- AtmosDheric Sounder 
(VAS) on the GOES-D in the 1980' s. A moisture 
analysis of 70 km resolution is combined with 
wind fields derived from tracking low-level cu- 
mulus clouds using SMS visible channel data. 

This combination allows for the computation of 
moisture convergence in the boundary layer for 
the pre-storm environment. 

Several researchers, including Hudson (1971), 

Newman (1972) and Sasaki (1973) computed surface 

moisture convergence and found good agreement 
between those patterns and reports of severe 
weather several hours later. Intense moisture 
convergence (~ 10 - 3 g kg~l s~F) is one condition 
found favorable for tornadic storm development 
(Sasaki, 1975). 

The day chosen as a case study was 24 April 
1975, when several severe thunderstorms developed 
during the late afternoon from north Texas to 
Missouri. 

This day is well represented in the litera- 
ture. Whitney (1977) examined the role of the 
subtropical and polar jets on this day and con- 
cluded that upper-level diffluence ahead of the 
advancing trough enhanced the convective insta- 
bility north of the sub-tropical jet. Maddox 
et al (1976) compared satellite sensed radiative 
temperature fields to corresponding radar echo 
characteristics in an attempt to quantify the 
detection of severe storms by satellite. Negri 
et al (1976) examined growth rates m SMS visi- 
ble channel data to determine an operationally 
significant parameter for detection of storm 
severity Purdom (1975) examined the low-level 
inflow into these storms from both an earth 
relative frame and a thunderstorm relative 
frame and concluded that a relative low-level 
southerly flow into a storm is important 
in that storm’s ability to produce a 
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tornado This paper will follow the latter ap- 
proach in using low-level cumulus velocities 
to compute divergence and subsequently moisture 
divergence. 

2 24 APRIL 1975 

2 1 Synoptic Situation 

Numerous severe storms occurred m a region 
stretching from southwest Oklahoma to eastern 
Tennessee during the afternoon and evening hours 
this day. A destructive, killer tornado struck 
at Neosho, Mo. and baseball-sized hail pounded 
the small town of Wewoka, Ok. Figure 1 is the 
surface analysis for 1600 GMT, near the time of 
the N0AA-4 satellite pass Features to note are* 

1") the intense moisture gradient across the 
dry line in western Oklahoma 

2) very moist air in eastern parts of Texas 
and Oklahoma 

3) cold frontal boundary m the Oklahoma 
panhandle, becoming stationary out Cowards 
the northeast 

4) squall line from previous days convective 
activity running through Oklahoma and 
Arkansas . 



Fig. 1. Surface analysis for 1600 GMT, 24 April 
1975. 


'''Visible and Infrared Spin-Scan Radiometer 
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Figure 2 is the 1500 GMT 500 mb analysis of data 
collected during NASA's Atmospheric Variability 
Experiment (AVE) on this date. The weak short 
wave moving through the midwest increased the 
instability over the region 


radiances for the CO 2 channels reached some limit- 
ing noise value (the instrumental noise level for 
the VTPR instrument) . 


.HOflA VTPR CHANNELS 



Fig 2. 500 mb analysis for 1500 GMT. Data 

were collected during NASA's Atmospheric 
Variability (AVE) experiment. 

2.2 Mesoscale Precipitable Water Analysis 

High horizontal resolution moisture infor- 
mation was obtained at 1615 GMT on 24 April 
1975. This moisture analysis was available in 
terms of total precipitable water (PW) values 
at a resolution of about 70 km from the VTPR 
on NOAA— 4. Table 1 shows the spectral charac- 
teristics of Che VTPR instrument. Clear column 
temperature profiles were retrieved through use 
of an iterative temperature retrieval algorithm, 
(Hillger and Vonder Haar, 1977). As an initial 
guess profile the retrieval program used a com- 
posite temperature and moisture sounding con- 
structed from 1800 GMT AVE radiosonde profiles. 
The temperature profile was then iterated until 
the RMS residual between observed and calculated 


APPROXIMATE 


CHANNEL 

HUMBER 

WAVELENGTH 
(micron? ) 

WAVENUMBER 

U"' 1 ) 

PEAK LEVEL OF 
WEIGHTING FUNCTION 
(millibars) 

C0 2 channels {15 micron CO^ absorption band) 

i 

14.96 

668 5 

30 

EQ branch) 
2 

14 77 

677 5 

50 

3 

14 38 

695 0 

120 

4 

14 12 

708 0 

400 

S 

13 79 

725 0 

600 

6 

13 38 

747 0 

surface 

HjO channel 

(Rotational water vapor absorption band) 

7 

18 69 

535 0 

700 

Window channel 

(Atmospheric window region) 

8 

11.97 

833 0 

surface 


Table 1 


Table 1. Spectral characteristics of the Vertical 
Temperature Profile Radiometer. 


After the temperature profiles were retrieved 
using the six CC >2 channel radiances, a comparison 
was then made between the observed and calculated 
radiances for the VTPR H,0 channel The observed 
t^O channel radiance is dependent upon both the 
temperature and moisture structure of the sound- 
ing column, whereas the CO 2 channels are in gen- 
eral less dependent on moisture than temperature. 

The calculated radiance is an integrated value 
obtained by applying the radiative transfer equa- 
tion to the retrieved temperature profile with 
the initial guess PW amount. If the retrieved 
temperature profile is sufficiently accurace, then 
the only difference between observed and calculated 
radiances should be due to the difference between 
the observed and initial guess precipitable water 
amounts. The initial guess PW amount is not 
changed in the iterative process and does cause 
Some errors in the derived temperature profile 
which are of smaller magnitude than the PW effect 
on the H 2 O channel. Since the initial guess PW 
amount is a constant value, the H^O radiance re- 
sidual should therefore be proportional to the 
actual PW amount in the sounding column. 

For the area of the satellite pass over the 
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Great Plains, a correlation between the H-0 radi- 
ance residuals and the actual PW values a£ AVE 
radiosonde launch sites is shown in Figure 3. 

The rather high correlation of -.88 shows that 
much of the moisture distribution on the raeso- 
scale was explained by the radiance residuals. 

All the values except Fort Sill, Oklahoma (FSI) 
are within the absolute error bias of +0.5 cm of 
HO which is considered as a maximum error for 
this data set. 



Fig. 3. Observed precipitable water values versus 
H^O channel radiance residuals. 

A least squares linear fit was made to the indi- 
vidual points to derive PW values at the places 
other than radiosonde sites which were used for 
calibration. The resultant PW field is shown in 
Figure 4 for the satellite pass time. One region 
in question is in central Kansas where high clouds 
caused the PW values to be low. The clouds there 
extended above 500 mb and the values obtained are 
not reliable. Otherwise, the individual values 
should be accurate to within + 0.5 cm of PW. 



Fig. 4. The precipitable water field (cm) at 
1615 GMT. 


2.3 Wind Fields 


Five minute interval SMS visible channel data 
were used to track cumulus clouds and derive low- 
level wind fields from their motion. Winds were 
tracked on NASA’s Atmospheric and Oceanographic 
Information Processing System (AOIPS) by NASA 
scientists and the data furnished to CSU. Cloud 
location is determined by means of a cursor, which 
is moved about the screen by means of a joystick. 
Prior navigation of the images to within one pic- 
ture element (- 1 km) allows for cloud tracking 
in an earth relative frame, [a description of the 
AOIPS and cloud tracking techniques was presented 
by Billingsley (1976).] Figure 5 is one such 
satellite image used. Cumulus clouds were tracked 
in an area approximately 11 degrees on a side 
centered near Oklahoma City, Oklahoma. Figure 5 
shows the original wind field, with a vector of 
length one degree representing a wind speed of 
17.8 ms -1 u velocity and 22.2 ms-1 v velocity. 
Hasler , et.al. (1976) showed that for tropical 
cumulus, excellent agreement existed between cloud 
motion and the wind at cloud base ( v c ^ ou( j " V u i n d 
<_ 1.3 ms - ^) . 

Fujlta et.al, (1975) used ATS Images to infer 
environmental winds and concluded that cumulus 
turrets 0.3 to 2 miles In size appear to be the 
best targets to infer mean winds in the sub-cloud 
layer. Accuracy in velocity computation was shown 
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to he t 1 ns speed and + 5° direction. More 
recently, Suchraan and Martin (1976) explored the 
accuracy and representativeness of tracer winds 
in the GATE area. They found that ship winds 
used as "ground" truth differed from satellite 
winds by less than 3 ms _i . Using 30 min%te in- 
terval ^magery, they found a reproducibility of 
1.3 ms for cumulus level winds. Cumulus 
tracked winds would seem to better indicate 
boundary layer motion than do the corresponding 
surface winds because: 

1) a 20 minute "average" of mean flow is 
obtained, and 

2) surface winds are influenced by topog- 
raphy and surface stress. 

Dosvell (1976) recognized the Importance of using 
"filtered" surface data in computations of low- 
level moisture divergence. 


ORIGINAL PAGE IS 
OF POOR QUALITY 



Fig. 5. SMS visible channel image at 1757 GMT. 


3. RESULTS 

3.1 Dynamic Parameters 

The wind field of figure 6 was objectively 
analyzed using a cubic spline technique developed 
by Frltsch (1971). Figure 7 shows the objectively 
analyzed wind field, using a grid interval of 0.5 
degrees, with lengths scaled as in figure 5. The 


divergence field, figure 7, shows Che areas of 
strongest low-level convergence located in cen- 
tral Kansas, southwestern Oklahoma and southern 
Missouri. Magnitudes of this convergence are 
10-20 x 10 -5 s _ L A final computron Is done with 
the addition of the satellite derived PW fields. 
The moisture field of figure 8 was objectively 
analyzed onto the same grid as the winds and the 
moisture divergence (7 • PW V) was calculated. 

This field is illustrated in figure 9. and is 
similiar in structure to the divergence field 
alone. The maximum moisture convergence was 
'30 x 10 - ^ cm H.,0 s -1 centered in eastern Kansas. 
Local maxima are also found in southern Missouri, 
southwestern Oklahoma, central Texas and northern 
Louisiana. The radar summary at 1835 (Fig. 10) 
shows some convective activity ooccurring in 
Kansas, however these storms were not severe. 



Fig. 6. Original wind field at 1800 CMT from 
satellite Inferred cumulus velocities. 


Because of the rather large time discrepancy 
between the moisture field (1600), the wind field 
(1800) and the time of initial storm genesis (2200), 
emphasis has been placed on the description of a 
technique rather than on demonstrating the possible 
correlation between moisture convergence patterns 
and the occurance of severe local storms. Maddox 
et.al. (1977) coupled 2100 GMT derived winds with 
corresponding surface mixing ratios on this day 
to determine possible correlations. The authors 
are well aware of some of the shortcomings of the 
technique described. A satellite wind which was 


noc accurately located in the vertical has been 
multiplied by an esciroition of the total atmo- 
spheric water content. Hie wind fields are sen- 
sitive to the density and spatial configuration 
of the cumulus clouds used to derive them, as 
well as to the objective analysis scheme used 
Total PW concent is likewise noc retrieved per- 
fectly. It is hoped chat by Che time VAS is 
launched there will be improved algoritnms for 
determining H^O content and its vertical dis- 
tribution 
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Tig 7. Objectively analyzed wind field ac 1800 
CUT Orid interval is 0 5 degrees 



tig 10. hWS radar depiction chart for 1835 (SIT. 



Fig. 8. The 1800 GMT divergence field Units 
are 

gence) have been contoured. 


10 - -* s 1 and only negative values (conver- 



Fig 9. The moisture divergence field. Units are 
10 -5 cm 1^0 s-1 Only negative values (moisture 
convergence) have been contoured 
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4 . SUMMARY 

The major results and limitations of this 
preliminary study of combining two types of sat- 
ellite data included: 

1) Using data from satellite sensors only, 
wind and moisture fields can be con- 
structed and then combined to give meso- 
scale horizontal moisture flux infor- 
mation. 

2) Satellite inferred cumulus velocities 
provide a usable representation of 
boundary layer flow, with the limitation 
of not being exactly located in the ver- 
tical Subsequent calculations may be 
dependent on the number and density of 
wind vectors and on the type of objec- 
tive analysis used. 

3) While moisture profiles cannot be re- 
trieved perfectly, total precipitable 
water can be recovered through its large 
negative correlation with radiance re- 
siduals. Most of this PW resides in 
the boundary layer. 

4) Moisture divergence fields reflect im- 
portant mesoscale forcing mechanisms of 
severe local storms, notable in the ab- 
sence of a well defined synoptic 
pattern. 

5) Because of the early NQAA-4 pass time, 
and a 2 h discontinuity between the 
wind and moisture fields, no inferences 
can be made to the correlation of low- 
level moisture convergence and the oc- 
currence of severe local storms for this 
case 

The proposed VISSR Atmospheric Sounder, with 
its high spatial and temporal resolution, will 
provide a vital tool for the implementation and 
evaluation of this proposed technique. 
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